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j^ ^ Abstract. We establish sharp pointwise Green's function bounds and consequent 

linearized stability for smooth traveling front solutions, or relaxation shocks, of gen- 
eral hyperbolic relaxation systems of dissipative type, under the necessary assump- 

C^^ ■ tions ([G,Z.1,Z.4]) of spectral stability, i.e., stable point spectrum of the linearized 

operator about the wave, and hyperbolic stability of the corresponding ideal shock of 
the associated equilibrium system, with no additional assumptions on the structure or 

Qb^ I strength of the shock. Restricting to Lax type shocks, we establish the further result 

■^r ■ of nonlinear stability with respect to small L^ Ci H^ perturbations, with sharp rates 

of decay in L^, 2 < p < cx), for weak shocks of general simultaneously symmetrizable 
systems; for discrete kinetic models, and initial perturbation small in W^'^ n W^'°° , 
we obtain sharp rates of decay in L^, 1 < p < oo, for (Lax type) shocks of arbitrary 
strength. This yields, in particular, nonlinear stability of weak relaxation shocks of 
the discrete kinetic Jin-Xin and Broadwell models, for which spectral stability has 

^h . been established in [HL,JH] and [KM], respectively. 

^ ' Our analysis follows the basic pointwise semigroup approach introduced by Zum- 

r^ . brun and Howard [ZH] for the study of traveling waves of parabolic systems; however, 

significant extensions are required to deal with the nonsectorial generator and more 
singular short-time behavior of the associated (hyperbolic) linearized equations. Our 

(^ ' main technical innovation is a systematic method for refining large-frequency (short- 

time) estimates on the resolvent kernel, suitable in the absence of parabolic smooth- 

^D ' ing. This seems particularly interesting from the viewpoint of general linear theory, 

r^ . replacing the zero-order estimates of existing theory with a series expansion to ar- 

bitrary order. The techniques of this paper should have further application in the 
a closely related case of traveling waves of systems with partial viscosity, for example 

in compressible gas dynamics or MHD. 
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c^ ■ Section 1. Introduction 

A variety of nonequilibrium processes in continuum mechanics can be modeled 
by hyperbolic relaxation systems of general form 

(1.1) f"U(?'1) =(-M , 

u, f E M"-, V, g, q E W , where 

(1.2) Rea{q^{u,v*{u))) <Q 
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2 STABILITY OF RELAXATION SHOCKS 

along a smooth equilibrium manifold defined by 

(1.3) q{u,v*{u)) = 0, 

and r is a (usually small) parameter determining relaxation time: for example, 
non-thermal equilibrium gas dynamics [W,MR,Lev], traffic dynamics [LW,AR,Li], 
and multiphase flow [BGDV,BV,NT,MaP]; for a general discussion of relaxation 
models, we refer the reader to [L.2,CLL] or to the surveys [N,Yo.4]. 

A particularly natural class of relaxation systems are those obtained by discretiz- 
ing continuous kinetic models such as Boltzmann or Vlasov-Poisson equations, for 
which the lefthand side of (1.1) reduces to a simple linear transport operator: for 
example, Broadwell and other lattice gas models [PI]. Here, loosely following Na- 
talini [N], we define the class of discrete kinetic models as systems of form (1.1) 
for which / and g are linear, constant coefficient. These include in particular the 
subclass of semilinear relaxation schemes, such as the Jin-Xin [JX] and Natalini [N] 
models, that have enjoyed considerable recent popularity as a method for numerical 
approximation of hyperbolic systems. 

The relaxation system (1.1) can be viewed "to zeroth order" as a regularization 
of the associated equilibrium, or "relaxed" system of conservation laws 

(1.4) ut + r{u)^ = 0, 

f*{u) := f{uj v*{u)), which system may or may not be hyperbolic (see [JX,BGDV,BV] 
for examples of nonhyperbolicity) . To "first order," it may be approximated, at least 
formally, by an associated parabolic system of conservation laws 

(1.5) ut + r{u)^ = r{B*{u)u^),, 
where 



(1.6) 



B*{u):=-f.q-\g:-v:f:) 

= -fvQv^{9u - 9vqy^qu + Qy^quifu - fvqy^Qv 



g* := g{UjV*{u)), is determined by Chapman-Enskog expansion as described in, 
e.g., [Wh,L.2,Z.l]. Here, the order referred to is with respect to relaxation time r; 
however, as pointed out by Liu in his seminal work [L.2], this formal approximation 
has a more readily justifiable interpretation in terms of large time, rather than small 
r behavior. It is this point of view that is relevant for the present work, in which we 
investigate large time behavior for a fixed relaxation time. From now on, we take 
without loss of generality r = 1, and suppress the parameter r. 

Note that, in the frequently occurring case that r < n, the n x n matrix B* 
is necessarily singular, since it factors through the r x r matrix q~^. Thus, the 
proper analogy is to "real viscosity," or partially parabolic systems of conserva- 
tion laws. Indeed, this analogy holds even in the case that B* is nonsingular. 



CORRADO MASCIA AND KEVIN ZUMBRUN 3 

since the high frequency behavior of the original system (1.1) is by necessity hyper- 
bohc, irregardless of formal expansion; see [Ro, SI. 1-2, Ma] for interesting discussions 
on the validity/meaning of the Chapman-Enskog expansion in the high-frequency 
regime. The analogy is in fact much deeper than this. As pointed out by Zeng 
[Ze.2], the two problems (relaxation/real viscosity) are essentially dual in the linear 
constant-coefficient case. A similar duality, at the structural level, may be seen in 
the variable-coefficient case [Z.4] 

A particularly interesting phenomenon, suggested by (1.4), is the existence of 
smooth traveling front solutions 

(w, v) (x, t) = (w, v) {x — st) 
lirn {u,v) = {u±,v±), 

z—>±oo 

of (1.1), where by necessity v± = v*{u±) and u± corresponds to a shock solution 
of (1.4). See [L.2] for a treatment of existence in the general case n = r = 1, 
and [YoZ,FZe] for generalizations in case n or r > 1; further results/discussion 
are given in Section 1.1, below. Such traveling waves are known as relaxation 
shocks or relaxation profiles. The question of their stability was investigated in 
[L.2] in the general case n = r = 1, for which the equilibrium system is scalar, 
under the assumption of "weak," or small-amplitude shocks. This analysis has 
been generalized to strong shocks in the special case of the 2x2 Jin-Xin model 
[MaN]. For the 3x3 Broadwell model, a standard discrete kinetic model for which 
n = 2, r = 1, Szepessy and Xin [S,SX.2] have announced the result of linear and 
nonlinear stability of weak shocks, extending partial results of [KM,CaL]: to our 
knowledge, the only such result for the system case n > 1 occurring in most physical 
applications. However, to date there is no complete analysis of stability for any 
other case, in particular for the case n > 1, r >> 1 arising through approximation 
of Boltzmann or Vlasov-Poisson equations by moment closure or discretization. For 
strong shocks, there is no complete analysis for any system such that n or r > 1. 

Useful necessary conditions for stability have been obtained for shock profiles of 
the Jin-Xin relaxation model in [Z.1,G] and for shock profiles of general relaxation 
and real viscosity models in [Z.4], generalizing corresponding results of [GZ,BSZ,ZS] 
in the strictly parabolic case. Strengthened versions of the one-dimensional hy- 
perbolic stability criterion of Erpenbeck-Majda [Er,M.l-3], these conditions yield 
interesting results of spectral instability (hence linearized and nonlinear instability) , 
similarly as in the strictly parabolic case (see e.g. [GZ,FreZ,ZS,Z.4]). Notably, the 
above-mentioned analyses apply to shocks of arbitrary strength and type (i.e.. Lax, 
under-, or overcompressive) . 

However, to our knowledge, for n > 1, the only positive stability result for 
relaxation shocks other than the ones of [KM,CaL,S] is a partial result recently 
established by H. Liu [Liu] of linear and nonlinear zero-mass stability of weak shocks 
of the Jin-Xin system, i.e., stability with respect to perturbations whose integral in 
the u variable is zero; see also a subsequent and closely related proof by Humpherys 
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[Hu] of spectral stability of small-amplitude Jin-Xin shocks, defined as stability 
of the point spectrum of the linearized operator about the wave. (Note: spectral 
stability is implied by zero-mass stability, see discussion, [ZH] or [Hu]). Being based 
on energy methods, none of these previous analyses for systems ([KM,S,Liu,Hu]) 
generalize to the case of strong (i.e., large-amplitude) shocks. Moreover, they each 
depend to some degree on special structure of the respective systems treated; hence, 
the extent to which they generalize to other systems is not clear. 

In the present paper, following the general philosophy set out in [ZH] for the 
study of viscous shock waves, we investigate the linearized and nonlinear stability 
of relaxation shocks of general systems (1.1), with arbitrary n, r, but assuming 
spectral stability. As in the parabolic case, what makes this enterprise nontrivial is 
that the essential spectrum of the linear generator accumulates at the imaginary 
axis, so that linearized decay cannot immediately be concluded from spectral sta- 
bility, and moreover can be at most at time-algebraic and not exponential rate. 
Indeed, much of the technical difficulty of the problem lies in this usually trivial 
step. Spectral stability can then be investigated separately as an ODE problem, 
or checked numerically as in [Br. 1-2]. An important advantage of this approach is 
that it is completely general, applying in principle to shocks of all types (i.e.. Lax, 
under-, or overcompressive) and strengths.^ 

More precisely, we study linearized stability for systems (1.1) satisfying appro- 
priate dissipation conditions as in [Yo.4] or [Ze.2], and nonlinear stability for weak 
Lax shocks of (general) symmetrizable systems (defined Section 7.2), or strong Lax 
shocks of discrete kinetic models. Our main example will be the Jin-Xin model, for 
which calculations are particularly transparent. In this case, combining our results 
with the existing spectral stability results of [Liu,Hu], we obtain the result that Lax- 
type Jin-Xin shocks of sufficiently weak strength are nonlinearly stable under small 
L^ n L°° perturbations, with no additional hypotheses: the first such result for a 
relaxation shock in the case n, r > 1. Likewise, using the zero- mass stability result 
of [KM] (recall: zero-mass implies spectral stability), we recover a more detailed 
version of the result of [S,SX.2] that weak Broadwell shocks are nonlinearly stable, 
with the additional information of sharp rates of decay in all L^. (Note: No rates 
of decay are given in [S,SX.2].) In the course of our analysis, we develop extremely 
detailed bounds on the Green's function of the linearized evolution equations about 
the traveling wave, of interest in their own right. These generalize bounds obtained 
in [ZH,Z.2-3] for shock profiles of strictly parabolic systems, and in [Ze.l,LZe], [Ze.2] 
for constant solutions of partially parabolic and relaxation systems, respectively. 

Our results are of interest from the physical point of view as some of the first com- 
plete shock stability results for systems with realistic dissipative mechanism such 
as relaxation or partial viscosity: in particular the first positive stability results for 
strong relaxation shocks (namely, the reduction to the numerically verifiable spec- 
tral stability condition {T>)), the first for quasilinear systems, and the first results 
for relaxation shocks of any strength to give rates of decay. From the technical side. 



^However, see Remark 2 just below, regarding subshocks. 
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they are of interest as a blueprint for the apphcation of the pointwise semigroup 
methods of [ZH] to the case that the associated semigroup is C^ and not analytic: 
in particular, when the Green's function G{x,t;y) is a distribution, lying outside 
the class of smooth functions for t > 0. 

With regard to the latter, technical issue, we mention related analyses carried 
out by Dodd [Do] and Howard-Zumbrun [HZ. 2] for shock profiles of dissipative- 
dispersive conservation laws, for which the semigroups associated with the linearized 
equations are likewise merely C^. The analysis of [HZ. 2] in particular proceeds from 
a somewhat similar point of view to that followed here, and could be considered 
as an intermediate point between the analysis of [ZH] and that carried out here. 
The main technical innovation distinguishing the present treatment is a systematic 
method for refining high-frequency bounds on the resolvent kernel to any desired 
accuracy; for comparison, the zeroth order step in this algorithm gives bounds at 
the level of [ZH], and the first order approximately at the level of [HZ. 2]^, whereas 
here we require bounds to the second order. In effect, we generate an arbitrary order 
asymptotic expansion, limited only by smoothness of the coefficients of the linearized 
operator, of the resolvent about an essential singularity at infinity, with strong error 
bounds valid almost up to the boundary of analyticity/essential spectrum. This, 
and other techniques of the paper should apply also to systems with incomplete 
parabolicity, such as the equations of compressible gas dynamics or MHD. 

We now describe our results in somewhat greater detail. Consider a smooth 
traveling wave solution (1.7) of relaxation system (1.1). Besides (1.2)-(1.3), we 
make the following, rather standard assumptions (see, e.g., [Yo.l-3,YoZ,Z.4]): 

(HO) /, (7, Q G C^ in the general case, f,g,qE C^ for discrete kinetic models. 

(HI) a {{df, dgY{Uj f)) real, semi-simple with constant multiplicity, and different 
from s. 

(H2) a {df*{u±)) real, distinct, and different from s. 

(H4) The set of solutions of (1.7) forms a smooth manifold («'', v^)^ 5 E U eM.^. 



With the possible exception of the "noncharacteristic condition," s ^ cr^df^ dgY, 
these are precisely analogous to the hypotheses given in [ZH,Z.4] for the viscous case. 
Conditions (H1)-(H3), with 9 set to zero in (H3), and relaxing the requirement of 
distinct eigenvalues in (HI), are essentially the standard set of hypothesees proposed 
by W.-A. Yong for general relaxation systems [Yo.1-3]; as discussed in [Y.4], these 
are satisfied for most known relaxation sytems, and indeed may be considered as 
criteria for a "reasonable" model. A slight difference here is that we have assumed 
strict hyperbolicity/nonsonicity of the equilibrium system, (H2), and stability, (H3), 



■^In fact, we make also key improvements generalizing the approach of [HZ. 2] to the system 
case, and to operators with nonconstant principal part, see Remarks 8.16 and 4.10 below. 
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only at the endpoints {u±,v±) of the shock, and not everywhere along the profile; 
this allows applications to interesting nonhyperbolic situations such as are discussed 
in [JX,BGDV,BV]: in particular, shocks of nonclassical, under- and overcompressive 
type are admitted in our treatment, along with the classical. Lax type. Existence of 
weak (necessarily Lax-type) relaxation profiles has been shown under the hypotheses 
of [Y.1-3] together with some additional mild nondegeneracy conditions (satisfied 
for most physical systems) by Yong and Zumbrun [YoZ]; see [FreZ,BR] for related 
studies. In Appendix Al, we give a simplified treatment of existence, showing that 
existence of weak relaxation profiles is in fact guaranteed by (H0)-(H3) alone. 

The slightly strengthened condition (H3) {9 > 0) may be recognized as the 
dissipativity condition of Zeng-Kawashima [Ze.2,Kaw], and is likewise satisfied for 
a variety of physical models, though not all. The degenerate case 6* = is quite 
delicate, and requires additional hypotheses on the nonlinear structure of (1.1); for 
a discussion of this interesting situation and related phenomena, see [Ze.2]. 

Remarks. 1. The noncharacteristic condition s ^ ^{df, dgY is necessary in order 
that the traveling wave ODE be of nondegenerate type, a natural assumption in the 
context of stability of relaxation profiles. (See Section 1.1, below, for a statement 
of the traveling wave ODE). 

2. In (H0)-(H4), we have made no assumptions on the structure or amplitude of 
the profile {u, v). However, note that the assumption of a smooth profile is in effect 
a limitation on the strength of the shock, since, as is well known, profiles of suffi- 
ciently large amplitude will develop "subshocks," or jump discontinuities [Whi,L.2]. 
More precisely, this is a restriction on the amplitude with respect to the number of 
resolving modes in the system [BR, MR]: For example, the allowable amplitude for a 
moment closure system is a strictly increasing function of the number of moments, 
apparently approaching infinity in the limiting (infinite-dimensional) case of the 
Boltzmann equations. The stability of profiles containing subshocks is a very inter- 
esting problem for future investigation, and likewise the stability (or instability) of 
strong Boltzmann profiles. Of course, stability of even weak Boltzmann profiles is 
a fundamental open problem. 

Example 1.1. The Jin-Xin model. 



(1-8) ,. , 2 



Ut+Vj: = 0, 

Vt + a^Ux = h{u) 



u, V & M"^, was introduced in [JX] as a numerical scheme approximating solutions 
of the equilibrium system ut + h{u)x = 0. System (1.8) is in general nonstrictly 
hyperbolic, with a{df, dgY consisting of the two eigenvalues ±a, each with constant 
multiplicity n. It is readily checked that (H0)-(H3) hold if and only if: (i) h E C^; 
(ii) a j^ s; (iii) the equilibrium system ut + hiu)^ = is (nonstrictly) hyperbolic at 
the endstates u±, with s ^ a(dh{u±))] and, (iv) there holds at the endpoints of the 
shock the subcharacteristic condition a^ > a[dh{u±Y). 
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Likewise, one finds that the travehng wave ODE for a relaxation shock of speed 
s reduces to that for an "effective" viscous system, 

(1.9) ut + h{u)j, = {a'^ - s'^)uxj:, 

with constant, scalar viscosity a^ — s^, so that (H4) reduces to the corresponding 
requirement in the viscous case. Note that (iv) above implies that the effective 
viscosity a^ — s^ is positive whenever a profile exists, since in that case s must lie 
between the extreme values of the spectra of dh{u±). This furnishes a rich set of 
examples of relaxation shocks falling under our framework, through the existence 
theory for viscous traveling waves. Note that (1.8) is a discrete kinetic model, since 
its lefthand side is linear. | 

1.1. Structure/Classification of Profiles. We begin by presenting some 
basic, but apparently new results relating the structure of relaxation profiles to 
the characteristic structure of the corresponding equilibrium shock in the relaxed 
system (1.4). These generalize observations of Majda and Pego in the viscous, 
strictly parabolic case [MP]; analogous results hold for viscous profiles with real (i.e., 
partially parabolic) viscosity, as can be seen by essentially the same arguments used 
here (see, for example, the unified framework for relaxation/real viscosity models 
given in [Z.4], Appendices A1/A2). 

Without loss of generality taking s = 0, by a linear change in the flux functions 
if necessary, we may write the traveling wave ODE in the form: 

(1.10) , 

g{u,v) = q{u,v), 

where ' denotes d/dx. From (1.10), we immediately see that q{u±,v±) = and 
f{u) = constant, so that endstates {u±,v±) lie on the equilibrium manifold, and 
satisfy the Rankine-Hugoniot jump condition 

(RH) r{u+) - r{u-) = f{u+, v+) - fiu-,v-) = 

for the equilibrium system (1.4): that is, {u-,u+) corresponds to a shock solution 
of (1.4). Let us now consider the behavior of solutions of (1.10) near the rest points 
{u±,v±), for definiteness {u^,v+). 

By assumption (HI) together with the Inverse Function Theorem, {u, v) -^ (/, g) 
may be inverted in a neighborhood of {u+, v+), reducing (1.10) to the r-dimensional 
ODE 

(1-11) g' = q{u{f+,g),v{f+,g)), 

where / is held constant at value /_|_. Linearizing about the rest point g+, we obtain 
after a brief calculation the linearized equation 

(1-12) <7' = (?.,?.) ff" M (^r)g. 

where the righthand side is understood to be evaluated at {u+jV+). 
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Lemma 1.2. Under assumptions (H0)-(H3), the matrix 

Ju Jv 



(1-13) {qu:qvj . ^ I > . 

yu yv / \ -'r / _|_ 

has no center subspace, i.e., {u±,v±) are hyperbolic rest points of the reduced trav- 
eling wave ODE (1.11). In particular, traveling wave solutions (1.7) satisfy 

(1.14) \{d/dxf{{u{x),v{x))-{u±,v±))\<Ce-'^\''\, A; = 0, . . . , 4, 

as X ^- ±oo. 

Lemma 1.2 verifies the hypotheses of the Gap Lemma, [GZ,ZH]^, on which we 
shall rely to obtain the basic ODE estimates underlying our analysis in the low- 
frequency/Iarge-time regime. A proof of Lemma 1.2 is given in Appendix Al, a 
description and proof of the Gap Lemma in Appendix A3. 

We shaU classify relaxation shocks mainly according to their hyperbolic type, i.e., 
the type of the corresponding hyperbolic shock for the equilibrium system (1.4). 
Specifically, let i-\- denote the dimension of the stable subspace of df*{u^), i- 
denote the dimension of the unstable subspace oi df*{u-), and i := z_|_ + z_. Indices 
i± count the number of incoming characteristics from the right/left of the shock, 
while i counts the total number of incoming characteristics toward the shock. Then, 
the relaxed (hyperbolic) shock («_,«_!_) is called: 

Lax type if i = n + 1, 

Undercompressive (u.c.) if i <n, 
Overcompressive (o.c.) if i>n + 2. 

In case all characteristics are incoming on one side, i.e., i+ = n or i_ = n, a shock 
is called extreme. 

Similarly as in the parabolic case, [MP,ZH,ZS,Z.4], there is a close connection 
between the hyperbolic type of a shock profile, and the nature of the corresponding 
connection in the phase portrait of the r-dimensional reduced traveling wave ODE 
(1.11) on the n-dimensional submanifold 

/(w, v) — SU = f{U-, V-) — SU- 

of ]R"+^, obtained by integration of the -u-equation ((/(-u, t;) — sw) = in the full 
traveling wave ODE. Define (i+ to be the dimension within this submanifold of 
the stable manifold at («+, Vj^), and d- the dimension of the unstable manifold at 
(•u_, V-), and d := d- + dj^. Then, we have: 



^See also the version established in [KS], independently of and simultaneously to that of [GZ]. 
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Lemma 1.3. Under assumptions (H0)-(H3), there holds 

(1.15) d — r — i — n. 

Lemma 1.3 generalizes an analogous relation pointed out in the viscous case by 
Majda and Pego [MP] (see, e.g., Lemma 1.1 [Z.4]). However, the proof is different 
from the essentially algebraic one given in [MP], following instead an alternative 
argument suggested in Remark 2.3, [ZH], that is more closely linked to the dynamics 
of the traveling wave; for an exposition of this proof, see Appendix Al. 

A complete description of the connection, of course, requires the further index 
£ defined in (H4) as the dimension i of the connecting manifold between {u±,v±) 
in the traveling wave ODE. Generically, one expects that i should be equal to the 
surplus d — r — i — n. In case the connection is "maximally transverse" in this sense, 
i.e.: 

f 1 undercompressive or Lax case, 

(1.16) i=\ 

{i — n overcompressive case, 

we call the shock "pure" type, and classify it according to its hyperbolic type; oth- 
erwise, we call it "mixed" under/overcompressive type. Throughout the paper, we 
shall assume that (1.16) holds, so that all relaxation profiles are of pure, hyperbolic 
type. This holds in particular for the weak profiles whose existence was established 
in [YoZ] , which are always of pure Lax type. 

Remark. In the case of the Jin-Xin model (1.8), relation (1.15) follows directly 
from the corresponding parabolic result of [MP], see discussion. Example 1.1. 

1.2. Spectral vs. Linearized Stability. We next discuss spectral and lin- 
earized stability of profiles, and the relation between them: specifically, we show 
that spectral stability, appropriately defined, is necessary and sufficient for lin- 
earized stability. Linearizing (1.1) about the stationary solution {u,v), we obtain 
the linearized equations 

(1.17) Ut = LU:=-{AU)^ + QU, 
where 

and 

(1.19) U = (u, vY, ueW, ve W. 

The appropriate notion of stability of shock profiles is orbital stability, or con- 
vergence to the manifold of stationary solutions {{u^ ,v^)} defined in (H4). Like- 
wise, the relevant notion of linearized stability is that of linearized orbital stability, 
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defined as convergence to the tangent stationary manifold, consisting of station- 
ary solutions {{d / ddj){u^ ,v^)} of tfie linearized equations (1.17) [ZH]. An evident 
necessary condition for linearized orbital stability, (and thus for nonlinear orbital 
stability, appropriately defined) is 

(Dl) a{L) n {Re A > 0} = {0}, 

where a refers to the spectrum of L with respect to some L^, say L? for simplicity. 
We shall refer to condition (Dl) as strong spectral stability. 

A standard result of Henry [He] gives that the rightmost boundary of the essential 
spectrum aess{,L) is given by the right envelope of the union of the (purely essential) 
spectra of the limiting, constant-coefficient operators L± as x ^ ±oo, which by 
(H3) intersect the non-strictly-stable half-plane Re A > precisely at the origin 
A = 0. Thus, (Dl) is equivalent (as often the case for stability of traveling waves) 
to nonnegativity of the point spectrum: 

(1.20) Re ap(L) n {Re A > 0} = {0}. 

One might conjecture that sufficient conditions for linearized orbital stability would 
be (1.20) augmented with the requirement that A = be an £-fold, semi-simple 
eigenvalue of L, where I is the dimension of the tangent stationary manifold. How- 
ever, the lack of spectral separation associated with accumulation of the essential 
spectrum of L at the stationary eigenvalue A = prevents us from drawing imme- 
diate conclusions in this direction. Indeed, it is not a priori clear how one should 
properly define eigenspace at a point of accumulation of essential spectrum. 

As described in [ZH], such technical difficulties may be conveniently resolved 
using an extended, or "effective" spectrum based on the Evans function of [E.l- 
4,AGJ,GZ], defined precisely in Section 3, below. An analytic function playing a role 
for differential operators analogous to that played for finite-dimensional operators 
by the characteristic polynomial, the Evans function -D(A) has the property that, 
on the resolvent set of L, its zeroes agree (in both location and multiplicity) with 
the eigenvalues of L [AGJ,GJ.l-2]. Taking this correspondence as the definition 
of effective point spectrum, we can thus extend the notion of eigenvalue past the 
essential spectrum boundary of L, via analytic extension of the Evans function 
[GZ,KS]. The associated effective eigenprojection may then be defined in terms of 
the induced meromorphic extension of the resolvent kernel, via calculus of residues 
[ZH]. For details, see Section 5.4, below. 

The first main result of this paper, generalizing the corresponding result estab- 
lished for viscous, strictly parabolic, shocks in [ZH], is: 

Theorem 1.4. Under assumptions (H0)-(H3), shocks of "pure" hyperbolic type 
are L^ D L^ -^ L^ linearly orbitally stable for p > 1 if and only if L has precisely 
£ effective eigenvalues in {Re A > 0} (necessarily at X = 0), or, equivalently, there 
holds the Evans function condition: 

(V) D{-) has precisely £ zeroes in {Re A > 0} (necessarily at X = 0). 
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As pointed out in [ZH], the apparently simple criterion {V) in fact hides consid- 
erable subtlety in behavior, in particular in the critical case that (Dl) holds, but 
(V) does not: see [ZH], section 12 for examples of possible behaviors in this case. 

Evidently, (V) is equivalent to the pair of conditions: 

{VI) D{-) has no zeroes in {Re A > 0} \ {0}. 

{V2) (ii) {d/dXYD{0) ^ 0. 

The first condition may be recognized as (Dl), above. The second has been shown 
in [Z.4], Appendix Al to be equivalent to 

(D2) A := 7A ^ 0, 

where 7 is a coefficient measuring transversality of the connecting manifold be- 
tween {u±j v±) (i.e., 7 7^ equivalent to transversality of the intersection of unsta- 
ble/stable manifolds at {u-,V-)/{u+,v+)), and A is the one-dimensional "hyper- 
bolic stability coefficient" defined in [ZS,Z.4], appropriate to the type (Lax, under- 
compressive, or overcompressive) of the wave. For example, in the Lax case, A is 
just the Liu-Majda determinant 

(1.21) A :=det(rr,...,r:i,_,r;;+„..., <+,[«]), 

for which A 7^ is equivalent to hyperbolic stability of the corresponding ideal shock 
of the equilibrium system (1.4) [M.1-3]; here r* are the eigenvectors of df*{u±), 
in order of increasing eigenvalue, and i± are as defined above (Section 1.1). For 
the generalizations of (1.21) in the undercompressive and overcompressive cases, 
and their interpretations with respect to hyperbolic theory, see [ZS] or [Z.4]. (Note: 
the overcompressive version is cited also in (6.54)-(6.55) below). In the physical 
literature, transversality, 7 7^ 0, is referred to as structural stability (i.e., stability of 
inner shock structure), and evolutionary (hyperbolic) stability A 7^ as dynamical 
stability [BE]. 

Corollary 1.5. Under assumptions (H0)-(H3), linear L^ n L^ ^ L^ orbital sta- 
bility of pure shocks, for p > 1, is equivalent to (D1)-(D2), or equivalently: strong 
spectral stability plus classical structural and dynamical stability. 

Remark 1.6. The necessary conditions obtained in [Go,Z.l,Z.4] were signed 
conditions 

(1.22) 7A > 0, 

with appropriate normalization of 7, A. For 7A 7^ 0, the sign of 7A is associated 
with a stability index determining the parity of the number of unstable eigenvalues 
Re A > of L, hence yields information intermediate to that given in (Dl) and 
(D2). The advantage of (1.22) is that in the extreme. Lax shock case, the sign of 
the quantity on the lefthand side is explicitly computable as the product of (1.21) 
and a Wronskian determinant evaluated entirely at one infinity: two linear-algebraic 
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quantities. For history/applications/further discussion of the stability index, see e.g. 
[GZ,BSZ,ZS,Z.4]. 

Note, with the result of Theorem 1.4, that (1.22) is strengthened to 

(1.23) 7A > 0. 

Remark 1.7. The twin conditions (D1)-(D2) reflect multiple scales present 
in the problem, (Dl) corresponding to high-frequency and (D2) to low-frequency 
scales. As discussed in [Z.4], introduction, this dichotomy is natural from the point 
of view of formal, matched asymptotic expansion. For a generalization of (D1)-(D2) 
to multi-dimensions, see [ZS,Z.4]. 

Remark 1.8. In the case discussed in [YoZ], of weak shock profiles arising in 
regions where the equilibrium system is strictly hyperbolic, profiles are necessarily 
of Lax type, and transversal, 7 7^ 0; moreover, the Liu-Majda condition A 7^ 
is automatically satisfied, from standard hyperbolic considerations [M.1-3]. Thus, 
linearized stability in this case is equivalent to strong spectral stability, (Dl). As 
remarked earlier, strong spectral stability is equivalent to zero mass stability, so this 
is indeed a substantial reduction of the problem. For example, zero-mass stability 
has been shown in the case of weak Jin-Xin profiles [Liu,Hu], whereas linearized 
stability has not. The verification of (Dl) for weak relaxation profiles in this general 
setting we regard as an extremely interesting open problem in the one-dimensional 
theory. 

1.3. Pointwise Green's Function Bounds. Theorem 1.4 is obtained as a 
consequence of detailed, pointwise bounds on the Green's function G{x,t;y) of the 
linearized evolution equations (1.17) (more properly speaking, a distribution), which 
we now describe. For simplicity, we first restrict to the case that {df, dgY is strictly 
hyperbolic, indicating afterward the extension to the general case. 

Let aj{x), J = 1, . . . (n + r) denote the eigenvalues of A{x) = (df, dgY{u{x)), and 
lj{x) and rj{x) be smooth families of associated left and right eigenvectors, respec- 
tively, normalized so that /*rfc = d^.. These are well-defined, by strict hyperbolicity. 
Eigenvalues aj{x), and eigenvectors /j, rj correspond to hyperbolic characteristic 
speeds, and modes of propagation of the relaxation model (1.1). Likewise, let a 



*± 
j 



j = 1, ...,n denote the eigenvalues of df*{u±), and Ij , r* associated left and 
right eigenvector, corresponding to hyperbolic characteristic speeds and modes of 
propagation of equilibrium system (1.4) at the endpoints u± of the relaxation profile. 
Define local, scalar dissipation rates 

(1.24) r]j{x):=-llQrj{x), j = 1, . . . , n + r, 
and asymptotic, scalar diffusion rates 

(1.25) /?;±:=(/fSV;)^, j = l,...,n 



1 ' "•! 
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where B"^ := B*{u±) is the effective diffusion defined in (1.6), predicted by formal, 
Chapman-Enskog expansion. As described in Appendix A2, these quantities arise 
in a natural way, through Taylor expansion of the (frozen-coefficient) Fourier symbol 

(1.26) -iCA{x) + Q{x) 

of the linearized operator L about ^ = oo and ^ = 0, respectively. As a consequence 
of dissipativity, (H3), we have that 

(1.27) rif := rij{±oo) > 0, i = l,...,n + r, 
and 

(1.28) Pf>0, j = l,...,n. 

Let us define also aj± := aj{±oo). 

The important quantities rjj, 13* were identified by Zeng [Z2] in her study of 
decay to constant (necessarily equilibrium) solutions {u,v) = {u±jV±) of relaxation 
systems, corresponding at the linearized level to the study of the limiting equations 

(1.29) Ut = L±U := -A±U^ + Q±U 

as x -^ ±oo of the linearized evolution equations (1.17). Using Fourier Transform 
techniques, Zeng established the following sharp, pointwise bounds on the Green's 
function of (1.29), stated for reference in our own somewhat different notation: 

Proposition 1.9(Ze.2). Assuming (H0)-(H3), plus strict hyperbolicity of {df, dgY , 
the Green's function G{x,t;y) associated with the linearized, constant- coefficient 
evolution equations (1.29) may he decomposed as 

(1.30) G{x,t;y) = H + S + R, 

where 

H{x, t; y) := J^ e"''^"*5(x -y- aft)rfl 



n-\-r 

±l±t 
J 

i=i 



and 



S{x,t;y) := J^ (47r/3;±(l + t))"'el^-^-". ^''/^^i (^+*)i?f Lf *, 
i=i 

denote hyperbolic and scattering terms, respectively, and R denotes a faster decaying 
residual. Here, a^ , a* , Ij , I* , r^ , r* , rjj and /3* are as defined just above, 
andR*^, L*± by 



J ' J 



^*± \ / r*±t 



(^■'1' «" - -oA.--) ' ^" - \ 



-Qv QuT. 
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Note, as suggested by the derivation of constants ?]• , /3 in Appendix A2, that 
short-time (~ large-frequency) behavior is dominated by hyperbohc term H, while 
large-time (~ small-frequency) behavior is dominated by scattering term S. More- 
over, the latter term is precisely that identified in [LZe] as describing large-time 
behavior for the Chapman-Enskog approximation (1.5), in agreement with the gen- 
eral philosophy of [L.2]. 

Our second main result, and the principal result of this paper, is the following, 
nonconstant-coefficient analog of Proposition 1.9: 

Theorem 1.10. Assuming (H0)-(H4), plus strict hyperbolicity of {df,dgY, and 
the spectral stability criterion (V) (equivalently, (D1)-(D2)), the Green's function 
G{x,t;y) associated with the linearized evolution equations (1.17) may in the Lax 
or overcompressive case he decomposed as 

(1.32) G{x,t;y) = H + E + S + R, 
where, for y < 0: 

n-\-r 

(1.33) H{x,t;y) := ^ e-^^*(5._,^.t(-2/)r,(a:)/*(y), 

(1.34) 

E{x,t;y) := 



and 
(1.35) 

Six,t;y):=X{t>i} ^ K' Lr\47r(3-t)-'/'e-^^-y-<-'^"/^^'^-' 

a*-<0 



i/2-(x-z;-f/4p;-t I e 



e^ -|- e" 



+x„>n E \<ir\Fq-^^''Mmt)-"'<r'--'■i'^''^'i' (i4i^) 



al >0,a*+>0 



denote hyperbolic, excited, and scattering terms, respectively, and R denotes a faster 
decaying residual (rates given in Proposition 6.1, below). Symmetric relations hold 
for y > 0. 
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Here, the averaged dissipation and convection rates 'f]j{y,t) and aj{y,t) in (1.33) 
denote the time-averages over [0,t] of r]j{x) and aj{x), respectively, along charac- 
teristic paths Zj = Zj{y,t) defined by 

(1.36) dzj/dt = aj{zj), Zj{0) = y, 
while, in (1.35), 

(1.37) 4^^'"^:=«f f^- 
and 






.±1 u,i /.± 



2 



(1-38) P]^{x,t;y) := ^jf + ^ h^ I ^^ ' 

represent, respectively, approximate scattered characteristic paths and the time- 
averaged diffusion rates along those paths. In all equations, aj, a* , Ij, L* , 

rj, Rj , rjj and /3* are as defined just above, and scattering coefficients [c;^''_], 
i — —jOj -\-, are constants, uniquely determined by 



r> + 00 



{d/d6,)u'{s)d. 

a--<0 a;->0 J = ^ 

for each k = 1, . . .n, and satisfying 



^k 



a->0 a+<0 

for each j = 1, . . . , £, where the constant vector tvj is the left zero effective eigen- 
function associated with the right eigenfunction {d/d5j){u^,v^). Similar, but more 
complicated formulae hold in the undercompressive case (see Remark 6.9, below). 
In addition, G vanishes identically outside the hyperbolic domain of influence 

(1.41) Zl{y,t) <X< Zn+riy:t), 

Zj defined as above. 

Theorem 1.10 will be established by Laplace Transform (i.e., semigroup) tech- 
niques generalizing the Fourier Transform approach of [Ze.l-2,LZe]. The exten- 
sion of Fourier Transform techniques to the variable-coefficient case is of course 
a central issue in mathematical physics and analysis, associated with such topics 
as semigroup, Sturm-Liouville, and scattering theory [Ti,RS,LP]; however, many 
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unresolved issues remain, particularly in connection with higher order or nonnor- 
mal differential operators. From this, general theoretical point of view, the results 
of Theorem 1.10, and the methods used to obtain them, seem of wider interest, 
independent of their application to stability of shock profiles. 

Remarks. Just as in the constant-coefficient case, short-time/large-frequency 
behavior is dominated by hyperbolic term H, while large-time/small-frequency be- 
havior is dominated by terms E + S agreeing with those predicted for the partially 
parabolic Chapmann-Enskog approximation (1.5), or, equivalently, the "effective," 
strictly parabolic, system defined by Kawashima [Kaw]. 

Terms E and S are essentially those derived for strictly parabolic systems in 
[ZH,Z.2-3], combined with projection onto the equilibrium manifold. As discussed 
in those references, they may be understood by the caricature of Gaussian sig- 
nals propagating for x ^ according to the limiting, constant-coefficient equations 
at dzoo (of the associated effective, strictly parabolic system), interacting with an 
"ideal" shock layer at x = according to the rules encoded by scattering coeffi- 
cients [c;^'!j_], i = — ,0, +: That is, a signal initiating as a delta-function at y < 
first decomposes into n Gaussian signals traveling in each of the equilibrium char- 
acteristic modes r^~ with approximately equilibrium characteristic speed a^~; each 
of these signals then travels until it reaches the shock layer x = 0, whereupon it is 
refiected/transmitted in each of the outgoing characteristic modes r* , thereafter 
moving with approximate equilibrium characteristic speed a* ^0 along path z*j^ ; 
at the same time, the incoming signal excites stationary modes, or "bound states," 
{d/d5j){u^jV^). Incoming and outgoing Gaussian signals comprise the scattering 
term S, excited stationary modes the excited term E. 

The "parabolic" terms E + S are considerably more complicated in the variable- 
coefficient than in the constant-coefficient case; by contrast, the hyperbolic term H, 
since localized, has an appealingly simple generalization to the variable-coefficient 
case, consisting in the "frozen-coefficient" evolution described by (1.36). 

The nonstrictly hyperbolic case. The result of Theorem 1.10 holds also in 
the case that {df, dgY) has multiple characteristics, with the following modification. 
If aj has multiplicity rrij, j = 1, . . . , J, then /j(x), rj{x) are defined instead as 
(n+r) X nij matrices of left and right eigenvectors of (d/, dgYiu, v) [x] , normalized so 
that /yTj = Ijrij ■ In addition to this "static" normalization, we make the "dynamic" 
normalization 

(1.42) /5(^^.y(^)^0. 

This is always possible, as shown in Lemma 4.9 below; for discrete kinetic models, it 
may be accomplished simply by taking Ij, Tj constant. Then, defining the rrij x rrij 
matrix 

(1.43) ?7,(x):=-/;Qr,(x), 

and augmenting characteristic equation (1.36) with the dissipative flow 

(1.44) dCj/dt = -Vj{zj)Cj, Cj{y) = /m,, 
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governing the evolution of dissipation matrix Q = (j{y,t) G M"^j^"^j, we have: 

Theorem 1.11. Assuming (H0)-(H4) and the spectral stability criterion (V) (equiv- 
alently, (D1)-(D2)), decomposition (1.32) holds also in the general (nonstrictly hy- 
perbolic) case, with (1.33) replaced by 

J 

H{x,t;y) := J]rj(x)Cj(y, t)(5^-a^t(-y)/5(y) 

(1.45) 7^ 

= ^r,(x)0(e-^"*)5._,^.,(-y)/5(y), 
j=i 

for some uniform r]o > 0, with Ij, rj, and Q as defined just above. 

That is, the trivial, scalar flow ({t) = e-lo ~''j(^(«))^'^ = e"''-'* must be replaced 
in the case of corresponding characteristics by genuine, system dynamics (1.44). In 
the strictly hyperbolic case, (1.45) reduces to (1.33), as it must. Likewise, condition 
(1.27) must be replaced by 

(1.46) Rea{vf)>0, j = l,...,J. 

Remark 1.12. Similarly, we may relax the requirement of strict hyperbolicity 
of df*{u±), under a mild additional assumption. If a* has multiplicity m* , then, 
likewise, we define /* and r* as n x m* blocks, and 

(1.47) (3*^ := l*^^B*^r*^. 

Make the additional assumption that each (3j be diagonalizable, without loss of gen- 
erality diagonal. Then, (1.32) holds with (1.34)-(1.35) replaced by matrix versions 
as in (1.45) above. Diagonalizability of /?* follows, for example, from simultaneous 
symmetrizability of A, Q, which holds for most physical systems of interest; this is 
satisfied, in particular, for our main example, the Jin-Xin system [Yo.4]. 

1.4. Nonlinear Stability. Finally, we discuss our results on nonlinear stability. 
In this part of the analysis, we restrict our attention to weak Lax shocks of simul- 
taneously symmetrizable models (defined precisely in Section 7.2, below), or strong 
Lax shocks of discrete kinetic models (defined as above as the subclass of models for 
which /, g are linear, constant-coefficient). Simultaneously symmetrizable models 
of the dissipative type considered here admit favorable energy estimates [Kaw,Ze.2] 
allowing us to bound nonlinear derivative terms and close our nonlinear iteration; 
however, these energy estimates appear to require the weak shock assumption. Dis- 
crete kinetic models, on the other hand, do not have any such derivative terms to 
begin with, since nonlinear it ies enter only through the collision term q. 

With these restrictions, we obtain our third main result: 
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Theorem 1.13. Let {u{x),v{x)) he a Lax type relaxation profile of a general re- 
laxation model (1.1), satisfying (H0)-(H3) and the spectral stability criterion (V) 
(equivalently, (D1)-(D2)). If, also, {u,v) is sufficiently weak, in the sense that 
\{u',v')\l^ is sufficiently small relative to the parameters of (1.1), and (1.1) is si- 
multaneously symmetrizable (with coefficients q E C^ now also in the case of discrete 
kinetic models), then, provided that initial perturbation (uq^vq) — {u,v) is bounded 
by (o in Li cind H^ , for Co sufficiently small, the solution {u,v){x,t) of (1.1) with 
initial data {uo^vq) satisfies 

(1.48) \{u,v){x,t)-iu,v){x-5it))\L.<CCo{l + t)--^^'-'/p\ 

for all 2 < p < oo , for some d{t) satisfying 

(1.49) m\<cco{i+t)--^ 

and 

(1.50) \S{t)\<CCo. 

In particular, (-u, v) is nonlinearly orbitally stable from L^ fl H^ to L^ , for all p> 2. 
For weak shocks of simultaneously symmetrizable discrete kinetic models (q G C^) 
and data merely small in L^ HH^ , we obtain (1.48) forp = 2, but only boundedness 
|(w,f)(x,t) — {u,v){x — d{t))\Loo < C(o, for p = oo, and interpolated rates for p 
between 2 and oo. 

For discrete kinetic models, not necessarily simultaneously symmetrizable, with 
q G C^ , and data small in W^'^nW^'°^ , we obtain the rates (1. 48) for alll < p < oo, 
for shocks of arbitrary strength. 

Theorem 1.13 is obtained using a modified version of the argument used in [Z.2] 
to treat the strictly parabolic case, new complications arising from the need to 
control contributions of the singular component H of the Green's function. Similar 
results should be possible for nonclassical, over- and under- compressive shocks by 
more detailed, pointwise analyses in the spirit of [L.1,ZH] and [Z.5], respectively. 

Remark 1.14. As pointed out by W.-A. Yong [Yo.4], the condition of simul- 
taneously symmetrizability is satisfied for most systems of interest in applications: 
for example, the discrete kinetic models of Platkowski-IUner [PI] ; the numerical ap- 
proximation schemes of Jin-Xin [JX] and their generalizations by Natilini [N] ; the 
BGK models of Bouchut [B]; and, most notably, perhaps, for all of the extended 
thermodynamic models in the moment closure heirarchies of Levermore or Dreyer 
[Lev,Dre]; further examples may be found among the general relaxation models 
described in [CLL,Yo.l-3]. 

Remark 1.15. It is an interesting question whether the limitation to weak 
shocks for general relaxation models in Theorem 1.13 is only a technical artifact 
of our analysis, or a genuine difference in nonlinear behavior between quasilinear 
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and semilinear models. We point out that this hmitation arises very late in our 
argument, and at the (nonlinear) level of energy estimates only. On the other hand, 
discrete kinetic models are in some sense closer to original kinetic models such as 
Boltzmann equations than are quasilinear models based on moment closure assump- 
tions, hence it is conceivable that they truly give a more faithful approximation of 
behavior in the large-amplitude regime. 

Remark 1.16. Theorem 1.13 in particular shows that strong spectral stability 
implies nonlinear stability, for weak shocks of general simultaneously symmetrizable 
models, or for strong shocks of discrete kinetic models. As noted above, strong 
spectral stability follows for weak relaxation profiles of the Jin-Xin and Broadwell 
models from the results of [KM] and [Liu,Hu], respectively, thus yielding complete 
nonlinear stability results for these two cases. A very interesting open problem 
is spectral stability of weak shocks of general dissipative relaxation systems, as 
identified by Zeng/Kawashima [Kaw,Ze.2], or, likewise, spectral stability of strong 
shocks of special discrete kinetic models. 

Remark 1.17. By a more refined, pointwise analysis, it should be possible 
to show that d{t) approaches a limit as t -^ oo, which may then be explicitly 
determined using conservation of mass; however, it is clear that there can be no 
uniform rate of convergence for general L^nL°° perturbations, since interaction with 
the shock layer (and thus formation of the final asymptotic state) may be delayed 
for arbitrarily long time by moving the perturbation to spatial infinity. Likewise, it 
should be possible to show that the perturbed shock converges in L^ to a translated 
profile superimposed with diffusion waves in the far field: see, for example, related 
results obtained in [L.1,ZH] for localized data in the strictly parabolic case. 

As pointed out in [Z.2], the latter results, expressed in terms of the perturbation 
U{x, t) — U{x — 5(+oo)) from the time-asymptotic translate U{x — 5(+oo)) rather 
than the instantaneous translate U{x — S{t)), contain an additional error term 

(1.51) U{x - 6{+oo)) - U{x - 6{t)) ~ {6{t) - 6{+oo))U'{x) 

not found here, the so-called "coupled viscous wave" of [SX.1,S]. As observed from 
a different point of view in [S], the shape U' of this coupled wave depends on 
a "dynamical" effective viscosity induced by the traveling wave ODE, which is 
in general different from either of the effective viscosities B*{U±) predicted by 
Chapman-Enskog expansion at the endstates U± of the profile: a^ — s^, for example, 
vs. a^ — dh'^{u±), in the case of the Jin-Xin model (1.8). As discussed in [S], this 
is related to the modified Chapman-Enskog expansion of [L.2], in which derivatives 
are assumed small in directions of propagation of the solution, i.e., shock speed s 
in the shock layer vs. equilibrium characteristic speeds a* in the far fields: in the 
case of the Jin-Xin model, s vs. (j{dh{u±)). For the class of initial perturbations 
Uq — U r^ Co(l + |a^|)~"^^^ considered in [L.1,ZH], it can be shown that (5(t) — 
5(+oo)) ~ Co(l + t)~^^'^-, recovering the results of those papers. 
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Plan of the paper. In Section 2, we establish the basic spectral representation 
of G via inverse Laplace Transform of the resolvent kernel. In Section 3, we make 
a careful construction of the resolvent kernel; using this representation, we make 
careful resolvent kernel estimates in Sections 4 and 5, in high- and low- frequency 
regimes, respectively. Combining these results, we derive in Section 6 the Green's 
function bounds that comprise the main results of the paper, establishing Theo- 
rems 1.10-1.11 and a partial generalization yielding necessity of {T>) for linearized 
stability. Finally, in Section 7, we establish sufficiency of {T>) for linearized sta- 
bility, completing the proof of Theorem 1.4, and carry out the nonlinear iteration 
establishing nonlinear stability. Theorem 1.13. For completeness, we carry out in 
Appendices A1-A3 various auxiliary calculations/lemmas used throughout the pa- 
per. 

Note: Subsequent to the completion of our analysis, F. Rousset [R] has an- 
nounced the partial result of zero-mass stability of Jin-Xin relaxation shocks sat- 
isfying the spectral stability condition of [G,Z.1,Z.4] (the same one assumed here): 
that is, the extension to strong shocks of the zero-mass result of H. Liu [Liu]. How- 
ever, his analysis appears to be incomplete in the high-frequency regime; specifically, 
the Gap Lemma of [GZ,ZH] is incorrectly cited at a key step in the proof of his 
high-frequency resolvent estimates, see Remark 8.5, Appendix A3.1, below. This 
gap is remedied by our own, more precise bounds. 

Section 2. The Spectral Resolution Formula. 

The starting point for our analysis is the spectral resolution, or inverse Laplace 
Transform formula for the Green's function G{x, t; y) associated with the linearized 
evolution equations (1.17), defined by: 

(i) {dt — Lj:)G = in the distributional sense, for all t > 0; and, 

(ii) G{x, t; y) —r b{x — y) as t ^ 0. 
Here, G is to be interpreted in (i) as a distribution in the joint variables (x, y, t), and 
in (ii) as a distribution in (x,y), continuously parametrized by t. We shall see that 
G, uniquely so defined, by uniqueness of weak solutions of (1.17) within the class 
of test function initial data, in fact lies in the space of measures and not functions, 
so should properly be called a Green's distribution rather than a Green's function. 
Following the convention of [LZe,Ze.2], we shall use the two terms interchangeably. 

The Green's distribution may be conveniently constructed using standard semi- 
group theory, considering the linearized operator L as a closed, densely defined 
operator on L^, with domain H := {U : U e L'^, AU e H^}. 

Lemma 2.1. If A{x) is smooth and symmetrizable, in the sense that there exists a 
smooth, invertible S{x) such that SAS~^ is symmetric, and Q{x) is bounded, then 
operator L : LU :— —{AU)^^ + QU satisfies the bound 

(2.1) \U\l2 < |A - A*r^|(L - X)U\l2 
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for all U E TC, and real X greater than some value A* . 

Proof. Without loss of generality, we may take A symmetric, since we may 
achieve this by the coordinate change U -^ SU, without affecting boundedness of 
Q. Setting 

(2.2) (L - X)cP = /, 

for test functions (p G C^ , and performing an elementary energy estimate, we obtain 

(2.3) {cP\AcP) + {cP,QcP) = {cPJ) + X\cP\\ 

where (■, ■) and | ■ | refer to complex L^ inner product and norm, respectively. Taking 
real parts in (2.3), we have the Garding type estimate 

(2.4) -^((/>, A'c/)) + Re{(t), Qc/)) = Re{(t), f) + i?eA|(/>|2, 
or, rearranging: 

(2.5) {ReX-X.)\<pf<\<j>\\f\, 

for some real A* > sufficiently large. But, this gives the claimed bound 

(2.6) \(P\ <\X- A.rVl = |A - A.r^KL - X)(P\ 

for all A > Ah, on the real axis, for (p G Cq" within the class of test functions. 
Observing that C^ is dense within 7i (as may be seen by the straightforward 
computation \{AUy — Alf^lui ^ as e ^ 0, where f^ '■= f * ?7^, and rj^ G C^ 
denotes the standard moUifier, together with the standard facts that \{AU)'^ — 
AU\h'^ — > and jt/*^ — t/|^2 — ;> 0), we may pass to the limit to obtain the same 
bound ior (/) ETC. | 

Corollary 2.2. If the operator L : LU :— —{AU)nc+QU of Lemma 2.1 has in addi- 
tion the properties that A is nonsingular, and A and Q are asymptotically constant 
as X ^> ±oo, then L generates a C° semigroup e^* on L^ , satisfying |e^*|i2 < Ce^^ 
for some real u. 

Proof. The bound (2.5) applies also to the limiting, constant-coefficient oper- 
ators L± := —A±dx + Q± as X ^ ±oo, whence the spectra of these operators is 
confined to Re A < A#. Under the assumption that A is nonsingular, we can apply 
a standard result of Henry ([He], Lemma 2, pp. 138-139) to find that the essential 
spectrum of L is also confined to this set. Since (2.1) precludes point spectrum for 
real A > A^,, we thus find that all such A belong to the resolvent set p{L), with the 
resolvent bound 

(2.7) |(L-A)-^U2 < |A-A,r^ 

But, this is a standard sufficient condition that a closed, densely defined operator 
L generate a C° semigroup, with 1 6^*1^2 < Ce^^ for all real a; > A* (see e.g. [Pa], 
Theorem 5.3, or [Fr,Y]). | 

An immediate consequence (see, e.g. [Pa], Corollary 7.5) is: 
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Corollary 2.3. For L, e^* as in Corollary 2.2, there holds for all real rj greater 
than some rjo the inverse Laplace Transform (spectral resolution) formula 

/•r}-\-ioo 

(2.8) e^*(t) = P.V. e'^\L-X)-^(t)d\ 

Jij—ioo 

for all (j) G D{LP') C H^ , where domain D{L^) is defined as in [Pa], p. 1.^ 

Remark. Symmetrizability, in the one-dimensional context, is essentially equiv- 
alent to hyperbolicity: A has real, semisimple eigenvalues; in particular, it is implied 
by either strict hyperbolicity, or else nonstrict hyperbolicity with constant multi- 
plicity of eigenvalues, so that all cases considered here fall into this class. Thus, this 
is no real requirement other than the usual one of hyperbolicity of the principal part 
of the symbol (Note: this may be satisfied even when the subcharacteristic condi- 
tion is violated, or when the reduced system is not hyperbolic, as in applications of 
the Jin-Xin model [JX] to phase transition problems). 

With the result of Corollary 2.2, we can express the Green's distribution simply 
as 

(2.9) G{x,t-y):=e'^'5{x-y), 

where the evolution e^^(j){x, y) of a distribution (j){xj y) with the property that, for 
k G Co°, the distributional pairing ((/>(■, y), k{y))y is in L^ for each fixed x, is defined 
through its action on C^ test functions h{x)k{y): 

{e^*(j){x,y),h{x)k{y))^^y := {h{x),e^\(l){x,y),k{y))y)^, 

or 

(2.10) {G{x,t-y),h{x)k{y)),,y:= {h{x) , e"^' k{x)) , . 

The defining properties (i)-(ii) of the Green's distribution are then consequences of 
the corresponding properties defining a C° semigroup (see, e.g., [Pa], p. 1-4). 
Defining the resolvent kernel, similarly, through 

(2.11) Gx{x,y):={L-\)-H{x-y), 

where the righthand side is defined again through its action on test functions, 
((L - A)- V(x, y), h{x)k{y))^^y := (/i(x), (L - A)-i((/>(x, y), k{y))y)^, 

or 

(GA(a;,y),V'(a^,y))x,y := (/i(x), (L - A)"^/c(a:))^, 

we find through (2.8), as a tautology, the fundamental: 



^Note: From the definition in [Pa], we find that the domain D{L) of L satisfies H^ CH C D(L). 
Likewise, the domain D{L'^), consisting of those functions V such that LV G D(L), satisfies 
if2 c D{L'^). 
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Proposition 2.4. For L, e^* as in Corollary 2.2, and G, G\ as in (2.9)-(2.11), 
there holds the inverse Laplace Transform (spectral resolution) formula 

(2.12) G{x,t;y) = P.V. e'^^Gx{x,y) d\, 

J rj — ioo 

valid for all r] greater than some rjQ . 

It is this relation that we shall use to compute the time-evolutionary Green's 
distribution G. Note that, on the resolvent set p{L), the resolvent kernel G\ likewise 
has an alternative, intrinsic characterization as the unique distribution satisfying 

(2.13) {L^-\)Gx{x,y)=6{x,y) 

and taking f E LP' to {G \{x , y) , f {y)) y G L^ . This is the characterization that we 
shall use in computating G\. 

Remark 2.5. Note that a semigroup on even a still more restricted function 
class than L^ would have sufficed in this construction, since we are constructing 
objects in the very weak class of distributions. Later, by explicit computation of 
G, we will verify that L in fact generates a C^ semigroup in any L'p . Note also that 
(2.10) implies the expected, standard solution formula e^^/ = {G{x,t;y)j f{y))y, 
or, formally: 

(2.14) e'^'f = JG{x,t;y)f{y)dy, 

for / in any underlying Banach space on which e^* is defined, in this case on L^, 
I < p < oo. 

Remark 2.6. For general interest, we point out another, more concrete route 
to (2.12), (2.14), generalizing the approach followed in [ZH] for the parabolic case. 
Namely, we may observe that, in each finite integral in the approximating sequence 
defined by the principal value integral (2.12), we may exchange orders of integration 
and distributional differentiation, using Fubini's Theorem together with the fact, 
established in the course of our analysis, that Gx is uniformly bounded on the 
contours under consideration (in fact, we establish the much stronger result that 
Gx decays exponentially in \x — y\, with uniform rate), to obtain 

,,,,, (dt-L,) r ^ e^'Gxix,y)dX = 6ix-y) T ' e^' dX 

(^•J-OJ J-q-iK Jri-iK 

-.5{x-y)5{t) 

di,s K ^ oo, for all t > 0. So, all that we must verify is: (i') the Principal Value 
integral (2.12) in fact converges to some distribution G(x, t; y) as K ^> cxd; and, (ii') 
G{x,t; y) has a limit as t ^ 0. For, distributional limits and distributions commute 
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(essentially by definition), so that (i')-(ii') together with (2.15) imply (i)-(ii) above. 
Facts (i')-(ii') will be established by direct calculation in the course of our analysis, 
thus verifying formula (2.12) at the same time that we use it to obtain estimates 
on G. Note that we have made no reference in this argument to the semigroup 
machinery cited above. 

Section 3. Construction of the Resolvent kernel. 

We next derive explicit representation formulae for the resolvent kernel G\, using 
the classical construction (see, e.g. [CH]) for the Green's function of an ordinary 
differential operator, in terms of solutions of the homogeneous eigenvalue equation 

(3.1) Q = {L - \)W = -{AWy + {Q - \)W, 

matched across the singularity x = yhy appropriate jump conditions, in the process 
obtaining standard decay estimates on the resolvent kernel (see Proposition 3.4, 
below) suitable for analysis of intermediate frequencies A. Here, A and Q are as 
defined as in (1.17), W G M^ := M"^"'"^, and ' as usual denotes d/dx. In the details of 
our treatment, we follow [ZH,Z.4]; see also related analyses of [AGJ,K.l-2] regarding 
solution of the resolvent equation by variation of constants. 

3.1. Basic construction. First, consider the limiting, constant-coefficient 
eigenvalue equations 

(3.2) Q = {L±-\)W ■.= -A±W' + {Q±-\)W. 
Define 

(3.3) A:=UAf, j = l,...,n + r, 

where A^ denote the open sets bounded on the left by the algebraic curves A^- (^) 
determined by the eigenvalues of ^^(/)('u±,t'±) + (^ )(w±, v±) as ^ is varied along 
the real axis. 

Lemma 3.1. The limiting equations (3.2) have no center manifold on the set A. 
Moreover, assuming (HO)-(Hl), (H3), we have: 

("zj A C {A : Re A > -ryjlm A|/(l + |Im A|), r] > 0; and 

(a) the stable manifold at -I-cxd and the unstable manifold at — oo have dimensions 
summing to the full dimension N = n + r. 

Proof. The fundamental modes of (3.2) are of form e^^V, where /i, V satisfy 
the characteristic equation 

(3.4) {-nA± + Q±-X)V = 0. 

The existence of a center manifold thus corresponds with existence of solutions 
H = i^, V of (3.4), ^ real, i.e., solutions of the dispersion relation 

(3.5) i-iCA± + Q±-X)V = 0. 
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But, A G a{—i^A± + Q±) implies, by definition (3.3), that A lies outside of A, 
establishing the first claim. Bound (i) , likewise, follows easily from the bound (H3) 
on the dispersion curves A • . Finally, nonexistence of a center manifold, together 
with connectivity of A, implies that the dimensions of stable/unstable manifolds at 
+00/— oo are constant on A. Taking A -^ +00 along the real axis, we find from (3.4) 
that these dimensions are those of the unstable/stable subspaces of A+/A-; but, 
these must sum to A^, since, by (HI), the dimensions of unstable/stable subspaces 
of A{x) are independent of x. This establishes (ii), completing the proof. | 

For asymptotically constant-coefficient ordinary differential operators, the right- 
most component (A, in this case) of the set of A such that the limiting, constant- 
coefficient eigenvalue equations {L± — X)W = are of hyperbolic type, with di- 
mensions of unstable/stable manifolds summing to that of the full phase space is 
called in the Evans function literature the domain of consistent splitting [AGJ]. Its 
significance (see, e.g. [He], Lemma 2, pp. 138-139) is that it lies in the complement 
of the essential spectrum of the variable-coefficient operator L, with its boundary 
lying in the essential spectrum: that is, it is a maximal domain in the essential 
spectrum complement, consisting of the component containing real, plus infinity. 
Moreover, provided that the coeffients of L approach their limits at integrable rate, 
it can be shown (see, e.g. [AGJ,GZ,ZH]) that each connected component consists 
either entirely of eigenvalues, or else entirely of normal points, defined as resolvent 
points or isolated eigenvalues of finite multiplicity. The latter, and for our purposes 
more significant fact will be seen directly in the course of our construction. 

Recall, (1.14), that coefficients of L converge exponentially to their limiting val- 
ues in L-j-, in particular at integrable rate. In standard fashion (see, for example, 
[Co,AGJ,ZH,Z.l,ZS,Z.4]), we may thus conclude: 

Proposition 3.2. // (HO)-(Hl), (H3) hold, then, about any A in the domain of 
consistent splitting A, there exist locally analytic (in X) bases (pi , . . . , cp'^ {x; A) and 
(/)^_i_-i^, . . . , (/)~_^_^{x; A) respectively spanning the stable manifold at +00 and the unsta- 
ble manifold at —00 of solutions of the variable- coefficient equation (3.1). Moreover, 
these manifolds are tangent as x ^> —00, x -^ +00, respectively, to their constant- 
coefficient counterparts, decaying exponentially for x >Q, x < Q, respectively, with 
uniform rate Ce~^^^^, i] > 0. 

Likewise, there exist bases ip^, . . . , ip'j^{x] A) and V'ah-I' • • • ' '^n+ri^'i ^) respectively 
spanning the unstable manifold at +00 and the stable manifold at —00 of solutions of 
(3.1), that are locally analytic in X, tangent as x -^ —00, x -^ +00, respectively, to 
their constant- coefficient counterparts, and growing exponentially for x > 0, x < 0, 
respectively, with uniform rate Ce'''^', rj > 0. 

Proof. By spectral separation of stable/unstable subspaces, we may conclude 
the existence of globally analytic bases of solutions of the limiting equations (3.2), 
using a standard result of Kato ([Kat], pp. 81-83). The result for the variable- 
coefficient equations (3.1) then follows by standard, asymptotic theory of ODE 
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[Co,Z.l], using the integrable rate of convergence of coefficients to their hmiting 
values. Alternatively, using the fact that convergence is in fact exponential, we 
may apply the results of Appendix A2.1 (specifically. Remark 8.4) to obtain a self- 
contained treatment. | 

Observation 3.3. To show that A is in the resolvent set of L, with respect 
to any L^, 1 < p < cxd, it is enough to: (i) construct a resolvent kernel G\{x,y) 
satisfying (2.11) and obeying a uniform decay estimate 



(3.6) \Gx{x,y)\<Ce 



■r]\x-y\ 



and, (ii) show that there are no Lp solutions of {L — X)W = 0, i.e., A is not an 
eigenvalue of L. (For p < oo, these are necessary as well as sufficient, though we do 
not show it here [Z.l]). For, then, the Hausdorff- Young inequality yields that the 
distributional solution formula (L — A)~^/ := f G\{x,y)f{y)dy yields a bounded 
right inverse of [L — A) taking L^ to L^, while nonexistence of eigenvalues implies 
that this is also a left inverse. On the domain of consistent splitting. A, we shall 
see that (i)-(ii) are equivalent. 

Following the strategy of Observation 3.3, we are forced to to choose G(x, y), for 
fixed y and x ^ y from among the decaying solutions at +oo/— oo, respectively, of 
the homogeneous eigenvalue equations (3.1), in such a way the Gx satisfies ai x = y 
the jump condition 

(3.7) [Gx]^^y:=Gx{y + 0,y)-Gx{y-0,y) = -A-\y), 

where [h]x=y denotes the jump of the function h aX x — y: that is, according to 
the classical construction of [CH]. 

Implementing this strategy on the domain A of consistent splitting, we seek a 
solution of form 

. . . . f ^+(x;X)N+(y;X) x > y, 

3.8 Gx(x,y) = I 

where 

(3.9) $+(x;A) = (</>+(x;A)...</>+(a;;A))GM^x'= 
and 

(3.10) $-(x;A) = (<y^^+i(x;A) . . . (^^(x; A) ) G R^x^^-'^), 

(pj as defined in Proposition 3.2. Here, A^""" = N~^{y;X) and A^~ = N~{y;X) are 
k X N and (A^ — K) x A^ matrices, respectively, A^ := n + r, to be determined by 
imposing the jump condition (3.7), or equivalently 

(3.11) ($+(x;A) ^-{x;X))(^^^^Jy^^^^=-A-\y). 
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Inverting (3.11), we obtain the expression 

(3.12) {-N-ly%)=-^''^ ^-)-'iy;X)A-\y). 

for the matrices A^^, yielding a final expression for the resolvent kernel Ga, in terms 
of the matrices $^, of 

r-($+(a:;A) ) ( $+ $" )"' (y; A) A-i(y) x > y, 

(3.13) Gx{x,y)=l{0 $-(a:;A))($+ ^- y' {y; X) A-\y) x < y, 

valid when the matrix 

(3.14) $:=($+ $-) 

is invertible. If $ is singular, on the other hand, then A is clearly an eigenvalue, 
since there must be some function lying in the span of both $"*" and $~, hence 
decaying exponentially at both ±oo; thus, we see that (i) implies (ii) on A, partially 
verifying our claim that they are equivalent. 

It will at times be useful to represent the resolvent kernel in a more intrinsic, 
geometric fashion. Observing that 

(3.15) (0 $-(x;A)) = $(a;;A) (^° ^)=^x;X)^-\y;X){0 ^-{y;X)), 

and, similarly, 

(3.16) ($+(x;A) 0) = $(x;A)$-\y;A)($+(y;A) 0), 
we get the expression 

r -J^y^-U+A-\y) x>y, 

(3.17) Gx{x,y) = l J^v^-U-A-^y) x<y, 

where we have defined the solution operator from y to x of (3.1), denoted by J^y^^ , 
as 

(3.18) J^y^"" =^{x-X)^-^{y-X) 
and the projections H^ on the stable manifolds at ±oo as 

(3.19) n+ = ($+(y;A) Q)^-\y-X) and H" = ( $"(?/; A) ) ^-^(y; A). 

Note that in (3.15), there is no longer any reference to the specific choice of bases 
functions (j)- , but only to the (uniquely defined) stable and unstable manifolds at 
+CXD, — oo, respectively. 

With the above representations in hand, we readily obtain the following basic, 
intermediate-frequency result. More careful analyses will be required in the large- 
and small-freqency limits. 
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Proposition 3.4. With respect to any L^ , 1 < p < oo, the domain of consistent 
splitting, A, consists entirely of normal points ofL, i.e., resolvent points, or isolated 
eigenvalues of constant multiplicity. Moreover, the resolvent kernel is meromorphic 
on A, satisfying either of formulae (3. 13), (3. 15). Finally, on any compact subset 
K of p{L) n A, there holds the uniform decay estimate 

(3.20) \Gx{x,y)\<Ce-'^\-y\^ 

where C > and rj > depend only on K , L. 

Proof. As noted above, eigenvalues coincide witti zeroes of ttie (locally) analytic 
function 

(3.21) D{X) := det $, 

which are either isolated of finite multiplicity or else fill all of A. For the remaining 
claims, it is sufficient to verify (3.20) holds for the function Gx defined by formulae 
(3. 13), (3. 15). For, it then follows that this function, off of the set of isolated eigen- 
values, indeed represents the resolvent kernel, and that all noneigenvalue points 
are in the resolvent set. By inspection, (3.13) is meromorphic, with poles of finite 
order less than the multiplicity of the corresponding zeroes of D; it follows using 
the spectral resolution formula (see, e.g., [Kat,Y,ZH]) that the eigenprojection of 
L at each eigenvalue Aq, defined as the residue at Aq of the resolvent (L — A)~^, is 
well-defined, with range contained in the L^ kernel of (L — A)'^, where k is the order 
of the pole. But, the solution set of (L — X)^W = is at most nK, hence finite. 

Exponential decay is best seen using the geometric formula (3.15). First, observe 
that, for -D(A) bounded away from zero, the manifold of decaying solutions at — oo 
is uniformly transverse at x = to the manifold of decaying solutions at -f-oo, hence 
is uniformly exponentially growing as x -^ +oo: that is, away from eigenvalues of 
L, this manifold is globally uniformly exponentially decaying as x decreases, on the 
whole real line — oo < x < -f-oo. Likewise, the manifold of decaying solutions at 
+00 is globally uniformly exponentially decaying as x increases. A side-consequence 
is that the two manifolds are uniformly transverse as a; ^ ±oo, approaching the 
unstable/stable manifolds of the limiting coefficient matrix of (3.1); by continuity, 
they are uniformly transverse, then, on the whole real line. Thus, projections H^ 
are bounded in (3.15), and fiows J^v^^ are uniformly exponentially decaying, giving 
the result (3.20). g 

The transmission function -D(-) defined in (3.21) is known as an Evans function. 
It is possible, with more care, to define a globally analytic Evans function on A, 
and thereby to determine information on location of eigenvalues via topological 
considerations, see e.g. [E.1-4,J,AGJ,PW,GZ]. However, for our purposes local 
analyticity is sufficient. Though there exist as many versions of the Evans function 
as there are analytic bases (j)j ^ it is evident that these differ, locally, only by a 
nonvanishing analytic factor. Thus, location and multiplicity of zeroes agree for all 
versions of the Evans function, as the following proposition of Gardner and Jones 
shows that they must. 
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Proposition 3.5[GJ.l— 2]. On A, the zeroes of (any version of) the Evans func- 
tion D{-) agree in both location and multiplicity with the eigenvalues of L. 

Proof. It is evident tliat tliey agree in location, whence the result follows if one 
can show that zeroes and eigenvalues both split under regular perturbation of L; 
the existence of such a perturbation (a nontrivial assertion) was shown by Gardner 
and Jones in [GJ.l]. Alternatively, they showed in [GJ.2] that the correspondence 
may be seen by explicit computation in a convenient Jordan basis. | 

Note that Proposition 3.5 gives an alternative proof that eigenvalues of L are of 
finite multiplicity on A. 

Remark 3.6. In the constant-coefficient case, the matrices $^ above can be 
calculated (almost) explicitly. To accomplish this aim, we have to find the decaying 
modes at ±oo, that is we have to find independent solutions to 

(3.22) -AW' + BW-XW^O ^^ W ^ A'^B - XI)W, 
such that W{±oo) = 0. So, let n^ = ^^ (A) with i = 1, . . . ,N he the solutions of 

(3.23) p(;u) =det (A-i(S-A/)-;u/) =0, 

with Re /U^ < < Re /U~ for any z = 1, . . . , /c and j = k + 1, . . . ,N . Then, assuming 
/Uj 7^ /U for any i ^ j, (p^ = r^ e^i ^ where r^ is a right eigenvector of the matrix 
A~^{B — XI) relative to the eigenvalue n^ , so that 

(3.24) $(a:;A) = (r+e^^" ... r+e^^" r^+iC^^+i" ... r^e^^" ) G M^^^. 

Example 3.7. For the Jin-Xin model, the equations linearized at the stationary 
solution (-u, v) are 



(3.25) , , 2 



ut - su^ + t>^ = 0, 

vt + a^Ux — svx = a{x)u — v. 



where a{x) :— dh{u{x)). Hence 

(here any element in the matrices A and B corresponds to an n x n block). 

If a(x) = constant, then p((u) = H^^i [(a^ — s^)|U^ + [(1 + 2A)s — aj]|U — A(l + A)] , 
where aj denote the eigenvalues of a, so that there are n pairs of eigenvalues 

(3.27) 4W = «,-('+f)-T.,(A) _ 
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where crj{X) = ^/[aj — (1 + 2A)s]^ + 4A(1 + A)(a^ — s^). If the subcharacteristic 
condition holds, then a^ — s^ > and for any A > there holds /U^(A) < < /U~(A). 
The corresponding eigenvectors (for A > 0) are 

W^=C{ [{1 + X)s + (a^ - s^)ijf] r, , {ajS - a^X) r, )* (C ^ 0), 

where rj denote right eigenvectors of a associated with aj. Hence, 

(3.28) $(x;A) = (w^+e^^" ■■■ W^+e^»" W^ e^^"" ■■■ ^"6^""), 
and the resolvent kernel is given by 

(3.29) G^{x,y) = { ^^"^ ' ' _^ ^ 

where V- are row vector defined by 

(3.30) {W+ W-)-' = (^^''_y 

■ 

3.2. Dual formulation. For use in our later treatment of low- frequency behav- 
ior, we develop a modified representation of the resolvent kernel involving solutions 
of the adjoint eigenvalue equation 

(3.31) {L*-X*)W = 0, 
where 

(3.32) L*W := A*W' + Q*W 

denotes the adjoint operator for L, and * for a matrix or vector denotes matrix 
adjoint, or conjugate transpose. Following [ZH,Z.4], we first point out: 

Lemma 3.8. For any W, W such that {L* - X*)W = and {L - X)W = 0, there 
holds 

(3.33) (W, AW) = constant. 

Proof. By direct calculation, 

{W, AWy = {A*W', W) + {W, {AW)') 

= {-{Q*-X*)W,W) + {W,{Q-X)W) = 0. 
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From (3.34), it follows that if there are k independent solutions (pf^ . . . , (j)'^ of 
(L — XI)W — decaying at +oo, and N — k independent solutions 4>'i^j^i, • • • , 4>]^ of 
the same equations decaying at — oo, then there exist N — k independent solutions 
^^_^i^ ■ ■ ■ ,'i/'jv °^ (-^* ~ X*I)W = decaying at +C)0, and k independent solutions 
ifji , . . . , -i/;^ decaying at — oo. More precisely, setting 

(3.35) vl>+(x;A) = (^+^,(x;A) ■■■ ^'^(x; A) ) G M^^^^"^), 



(3.36) *-(x;A) = (^r(x;A) ■■■ V^(x; A) ) G M^^^ 
and 

(3.37) *(x;A) = (^-(x;A) *+(a;; A) ) G R^^^, 

where ip- are the exponentially growing solutions described in Proposition 3.2, we 
may define dual exponentially decaying and growing solutions tjjj and (p, via 

(3.38) (^ $)^A(* *)±=^- 

Recalling (see, e.g. [ZH], Lemma 4.2) the classical duality principle that the 
transposition G'^{yjx) of the Green's function G\{x,y) associated with operator 
(L — A) should be the Green's function for the adjoint operator (L* — A*), i.e., 

(L;-A*)Gl(a:,y) = 5(x-y), 

we may seek, alternatively, to represent the resolvent kernel in the form 

$+(a:;A)M+(A)§-*(y;A) x>y, 



(3.39) Gxix.y) 
^ ^ Av ,y; , $-(^._x)M-(A)^+*(y;A) x<y, 

where M+ = M+(A) and M" = M-(A) wekxk and (A^ - A;) x (iV - k) matrices, 
respectively, to be determined by imposing the jump condition (3.7). 
Solving, as before, we obtain 

/-M+(A) \ 

(3.40) M(A) := \=^-\z;\)A-^{z)^-^*{z;\), 

\ M-{\)l 

where 

(3.41) i|r:=(§- §+). 
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Note that the independence of the righthand side with respect to ^ is a consequence 
of the previous lemma. Thus, 

,^,^, ^, , r-($+(x;A) 0)M(A)(^-(y;A) 0)* x>y, 

(3.42) G\{x,y) = { 

1(0 $-(x;A))M(A)(0 vl/+(^;A))* x<y, 

with M as in (3.40). It is readily seen [ZH,Z.4] by duality, (3.33), that \E' is nonsin- 
gular if and only if $ is nonsingular, so that formula (3.42) is indeed valid on the 
same region as were the previous ones; see also Proposition 3.11, below. 

Finally, let us note as before that it is possible to represent the matrix G\ by 
means of intrinsic objects as solution operators and projections on stable manifolds 

-J^'^'=U+A-^{z)U-^'^y x>y, 

(3.43) Gx{x,y) = { T^^''U-A-^{z)U+^'^y x < y, 



where 

(3.44) ^^^y ^^-\z]X)^{y]X) 
and 

(3.45) n+ = ^-i(^;A)(^^+°.^^) and IIJ = ^"H^; A) ("^"^q^' ^^ ) . 

Remark 3.9. The exponential decay asserted in Proposition 3.4 is somewhat 
more straightforward to see in the dual formulation, by judiciously choosing z in 
(3.43) and using the exponential decay of forward and adjoint flows on a: ^ alone. 

3.3. Spectral decomposition. The matrix M(A) in (3.40) depends on the 
choice of functions (p- and ip- . For example, formula (3.42) may be considerably 
simplified in the constant-coefficient case by the choice of bases \1/ = $^, which 
yields 

(3.46) ^{x]X)A{x)^{x;X) = I 
and M(A) = /, or equivalently 

(3.47) G.(., .) = ( - ^-'« ^^^ '^^'''^'- " '"■'' 

[ Ej=i (pj (2^; ^)(pj *(y; ^) X <y, 

where (/> • , (^ • may be taken as pure exponentials 

(3.48) 4>f{x)~4>t%y) = e^^(^)(^-^V^±(A)V^±*(A). 



CORRADO MASCIA AND KEVIN ZUMBRUN 



33 



This appealing formula may be viewed as a generalized spectral decomposition. Note 
that, moving individual modes (/> (j), in the spectral resolution formula (2.12), us- 
ing Cauchy's Theorem, to contours /Uj(A) = i^ lying along corresponding dispersion 
curves A = Aj(^), we obtain the standard decomposition of e^* into eigenmodes of 
continuous spectrum: in this (constant-coefficient) case, just the usual representa- 
tion obtained by Fourier transform solution. 

Example 3.10. For the Jin-Xin model (3.25), the dual bases of (^^ can be 
obtained by inverting the relation \l/(y; A)A$(y; A) = /. Hence, recalling formula 
(3.28) and definition of V- in (3.30), we obtain 



^{y-\) 



+ 



fVi 



v+ 



V,-e- 



-mV V 



-M„y 



A*! y 



\ V- e-^~-y j 



A- 



(3.49) 



1 



(aV+-a)e-^+^ (;U+s - A)e-^+? 



(a^A - as)a{X) U^ _ a^^,_)e-^-v (A - ^,_s)e-^-y 



(we used (t(A) = {c? — s'^){fx- — /U+)). As noted before, for such a choice of a basis 
for the space of solutions to the adjoint equation, M(A) = / and the resolvent kernel 
Gx can be represented by 



(3.50) Gx{x,y) 

or, expanding, (3.47). | 



(0 $-(x;A))(0 ^+{y;X)r 
-($+(a;;A) 0){^-{y;X) 0[ 



X <y, 
x> y, 



In the asymptotically constant-coefficient case, the notion of continuous spectrum 
is in general not particularly useful outside of the self-adjoint case, see e.g. discussion 
of [OZ]. However, as pointed out in [ZH], the more primitive relation (3.47) persists 
in the form of a scattering decomposition encoding interactions of scalar modes at 
±oo. This decomposition gives a natural generalization of the Fourier transform to 
the asymptotically constant-coefficient case, one of the main contributions of [ZH] . 

The following result of [Z.4], a slight extension of Proposition 7.1 of [ZH], gener- 
alizes (3.47) to the asymptotically constant-coefficient case. 

Proposition 3.11 [Z.4]. On Ar\p{L), there hold 



(3.51) 



Gx{x,y) 



j,k 



M+{X)cP+{x;X)4^-{y;Xy 
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for y < < X, 

(3.52) GA(x,y) = 5^rf+(A)</>-(x;A)4-(y;A)*-5^^,-(x;A)^,-(y;A)* 

j,k k 

for y < X < 0, and 

(3.53) GA(x,y) = 5^rf-,(A)(/>-(x;A)^-(y;A)* + 5^</>-(x;A)<^-(y;A)* 

j,k k 

for X < y < 0, with 

(3.54) M+ = (-/,0)($+ $-)"^^- 
and 

(3.55) rf± = (0,/)($+ $-)"^^-. 

Symmetric representations hold for y > 0. 

Proof. The matrix M"*" in (3.40) may be expanded using duality relation (3.33) 
as 

M+ = (-/,0)($+ $-)-'yl-i(^- $-)*"' ("q 
(3-56) =(-/,0)($+ $-)-'(^- $")(q) 

^ I 2 

= (-/,o)($+ $-)"'^^, 

yielding (3.54) for x > y, in particular for y < < x. 

Next, expressing (p- (a;; A) as a linear combination of basis elements at — oo, we 
obtain the preliminary representation 

(3.57) G^{x,y) = Y,d%{X)<p-{x;X)^^{y;Xr + Y,e%^-{x;X)i^^{y;Xr, 

j,k j,k 

valid for y < X < 0. Duality, (3.33), with (3.17), and the fact that n_|_ = I — n_, 

gives 

f+)=(^- ^-)*An+vi>-,^ 

(3.58) =-($- ^-)"^[/-(0 $-)($+ $-)"^]^- 



K-l)-(o V)<*^ *-)-'*■ 
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yielding (3.52) and (3.55) for y < x < 0. Relations (3.53) and (3.55) follow for 
X < y < in similar, but more straightforward fashion from (3.33) and (3.17). | 

Remarks. 1. Note, in the various scattering decompositions given in Propo- 
sition 3.11, that all functions (/>■ , (p, , (p, , -0- are evaluated on the intervals 
[0, ±oo] on which their behavior is known. Indeed, by Proposition 3.2, they are 
well-approximated by corresponding solutions of the limiting, constant-coefficient 
equations at ±oo, i.e., exponential modes generalizing the continuous eigenmodes 
obtained by Fourier transform in the constant-coefficient case. 

2. Letting x, y ^ — oo in (3.53), we find that 

k k 

where (^^ , (pj^ denote solutions of the asymptotic constant-coefficient ODE of the 
eigenvalue and adjoint eigenvalue equations at — oo and G^ denotes the associated 
(constant-coefficient) resolvent kernel. Thus, Proposition 3.11 has the important 
and intuitively appealing implication that far-field behavior reduces to that of the 
constant coefficient case. 

3. The fact that d'^ = d~ may be deduced without calculation from the fact that 
the constant-coefficient formula (3.47) satisfies the resolvent kernel equation (2.11), 
but not the correct boundary conditions of decay at ±oo. The continuous d^ 
correction satisfies the homogeneous eigenvalue equation, serving only to enforce 
decay at -t-oo. As noted already, this term decays at — cx), as do the remaining, 
constant-coefficient terms J^k'^k^k- 

Section 4. High frequency expansion. 

We now turn to the crucial estimation of the resolvent kernel in the high frequency 
regime, as |A| — *> with Re A > —rj. According to the usual duality between spatial 
and frequency variables, large frequency A corresponds to small time t, or local 
behavior in space and time. Therefore, we begin by examining the frozen-coefficient 
equations at each value xq of x: 

(4.1) Ut = L,,U := -A{xo)U^ + Q{xo)U. 

As usual, we shall for clarity state first the result in the case that {df, dgY is strictly 
hyperbolic, indicating afterward the extension to the general case. 

Lemma 4.1. Let {df,dgY{u,v){xQ) have distinct, real eigenvalues. Then, for |A| 
sufficiently large, the eigenvalue equation {Lj;^ — \)W = associated with the frozen- 
coefficient operator L^^^ at xq has a basis of N = n + r solutions 
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that are analytic in 1/X, satisfying 4>j = e^^^'^'^°^^Vj{X,xo), with 

/Uj(A, xo) = -X/aj{xo) - r]{xo)/aj{xo) + C(l/A), 
^ ■ ' y,(A,xo) =r,(xo)/a,(xo) + 0(l/A), 

aj{-), r]j{-), and rj(-) as defined in Section 1.3. Likewise, the adjoint eigenvalue 
equation {L^^^—XyW = has a basis of solutions {(j)^ , . . . ,(j)j^ , 4>k+i^ • • • ? 4>n}{^') -^7 ^o) 
satisfying (pj = e~^^^^'^°^^Vj{X,xo), with 

(4.3) Vj{X,xo)=lj{xo) + 0{l/X), 

lj{-) as defined in Section 1.3. 

The expansions (4-2)-(4-3) hold also in the nonstrictly hyperbolic case, with (j)j, 

(pj and rj, Ij now denoting n x rrij blocks, and Hj, rjj denoting rrij x ruj matrices, 
J = 1, . . . , J, where rrij is the multiplicity of eigenvalue aj of A{xo). 

Proof. Similarly as in the proof of Lemma 3.1, this follows by inversion of the 
corresponding expansions about ^ = oo of the eigenvalues Aj(^) and eigenvectors 
Vj{^) of the frozen-coefficient Fourier symbol (1.26), using the basic relation Hj = i^ 
relating characteristic and dispersion equations (3.4) and (3.5). The expansions for 
the Fourier symbol are carried out in Appendix A2. | 

Example 4.2. (Jin-Xin model) For model (1.8), everything is almost explicit 
and expansion of the eigenvalues and eigenvectors can be computed directly. Re- 
calling that /U^ = 2^2- =F 20^ '^here cr(A) = y/a'^ + 4a^A + 4a^A^, a = dh{u{xo)) G 
^nxn^ we find that " 

I , , , A a(x) =Fa ^ /l\ , 

where n±, r]± are n x n blocks. Similar explicit expansions for Vj, Vj can be 
computed using the formulae V± = {n±, —A)* and V± = (±l/Aa(A))(— a^yU±, A)* of 
Examples 3.7 and 3.10. | 

By (3.47), Lemma 4.1 gives an expression for the const ant- coefficient resolvent 
kernel Gaxq at xq of, to lowest order: 
(4.4) 

ELi aj(xo)"^e(-^/"^-''^/"^)(^")(^-^Vj(xo)/5(xo) x<y, 



G 



\:x 







r^u 



EV=fc+i aj(xo)"^e(-^/"^-''^/"^)(^o)(^-^)rj(xo)/Kxo) x>y, 



-3=k+ 



Our main result of this section, and perhaps the main technical challenge of this 
paper, will be to establish on an appropriate unbounded subset of the resolvent set 
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p{L) that the variable coefficient solutions (j)- , cj)- "track" their frozen coefficient 
analogs to order 1/|A| in both direction and (exponential) growth rate, yielding 
order 1/|A| estimates for the resolvent kernel. More precisely, define 

(4.5) O := {A : -rji < Re A}, 

with ?7i > sufficiently small that O \ i?(0, r) is compactly contained in the set of 
consistent splitting A (defined in (3.3)), for some small r > to be chosen later; 
this is possible, by Lemma 3.1(i). 

Proposition 4.3. Assume that there hold (H0)-(H4), plus strict hyperbolicity of 
{dfj dgY . Then, for any r > 0, and rji = r]i{r) > chosen sufficiently small, there 
holds O \ -B(0, r) C A n p{L). Moreover, for R> sufficiently large, there holds on 
O\i?(0, i?) the decomposition 

(4.6) Gx{x,y) = Hx{x,y) + ex{x,y), 

(4.7) 

_ -Ef=.+ia.(y)-'e^^^-^/"^-''^/"^)(^)'^V,(x)/*(y) x > y, 

E-=i%(2/)-'e^'^"'/"^"'^/"^^^'^''r,(x)/*(y) x<y, 



(4.8) Qxix. y) = \-^B{x, y- A) + \-\x - y)C{x, y- A) + \-''D{x, y- A), 

where 

J2f=k+i e" ^"^ ^'"^ ^^^ '^^b't (3^' y) x>y, 



(4.9) B{x,y;X) 



Ej=i e ^^ ^/"^'(^^ '^%. (x, y) X <y, 



(4.10) 
C(x,y;A) 



EA'^ — f ^ \/ai(s)ds— f^ X/ai(s)ds 4- /■ \ ^ 

^,j=k+l^^^^ze[x,y]e Jy ' ''^ > J. / ^^ ; cj.{x,y]z), x>y, 

Ek — f ^ \/ai(s)ds— f^ \/ai(s)ds — /■ \ - 

. ,,.^;^mean^g[^^y]e ^v ' -^ > J. / ^ ^ c^-[x,y\z), x<y, 



with 



(4.11) \h%\cfA<Ce-^\''-y\ 



and 



(4.12) D(x,y;X) = { r.r. ., . .. 
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for some uniform 9 > independent of x, y, z. 

Likewise, there hold derivative bounds 
(4.13) 

{d/dx)ex{x, y) = (S°(x, y- \) + {x - y)Cl{x, y; A)) + A'^ {bI{x, y; A) 

+ {x- y)Cl{x, y; X) + {x - yfDl{x, y; A)) + X'^E^ix, y; A) 

and 
(4.14) 

{d/dy)ex{x, y) = (S;(x, y- X) + {x - y)Cl{x, y; A)) + A'^ (Sj(a:, y; A) 

+ {x- y)Cl{x, y; X) + {x- yfol^x, y- A)) + X'^Ey^x, y; A), 

where B^ and Cp satisfy hounds of form (4-9) and (4-10), Dl now denotes the 
iterated integral 
(4.15) 
D}j{x,y;X) = 

EN — f^ X/ah(s)ds— f ^ \/ai(s)ds— f^ \/ai(s)ds 

h,i,j=k+i^^^^y<w<z<xe Jy ' ''^ ' J™ / n ; j. / jkj 

^dh]j'j{x,y;z), x > y, 

Sy ^/<^h{s)ds-J^ X/ai{s)ds-J^ X/aj{s)ds 






^d1\Ax,y]z), x<y, 



with \d^\ I < Ce ^'^ ^1, and Ep satisfies a bound of form (4-12). 

The hounds (4-8)-(4-14) hold also in the nonstrictly hyperbolic case, with (4-7) 
replaced by 
(4.16) 

,, , , f -ELK+iaj(2/)"^e"-''« ^/"^■^^^'^Vj(x)Cj(x,y)/*(y) x > y, 

H\{x,y) = < ^ -T - f-x/ ( )d 

lJ2j=iajiy) e Jy 1"'^^' ^rj{x)Cjix,y)^{y) x < y, 

where Cj{x,y) G R'^j^'^j denotes a dissipative flow similar to Qj in (1.44), but with 
respect to variable x: i.e., 

(4.17) dCj/dx = r]j{x)Cj{x,y)/aj{x), Cj{y,y) = I, 

or, equivalently, (j{x,y) = Qjiy.r) for Zjiy.r) = x, where Zj as in (1.36) denotes 
the characteristic path associated with aj; rj, Ij now denote n x rrij blocks; and 
Hj, rjj denote ruj x ruj matrices, j = 1, . . . ,K , j = K + 1, . . . , J, where nij is the 
multiplicity of eigenvalue aj of A{xo), with 

ai < ■ ■ ■ < ax < < aK+i < ■ ■ ■ < aj. 

(Note that we have not here assumed (V)). 

Remark. The three error terms in (4.8) may be recognized as roughly analogous 
to the terms that would arise in the constant-coefficient case through expansions 
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(4.2)-(4.3); see also the related Fourier transform analysis of [LZe,Ze.l-2]. The 
complicated "mixed" form of the second, averaged term (corresponding to £j in 
Proposition 8.15, Appendix A3. 2. 3) comes from our search for minimal regularity 
hypotheses on f, g, q. At the expense of a further order of differentiability, allowing 
us to approximately diagonalize the eigenvalue equation to one higher order, we 
could obtain instead a term of form X~^{x — y)C{x, y; A), where, for example, 

k 

C{x,y-X) = y^e^;i-y<^-v,/a,Kz)dz 

(4.18) fr[ 

X (eS-(-'^)(-^)/^ - l)/((x - y)/X)r,{x)l'^{y) 

for X < y and symmetrically for x > y, precisely mimicking that of the constant- 
coefficient case. However, this distinction is unimportant for our analysis, in which 
the precise form of these terms plays no role; indeed, the (variational) form of (4.10) 
is somewhat more convenient for our purposes. 

The proof of Proposition 4.3, and its nonstrictly hyperbolic analog (stated at the 
end of the section) will be the work of the remainder of the section. 

4.1. General Framework. We have to solve, approximately, the eigenvalue 
equation 

{AWy = {Q-X)W 

as |A| ^ oo within the region O defined above: or, setting Z = AW, the equation 

(4.19) Z' = {Q-X)A-^Z. 

(As noted in [ZH] , such an invertible transformation is a trick to reveal the reduced 
regularity requirements on coefficients of divergence- form operators) . 

Following standard procedure (see, e.g., [AGJ,GZ,ZH,Z.4]), we perform the rescal- 
ing X = \X\x (suggested by the leading order behavior of /u(xo), (4.2)) to obtain the 
perturbation equation 

(4.20) Z' = {Ao{x) + \X\-^Ai{x))Z, 
where 

Ao:=-XA-WX\-^x), 

(4.21) Ai{x):=QA-\\X\-^£), 

X:=X/\X\. 

This expansion is to be considered as a continuous family of one-parameter pertur- 
bation equations, indexed by A G 5"^. We shall have to take some care, therefore, to 
ensure that all derived bounds depend continuously on the index A; this accounts 



40 STABILITY OF RELAXATION SHOCKS 

for the rather meticulous description of error terms in the relevant Proposition 8.15 
of Appendix A3. 2. 3, below. 

We seek to develop corresponding expansions of the resolvent kernel in orders 
of |A|~^, or, equivalently, by representation (3.17), expansions of stable/unstable 
manifolds of (4.20), and the reduced flows therein. This we will accomplish by a 
two-stage process. First, we will reduce by a series of coordinate changes to an ap- 
proximately block-diagonal system segregating spectrally separated (in particular, 
stable/unstable) modes, with formal error of a suitably small order. Second, we 
will convert the formal error into rigorous error bounds using a reflned version of 
the Tracking Lemma of [GZ,ZH,Z.1,Z.4]. 

Both of these steps will be carried out in a general setting suitable for applications 
to it arbitrary ordinary differential operators of nondegenerate type (i.e., for which 
the principal part of the symbol has flxed order and type) , even including operators 
that do not have asymptotically constant, but only smooth coefficients. Together, 
they give an algorithm for the systematic estimation of high-frequency behavior, to 
an order of accuracy depending only on the smoothness of the coefficients of the 
operator in question. 

4.2. Formal diagonalization procedure. More generally (fixing any indexing 
parameters such as A above), consider an ODE 

(4.22) W = A{ex, e)W (xeR), 
with small parameter e, where A has formal Taylor expansion 

(4.23) A{y,e) = J2e''A,iy) + 0{eP+'). 

k=0 

(Note that the case considered in (4.20) corresponds to linear dependence with 
respect to e.) Assume: 

(hi) A is CP+^ in y, e, or, if Aq, . . . , Ar are constant, then [d/deY+^A is C^"^ 
in y, e. In either case, all derivatives considered, in both y and £, are uniformly 
bounded for all y G M, £ < eo, for some Eq > 0. 

(h2) AQ{y) is block-diagonalizable to form 

(4.24) Doiy) = diag {do,i, . . . , rfo, J, 

dj G C"'^^"'^ , with spectral separation between blocks (i.e., complex distance be- 
tween their eigenvalues) uniformly bounded below by some constant 7 > 0, for all 
e < £07 3; G M. 

We begin by showing that (4.22) may be approximately block-diagonalized to 
order e^. 
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Proposition 4.4. Given (hl)-(h2), there exists a uniformly well- conditioned change 
of coordinates W = TW such that 

(4.25) W' = D{ex,e)W + 0{eP+^)W 

uniformly for x eM., e < Sq, with 

p 

(4.26) D{y,e) = J^e'Dkiy) + 0{eP+') 

k=0 

and 

(4.27) T{y,e) = J2e'niy) + 0{e^+'), 

k=0 

where Dj, Tj G CP~^^~^{y) for all < j < p. and each Dj has the same block- 
diagonal form 

(4.28) D,(y)=diagK- !,...,(/,• J, 

dj^k e C'"'=^"S as does Do . 

Our proof of this result relies on the following linear algebraic fact. 

Lemma 4.5. Let di G C"^i^"^i and d2 G C"^2xm2 j^^y^ norm bounded by Ci and 
respective spectra separated by I/C2 > 0. Then, the matrix commutator equation 

(4.29) diX - Xd2 = F, 

X G C"^ix"^2, is soluble for all F E C^i^^^, with \X\ < C(Ci,C2)|F|. 

Proof. Consider (4.29) as a matrix equation VX = T where X and JF are 
vectorial representations of the ra\ x m^ dimensional quantities X and F, and V 
is the matrix representation of the linear operator (commutator) corresponding to 
the lefthand side. The spectrum of T> is readily seen to be the difference 

{cti — 0^2 : OLj G o{dj)\^ 

between the spectra of d\ and d^^ with associated eigenvectors given by ri/2, where 
r\ is the right eigenvector associated with a\ and I2 is the left eigenvector associated 
with OL2- (Here as elsewhere, * denotes adjoint, or conjugate transpose of a matrix 
or vector). By assumption, therefore, the spectrum of T> has modulus uniformly 
bounded below by I/C2, whence the result follows. | 
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Proof of Proposition 4.4. Substituting W = TW into (4.25) and rearranging, 
we obtain 

W' = {T-^AT - T-^T')W 

= {T-^AT-eT-%)W, 

yielding the defining relation 

(4.30) {T-^AT - sT-^Ty) = D 

for T. 

By (li2), there exists a uniformly well-conditioned family of matrices Tq{x) such 
that Tq^AqTq = Dq, Dq as in (4.24); moreover, these may be chosen with the same 
regularity in y as Aq (i.e., the full regularity of A). Expanding 

(4.31) T-\y,e)=ll-{J2e'T-'T,) + {J2e'T,-'T,f-... \ T,-' 

\ J=2 j=2 / 

by Neumann series, and matching terms of like order e\ we obtain a hierarchy of 
systems of linear equations of form: 



(4.32) DoT,-'T, - T,-%Do - F, = D 



J' 



where Fj depends only on A^ for < k < j and T^, (d/dy)Tk for < j — 1. 
On off-diagonal blocks (/c,/), (4.32) reduces to 

(4.33) di[T-%Y'^''^ - [T-'T^]('^''^dk = Ff ''\ 

uniquely determining [Tq~ Tj]*^^''-*, by Lemma 4.5. On diagonal blocks (/c, /c), we are 
free to set [Tq~^Tj]*^^'^\ whereupon (4.32) reduces to 

(4.34) [Dj]^^^^'^ = -[Fj]^^^^\ 

determining the remaining unknown [Dj]^''''^'. 

Thus, we may solve for Dj, Tj at each successive stage in a well-conditioned way. 
Moreover, it is clear in the general case that Aq ^ constant that the regularity of 
Dj, Tj is as claimed, since we lose one order of regularity at each stage, through 
the dependence of Fj on derivatives of lower order Tk, k < j. In the case that 
Aqj . . . , As are constant in y, we may choose also To, . . . , T^ to be constant in y, 
again recovering the claimed regularity. | 

Remark 4.6. Proposition 4.4 may be verified in more indirect fashion without 
using Lemma 4.5, by expressing T instead as a product TqTi • ■ -Tp and solving the 
resulting succession of nearby diagonalization problems. 
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Remark 4.7. For relaxation systems, we will need to expand to formal error 
order £^, i.e., p = 1. Thus, we require A, Q E C^ in the general case, and C^ for 
discrete kinetic models, in which case Aq above is constant; this translates to our 
requirements of C^ /C'^ on /, g in (HO). 

Refined diagonalization procedure. Now, in addition to (hl)-(h2), assume: 

(h3) ^o(') decays exponentially in y. 
(Note: this is the case for operators whose coefficients approach constants at expo- 
nential rate as y = ex ^ ±oo, in particular for relaxation profiles). In this case, we 
point out that, by judicious choice of the initial diagonalizing transformation Tq, we 
may arrange that D is given to first order by the "frozen-coefficient" approximation 
To'^ATo. 

Proposition 4.8. Let (hl)-(h3) hold. Then, in Proposition 4-4? ^^ ^^ possible to 
choose To(-) in such a way that, also, Dq + eDi is determined simply by the block- 
diagonal part of 

(4.35) T^\Ao + eA^)To. 

Reviewing (4.30) and (4.34), we find that (4.35) is equivalent to the requirement 
that TQ^{d/dy)To vanish on diagonal blocks. The matrix Tq is given by 

(4.36) To = {Ru...,Rk). 

where the blocks Rj are made up of columns spanning eigenspaces of Aq corre- 
sponding to the various blocks. Let us write suggestively 

(4.37) T,-' = I : 

where LjRk is a zero block if j ^ k and an identity block if j — k. Then, Proposition 
4.8 follows from: 

Lemma 4.9. Let (hl)-(h3) hold. Then, given any initial prescription of Rj{0), 
Lj{0) at y = 0, there exists a unique choice of Rj, Lj in (4.36)-(4.37) such that 

(4.38) LjR'j = L'^Rj = 0, 

which, moreover, is uniformly smooth and bounded to the same order of differentia- 
bility as is Aq. 

Proof. Let Rj, Lj denote any smooth choice, and define 

(4.39) Rj--=RjCyj L, :=a-iL„ 
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where aj are nonsingular kj x kj matrices to be determined, kj denoting the dimen- 
sion of the jth block (i.e., Rj dimension A^ x kj). Defining aj to solve the linear 
ODE 

(4.40) «;. = -Lj R'^ a, , a, (0) = h^ , 

where ' in this single instance denotes d/dy, we obtain (4.38) by direct calculation: 
(4.41) 

L,R'^ = L,{R'^a, + R,a'^) = Lj{R'^aj - RjL.R'^a^) = L,{I - RjL,)R'^a, = 0. 

The other relation comes from LjRj = /, from which L'.Rj + LjR'- = 0. Uniform 
global bounds and smoothness then follow by linear theory, and the fact that Rj , Lj 
could be chosen in the first place so that LjR'j decays exponentially in y at spatial 
infinity. | 

Remark 4.10. The key result of Lemma 4.9 was first pointed out in the diag- 
onalizable case by Goodman, in [Go.l]. There, the careful choice of diagonalizing 
coordinates was essential in his analysis of stability by energy methods. It is es- 
sential for us here, as well, in this different context. Note that we do not need this 
result in the case that Aq is identically constant, corresponding to the case that the 
principal part of the symbol of the associated operator L has constant coefficients, 
for which we may simply choose L, R constant as well. This was done, for exam- 
ple, in the treatment of scalar dispersive-diffusive systems in [HZ. 2], for which the 
dispersion coefficient was assumed constant, and can be done also in the case of 
discrete kinetic (relaxation) models. Thus, the role of Lemma 4-9 in our analysis 
is to give sharp estimates in the case of a variable- coefficient principal part. 

Remark 4.11. An important observation for us is that, in both Propositions 
4.4 and 4.8, any additional regularity of A with respect to e translates directly to 
regularity of the error term 0{e^^^). More precisely, if A G C'^{e -^ C^^^{y)), q > 
p + 1, then the error term has regularity C"^~^~^(£ -^ C^{y)). In case A G C^(£ -^ 
CP^^(y)), (J > p -t- 1, as often happens (for example, for relaxation profiles), the 
error term has regularity C"^(£ -^ C^{y)). This observation is used in conjunction 
with Proposition 8.15 below; specifically, it gives regularity with respect to 6 (in 
this case e^) in the error term On. 

Remark 4.12. A second important observation is that even though we have 
used in (4.20) the (standard) homogeneous rescaling x = \X\x, which introduces 
terms such as A = A/|A| that are singular at A = oo, any associated singularities 
at infinity in the resulting expansion are only apparent. That is, they are coordi- 
nate singularities that disappear when the diagonalized system is converted back 
to original coordinates. More precisely, since \X\A is analytic at oo with respect to 
A~^ = e/A, we may choose the initial diagonalizing transformation Tq to be ana- 
lytic in A~^ as well (indeed, it may be chosen independent of A; see (4.21), with 
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the effect that |A|iI>o is also analytic in A~^. If we make this choice, then it is 
readily verified that all further Tj, Dj also enjoy these properties, since they are 
preserved by the (uniquely determined, from this point) diagonalization procedure. 
Clearly, this observation is not unique to the case of relaxation systems, but holds 
for general equations and rescalings. An immediate consequence is that the rigorous 
bounds obtained from the formal diagonalization series by the fixed point iteration 
described in Appendix A3. 2, when expressed in original coordinates, must also be 
analytic in A~^; for, the fixed point iteration, being scale invariant, may just as well 
be carried out in the original (i.e., unrescaled) coordinates to begin with, whereupon 
we obtain analyticity as usual from the uniform convergence of iterates. 

We point out that we could alternatively have performed the formal diagonaliza- 
tion to begin with in variable e := A~^, using the modified perturbation formulation 

W' = X{-A-^+X-^QA-^), 

and taking Aj to depend on {X~^ex,e). 

4.3. Error bounds/application to relaxation systems. We now return 
to our original problem, that of approximating the stable/unstable manifolds at 
+oo/— oo of ODE (4.22), and the reduced fiows therein. In Appendix A3. 2, we 
present a systematic method for generating such approximations, given a system in 
the approximately diagonal form (4.34). The main requirement in the theory is that 
the approximate fiows JF^^^ generated by the blocks dj := D^^'^' be "exponentially 
separated" to some order e^ , k < p, in the sense that, as x increases, they are either 
uniformly exponentially decaying or else uniformly exponentially growing, i.e.: 

(4-42) <^ ^ 

for X > y, with 

(4.43) 6*1 < < 02. 

Observation 4.13. In the simplest case, that Aq is diagonalizable, the verifica- 
tion of (4.43) is straightforward. For, in this case, the block flows are scalar, and 
exponential separation is equivalent to the integral spectral condition 

f r Re d,{z) dz < e'^Oiix - y) + 0(1), 

(4.44) < 

\j^Redj{z)dz>eH2{x-y) + 0{l), 

or, roughly, the satisfaction on average of a spectral gap Re dj < e'^Oi/'Re dj > e^92- 
Moreover, this is clearly equivalent to the satisfaction of the corresponding condition 

, , , i;Re{Elo^'di,,{z)dz<e'^9i{x-y) + 0{l] 

(4.45) ' ^ 



/; Re (Eto ^'di,j{z) dz > eH2{x - y) + 0{l) 
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at the kth stage of the diagonahzation process; thus, we need only check (4.45) for 
each 1 < k < p. In particular, it is sufficient that there hold the pointwise spectral 
gap 

(4.46) <^ ^^i=° '' - 

lRe(Eto£'^^j>£'^2, 

on all except a set of measure l/ct, where it is violated by at most a; for, the 
corresponding error can be absorbed in the multiplying constant C of (4.44)-(4.45). 
In the case of strictly hyperbolic relaxation systems satisfying (H0)-(H4), (4.46) 
holds on n for k = p = 1 for all \y\ > M, or equivalently \x\ > M/e, by Lemma 4.8 
together with Lemma 4.1, (1.27), and the exponential approach in y of coefficients 
to their constant states. Lemma 1.2, while on \y\ < M, or |a;| < e, it is violated 
by at most order e. Thus, the entire theory of Appendix A3. 2 applies, yielding the 
error bounds reported below. | 

Observation 4.14. When Aq is only block-diagonalizable, the verification of 
(4.43) is more involved. Since dj is now a matrix, (4.44) is neither necessary nor 
sufficient. Moreover, the expansion (4.26) must be considered as a matrix per- 
turbation problem, and so the dependence of spectrum on coefficients D^ is not 
necessarily monotone in k. Thus, (4.42)-(4.43) are in general not equivalent (as in 
the scalar case) to the corresponding conditions on the kth order approximate flow, 
unless there hold additional properties such as diagonalizability of X]i<fc-i ^''d'lj (^^ 
implied, in particular, by symmetry) or special structure of dij relative to dk-ij, 
for all k < I < p. 

On the other hand, suppose as in the present case that dij are scalar until step 
k, where splitting occurs. Then, the contribution of these lower order terms to the 
approximate flow induced by dkj may be split off as an exponential factor 

times the approximate flow induced by e^dkj alone, simplifying the computation 
somewhat. In particular, for nonstrictly hyperbolic relaxation systems satisfying 
(H0)-(H4), with /c = p = 1, we obtain just e~-'y '°'^^^' ^ times the flow induced 
by sr]{z), tj as deflned in Section 1.3. The latter flow may be seen to have an 
order e spectral separation, by (1.46) plus exponential convergence of rjj to tj- as 
y -^ ±oo^, while the former has exponential separation (in agreeing sense) of order 



^For example, we may apply the Gap Lemma of Appendix A3.1 separately to cases x, y £ 
(—00,0] and x, y £ [0, +oo), treating the case y < < a; by the semigroup property as the 
concatenation of cases y < = x and y = < x. Alternatively, we may proceed as in Appen- 
dix A3. 2 by applying a change of coordinates Wj = Pwj, P = I + 0{e^^^^), converting r]j to 
positive/negative definite form for all x > M/e). 
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Re A/mirijy \aj{y)\. Combining these observations, we find that there is spectral 
separation on aU of O. Thus, the theory of Appendix A3. 2 apphes again, yielding 
the desired error bounds. | 

Remark 4.15. More generally, we point out, in the case relevant to traveling 
waves that all Aj converge exponentially in y to constants A- as y ^ ±oo, that, 
provided that dij are scalar, and there is neutral separation to zeroth order 

Re dij > 0, 
Re dij < 0, 

for all y (as is then a necessary condition for spectral separation along rays pointing 
into the domain of consistent splitting) , then a sufficient condition for spectral sep- 
aration of order e^ is spectral gap of order e^ in the limiting blocks dj {±oo) (also a 
necessary condition, by the Gap Lemma of Appendix A3.1). For, the diagonalizing 
construction described in Proposition 4.4 clearly preserves the property of expo- 
nential decay to constant states. Thus, making a constant change of coordinates 
such that dj becomes positive/negative definite of order e^, we have that dj is pos- 
itive/negative definite to order e^ + 0{ee'^^), readily yielding (4.42). More precise 
accounting of the separation, as in the two observations just above, requires more 
detailed consideration of the flows of different orders. We remark that dij scalar 
is roughly equivalent to diagonalizability/constant multiplicity of Aq, since it may 
then be arranged that each block dj correspond to a single eigenvalue of Aq. Thus, 
the requirement of spectral separation is in practice no restriction for applications 
to stability of traveling waves. 

Proof of Proposition 4.3. We will carry out in detail the basic estimate (4.8), 
indicating the extension to derivative bounds (4.13)-(4.14) by a few brief remarks. 
First, observe that in the modified coordinates Z = AW, (3.17) just becomes 



A-i(x)J^^^^n- x<y 



(4.48) Gx{x,y) 



-A-^{x)J^y^''U+ x>y 



where II-i-, J^v^^ now denote the corresponding objects in Z-coordinates. 

Applying the formal diagonalization procedure of Proposition 4.4 to (4.20), and 
choosing the special initialization of Lemma 4.9, at the same time taking care as 
described in Remark 4.12 to preserve analyticity with respect to A~^ in original 
coordinates, we obtain the approximately diagonalized system 

z'j = {Xaj + er)j)zj + £^ ^ ©ifc^fc 

= |A|"^ I {Xaj + r]j)zj + X~^^Oj,kZk 
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where error terms Gj jt are bounded and analytic in A~^ as A goes to infinity, or, 
converting back to the original spatial coordinates: 

(4.49) z] = (Xaj + 7]j)zj + A'^ J] Qj^kZk 

\ j,k 

Moreover, the approximatedly diagonalizing transformation T is given by 

T = To + A-^Ti 

= (ri,...,rj)+A-iTi, 

where Tj are as described in Proposition 4.3. 

As described in Observations 4.13-4.14 just above, the approximate, block- 
diagonal flows Zj = ((Attj + 'r]j)zj are exponentially separated, with uniform rates 

f C^y^"" ^/aKiz)dz^-e\x-y\^ 

where a^ < < ax+i- Applying the fixed point construction of Appendix A3. 2. 2(ii), 
we find that the stable/unstable manifolds of (4.49) are uniformly transverse for 
A G O, |A| sufficiently large, with associated projections II-i- given by 

K 

(4.51) n+ = J] r^l* + X-^E+ + C>(A-2) 

and 

J 

(4.52) n_ = J] rjl* + \-^E- + C(A-2), 

3=K+1 

for some continuous matrices E±{x). 

From transversality, we obtain immediately that O \ -B(0, R) C p{L) n A, for 
R sufficiently large, as claimed. Moreover, approximating the reduced fiows on 
stable/unstable manifolds using Proposition 8.15 of Appendix A3. 2. 3, and recall- 
ing representation (4.48), we obtain after rearrangement the claimed estimates on 
the resolvent kernel G\, where error terms B correspond to error terms i?± in 
(4.51)-(4.52), C to term Sf in (8.61)-(8.62) of Proposition 8.15, and D contains aU 
remaining, 0(A~^) errors: for example, the third, 0{5'^ /rj^) — 0{\~'^) term in the 
righthand side of (8.61), the 0{\~'^) terms in (4.51)-(4.52), and cross-terms coming 
from the product of different order A~^ terms. 

Derivative estimates. To obtain derivative estimates (4.13)-(4.14), we may follow 
precisely the same procedure, but carrying out to one further order the fixed point 
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iteration at the stage of the reduced flow estimate described in Proposition 8.15: 
that is, to order A~^, with error term 0{X~'^. (Note: this leads to the second-order 
iterated integral (4.15)). This yields a Taylor expansion to the same higher order for 
the X and y derivatives of the reduced flow, within the resettled spatial coordinates, 
and thus, in original coordinates, to the desired A~^ order, with 0(A~^) error. | 

Section 5. Low frequency bounds. 

Our next object is the estimation of the resolvent kernel in the remaining critical 
regime |A| ^0, corresponding to large time behavior of the Green's function G, or 
glohttl behttvior in space and time. Here, the details of the system are not important, 
and we may follow essentially the same treatment as in the strictly parabolic case, 
combining aspects of [ZH,Z.2,Z.4]. From a global perspective, the structure of the 
linearized equations is that of two nearly constant-coefficient regions separated by 
a thin shock layer near x = 0; accordingly, behavior is essentially governed by the 
two, limiting far-field equations: 

(5.1) Ut = L±U, 

coupled by an appropriate transmission relation at a: = 0. 

Similarly as in the previous section, we begin by examining these limiting sys- 
tems. For simplicity of exposition, we make the provisional assumption: 

(Al) The r X r matrices dqA'^ (0,/^)* have distinct eigenvalues. 

(Recall, (1.12), that these matrices govern behavior near {u±, v±) in traveling wave 
ODE (1.10).) Since none of our estimates will depend on the distance between 
eigenvalues, this assumption is easily removed by a limiting argument. 

Lemma 5.1. Assume thttt there hold (H1)-(H3) ttnd (Al). Then, for |A| suffi- 
ciently snittll, the eigenvttlue equtttion {L± — X)W = ttssocittted with the limiting, 
constttnt- coefficient opertttor L± htts tt bttsis of n + r solutions 

(5.2) W^ = e^f^^^''Vj{X), 
/U • , Vj , ttnttlytic in X, consisting of r "fttst" modes 

(5.3) ^--rf^OiX,. 



V,± = 4±'S± + ©(A), 



Sj = (0, Sj *)* where 'jj , Sj ttre eigenvttlues ttnd ttssocittted right eigenvectors of 
dqA'^{Q,Iry (equivttlently, 7- , S- ttre nonzero eigenvttlues ttnd ttssocittted right 
eigenvectors ofQA~^), ttnd n "slow" modes 



(5.4) 



± '^'- ^''-*± + A2/3*±/a*±' + 0(A3), 
V,%^{X):=Rf + 0{X), 
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where a* , R* , and /3*± are as in Proposition 1.9. Likewise, the adjoint eigenvalue 
equation {L± — X)*Z = has a basis of solutions Wi = e~^J *■ '^Vj{X), with 

(5.5) V,{X)=f^ + 0{X) 
and 

(5.6) Vr+j{X) = L*+0{X), 

V analytic in X, where Tj are the left eigenvectors of QA^^ associated with the 
nonzero eigenvalues — /U • , and L* are as in Proposition 1.9. 



Proof. Equivalently, we must show that /U , V- describe the eigenvalues and 
associated eigenvectors of coefficient matrix A~^{Q — A)±, and (see Lemma 3.8) 
that V- describe the left eigenvectors of {Q — X)A'^ . At A = 0, there are r distinct 
nonzero eigenvalues fi^ = 'jj , by assumption (Al), and an n-fold eigenvalue /U = 
corresponding to the n-dimensional kernel of Q. By spectral separation, the former 
vary analytically as A is varied, with Taylor expansions (5.3) and (5.5). 

The latter may be seen to vary analytically, with expansions (5.4) and (5.6), by 
inversion of the expansions 

x,m = -iafc-pfe+oie) 

carried out in Appendix A2 for the the dispersion curves near ^ = 0, together 
with the fundamental relation n = i^ between roots of characteristic and dispersion 
equations (recall (3.4)-(3.5) of Section 3.1). Alternatively, one can carry out directly 
the associated matrix bifurcation problem as in [ZH] to obtain the same result. | 

Our main result of this section is then: 

Proposition 5.2. Assume that there holds (H0)-(H4) and (D) (equivalently, (Dl)- 
(D2)). Then, for r > sufficiently small, the resolvent kernel G\ has a mero- 
morphic extension onto B{0,r), which may in the Lax or overcompressive case be 
decomposed as 

(5.7) Gx = Ex + Sx + Rx 
where, for y < 0: 

(5.8) ExMy):=X-^ J^ [cfj^ (:;.'g ) i^re(V<^-^^^^'/"^"^^ 



(5.9) 

Sx{x,t;y):= J] [^.,+ ]^*+^*-*g(-A/a;-+A^/3r/<-)x+(A/<-A^/?r/<-^). 



*+ , i2o,+ , , + 3. 



a;->0,a*+>0 
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for y < < X, 
(5.10) 

Sx{x,t;y) 



/ ^ K K 



"fc >o 



+ V [cJ''-ii?*-L*-*e^-^/"r+^'/5r/"r')^+(Var-A'/3r/<"')y 

/ -J L /c, — J J k 

for y < X <0, and 
(5.11) 

S\{x,t;y) := 

<-<o 

+ V [cJ'-ii?*-L*-*e(-^/"r+^'/5r/"r')^+(Var-A^/3r/<"')y 

/ ^ L fc, — J J A: 

/or X < y < 0, anc? i?A = Rf + Rf is a faster- decaying residual satisfying 

(a/ax)^(a/ay)"i?f = C(e-^l"-^l) + J2 e-^l"le(^/<*-^'^^^/<*')^ 

(5.12) «r>o 

X ('A"-iC(e^(^')^ - 1) + A"-iC(e^(^')" - 1) + C(A")) , 

(5.13) 

{d/dxr{d/dyrRf{x,t;y) = 

a*->0,a*+>0 

X ('A^+«C(e^(^')^ - 1) + A^+"C(e^(^')" - 1) + A"C(e-^l"l) + C(Ai+^; 

/or y < < X, 

(5.14) 

(a/ax)^(a/ai/)"i?f(x,t;y)= J] e(-^/< +^'^^ /< ')("-^) 

«r>o 

X (^A'^+"C(e^(^')^ - 1) + A^+"C(e^(^')" - 1) + A"C(e-^l"l) + C(A^+^) 

ap>0,a*-<0 

X ('A^+"C(e^(^')^ - 1) + A^+"C»(e^(^')" - 1) + A"C(e-^l"l) + C(A^+^) 
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for y < X < 0, and 

(5.15) 

{d/dx)'^{d/dy)''Rf{x,t;y)= ^ e^-^/< +^'^fe K ')(=^-y) 



g(-A/a*-+A2/3;-/a*-')x+(A/<--A'/3r/<"')y 



X ('A^+°C(e^(^')^ - 1) + A^+"C(e^(^')" - 1) + A"0(e-^l"l) + C(Ai+" 

+ E 

a* >0,a* <0 

X (A^+"C(e^(^')^ - 1) + A^+"C(e^(^')" - 1) + A"C(e-^l"l) + C(A^+" 

/or a:<y<0, q; + /3<1, tyit/i each 0{-) term separately analytic in X, and where 
I'^kal' CI = — ,0, +, are constants to be determined later. Symmetric bounds hold for 
y > 0. Similar, but more complicated formulae hold in the undercompressive case 
(see Remark 6.9, below). 

5.1. Normal modes (behavior in x). We begin by relating the normal 
modes of the variable coefficient eigenvalue equation (3.1) to those of the limiting, 
constant-coefficient equations (5.1). 

Lemma 5.3. Let (H0)-(H3) hold. Then, for A G -6(0, r), r sufficiently small, there 
exist solutions W- (x; A) of (3.1), C^ in x and analytic in X, satisfying 

^^^g^ W^{x;X) = V^^{x;X)e-f^ 

{d/dX)''V^^{x- A) = {d/dX)''V^^{X) + 0(e-^>l \Vj^{X)\), x ^ 0, 

for any k > and < 6 < 6, where 9 is the rate of decay given in (I.I4), fJ^j (A), 
V- (A) are as in Lemma 5.1 above, and 0{-) depends only on k, 9. 

Proof. With the exception of the regularity C^ asserted in x, this is a direct 
consequence of Lemma 5.1 together with Proposition 8.1 of Appendix A3.1 (the Gap 
Lemma) . The additional regularity of C^ vs. the C^ asserted in the Proposition is a 
consequence of the divergence form of the eigenvalue equation, and may be deduced 
as in the previous section by the change of coordinates (4.19). | 

The bases (j)- , ip- of Section 3.1 may evidently be chosen from among W- , 
yielding an analytic choice of bases in A, with the detailed description (5.16). It 
follows that the dual bases (j)- , V' • defined in Section 3.2 are also analytic in A and 
satisfy corresponding bounds with respect to the dual solutions (5.5)-(5.6). With 
this observation, we have immediately, from Definition (3.21) and Proposition 3.11: 

Corollary 5.4. Let (H0)-(H3) hold. Then, the Evans function -D(A) admits an 
analytic extension onto B{Q^r), for r sufficiently small, and the resolvent kernel 
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G\{x, y) admits a meromorphic extension, with an isolated pole of finite multiplicity 
at X = 0. 

Moreover, we have close control on the {x, y) behavior of Gx through the spectral 
decomposition formulae of Proposition 3.11. For slow, dual modes, the bounds 

(5.16) (in particular, the consequent bounds on first spatial derivatives) can be 
considerably sharpened, provided that we appropriately initialize our bases at A = 
0. This observation will be quite significant in the Lax or overcompressive case. 
Likewise, fast-decaying forward modes can be well-approximated near A = by their 
representatives at A = 0, using only the basic bounds (5.16). These two categories 
comprise the modes determining behavior of the Green's function to lowest order. 

Specifically, due to the special, partially conservative structure of the underlying 
evolution equations, the adjoint eigenvalue equation (3.31) at A = admits an 
n-dimensional subspace of constant solutions 

(5.17) W = {wl, wiY = (constant, 0)*; 

this is equivalent to the observation [L.2] that integral quantities in variable u 
are conserved under time evolution for relaxation systems of form (1.1). Thus, at 
A = 0, we may choose, by appropriate change of coordinates if necessary, that slow- 
decaying dual modes (p, and slow-growing dual modes t/j, be identically constant. 
Note that this does not interfere with our previous choice in Lemma 5.3, since that 
concerned only the choice of limiting solutions VF of the asymptotic, constant 
coefficient equations at x -^ ±oo, and not the particular representatives W, that 
approach them (which, in the case of slow modes, are specified only up to the 
addition of an arbitrary fast-decaying mode). 

Remark 5.5. The prescription of constant dual bases just described requires 
that we reverse our previous approach, choosing dual bases first using the Gap 
Lemma, then defining forward bases using duality. Alternatively, rewriting the 
forward eigenvalue equation at A = as 

where Z — i J := AW, we may observe that fast-decaying modes satisfy 

i.e., the linearization of the reduced traveling wave ode (1.11). By duality relation 
(3.33), requiring slow dual modes to be constant is equivalent to choosing fast- 
growing (forward) modes as well as fast-decaying modes from among the solutions 
of (5.19). (Recall: though fast-decaying modes are uniquely determined as a sub- 
space, fast-growing modes are only determined up to the addition of faster decaying 
modes.) 
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Lemma 5.6. Assuming (H0)-(H3), with the above choice of bases at X = 0, and 

J 



for A G -6(0, r) and r sufficiently small, slow modes W, satisfy 



(5.20) W^{y; A) = e-^(^)^y^(0) + XQf{y; A), 



where 
(5.21) 



\{d/dy)ef\ < C|e-^(^)^|(|A| + e-^l^l), 



9 > 0, as y ^ ±oo, and Vj = constant.^ Similarly, fast-decaying (forward) modes 
W^ satisfy 

(5.22) W^{x; A) = W^{x; 0) + Xef{x; A), 
where 

(5.23) |ef I, \{d/dx)ef\ < Ce-^l"l 

as X ^' ±oo, /or some 9 > 0. 

Proof. Applying the Gap Lemma to the augmented variables 

W±(!,;A) := ( ^^i- ) fc A) =: €-"«»¥* (j,; A) ( J^,) (!,;A) 



--'"^•M .,.^%.-)to^^) 



and 



•±(x;A):=(^i,^(x;A)=:e^W-V±(x;A) 



we obtain bounds 



e"^^^^^ + +' +' 1 (y; A) 



Wt(a:;A)=Vt(x;A)e-^^-^^)" 
(5.24) ' ' 

{d/dX)''Yf{x; A) = (a/aA)'=V±(A) + 0(e-^l"l|V±(A)|), x ^ 0, 



J\7y \/ /JV^' V IJ 



^In particular, this includes \{d/dy)W ■ (y;A)| < C|Ae^J | in place of the general spatial- 

derivative bound \{d/dy)W^iy;X)\ < C{\fiV^^\ + \{d/dy)V^^\)e^f^^''y ~ C|e^*^^^^|. 
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and 

(5 25^ 



9 > 0, analogous to (5.16), valid for A G -6(0, r), where 
and 

Via^ — I ^ 



By Taylor's Theorem with differential remainder, applied to V, we have 

(5.26) 

1 

2 



^(y. A) = e-^^ (^)^ ( V±(y; 0) + \{d / d\)N'r {y- 0) + -A2(a/aA)2v±(2/; A. 



for some A* on the ray from to A, where, recall, {d/dX)Y- {y; ■) and (9/9A)^V • {y; ■) 
are uniformly bounded in L°°[0, ±oo] for A G 5(0, r). Together with the choice 
V- {y, 0) = constant, this immediately gives the first (undifferentiated) bound in 
(5.21). 

Applying now the bound (5.25) with /c = 1, we may expand the second coordinate 
of (5.26) as 
(5.27) 



{d/dy)wHy. A) = e-^^ (^)^ - ^^^vHy■. 0) + V^ (y; 0) 



= e-'^^^)^ [-\{{d/d\)iif{Q)V^{Q) + O(e-^l^l)) + 0{\ 

and subtracting off the corresponding Taylor expansion 

{d/dy){e-^f^^^yV^{y,Q)) 

= ^±(A)e-^^^)^V^±(y,0) 

= e-^^^)^ (^_^±(0)lA±(y;0) - A(a/aA)^±(0)V^±(y;0) + ^(A^ 

= e-^^(^)^ {~\{d/d\)^ifm^{y- 0) + 0(A^ 

we obtain 

(5.29) {d/dy)Qf{y, A) = e"^^^^^)^ (AC(e-^l^l) + 0{\ 



(5.28) 
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as claimed. 

Similarly, we may obtain (5.22)-(5.23) by Taylor's Theorem with differential 
remainder applied to W = e'^^^'^ , and the Leibnitz calculation 

{d/dX)W = {d^i/d\)xe^''V + e^""" {d / d\)V, 

together with the observation that |xe~^l^l < Ce"^'^'/^. | 

Remark 5.7. For a more explicit derivation of (5.20)-(5.21) in the scalar case, 
see Lemma 2.2, (2.6)-(2.7) in [HoZ.2]. (Note: in [HoZ.2], the one- dimensional case 
corresponds to ^ = 0). 

This leaves only the problem of determining behavior in A through the study of 
coefficients M, d^. To this end, we make the following further observations in the 
Lax or overcompressive case, generalizing the corresponding observation of Lemma 
4.30, [Z.4], in the strictly parabolic case: 

Lemma 5.8. For transverse f7 7^ Q) Lax and overcompressive shocks, with the 
above- specified choice of basis at X = 0, fast-growing modes ijjj~, ijj~ are fast- 
decaying at —00, +00, respectively. Equivalently, fast-decaying modes i/j^', -^J are 
fast-growing at —00, +00 .■ i.e., the only bounded solutions of the adjoint eigenvalue 
equation are constant solutions. 

Proof. Noting that the manifold of solutions of the r-dimensional ODE (5.19) 
that decay at either x -^ +00 or x ^ —00 is by transversality, together with (1.15), 
exactly (i_|_ + (i_ — £ = r, we see that all solutions of this ODE in fact decay at least 
at one infinity. This implies the first assertion, by the alternative characterization 
of our bases described in Remark 5.5 above; the second follows by duality, (3.33). 
Since the manifold of fast-decaying dual modes is uniquely determined, independent 
of the choice of basis, this implies that nonconstant dual modes that are bounded at 
one infinity must blow up at the other, hence bounded solutions must be constant 
as claimed. | 

Remark 5.9. As noted in [LZ.2,ZH,Z.2,Z.4], the property of Lax and overcom- 
pressive shocks that the adjoint eigenvalue equation has only constant solutions 
has the interpretation that the only L^ time-invariants of the evolution of the one- 
dimensional linearized equations about {u,v){-) are given by conservation of mass, 
see discussion [LZ.2]. This distinguishes then from undercompressive shocks, which 
do have additional L^ time-invariants. The presence of additional time-invariants 
for undercompressive shocks has significant implications for their behavior, (see dis- 
cussions [LZ.2] Section 3 and [ZH], Section 10): in particular, the time-asymptotic 
location of a perturbed undercompressive shock (generically) evolves nonlinearly, 
and is not determinable by any linear functional of the initial perturbation. By 
contrast, the time-asymptotic location of a perturbed (stable) Lax or overcompres- 
sive shock may be determined by the mass of the initial perturbation alone. 
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5.2. Scattering coefRcients (behavior in A). We next turn to the estimation 
of scattering coefRcients Mjk, d-^^. Consider coefficient Mjk- Expanding (3.56) using 
Kramer's rule, and setting 2; = 0, we obtain 

(5.30) M% = D-'C%, 

where 



(5.31) C+:=(/,0)(|; !_;) 



vl>-' 



and a symmetric formula holds for C~ . Here, P^'^J denotes the adjugate matrix of 
a matrix P, i.e. the transposed matrix of minors. As the adjugate is polynomial in 
the entries of the original matrix, it is evident that |C^| is uniformly bounded and 
therefore 

(5.32) \M%\ < Ci\D-^\ < C2\-^ 

by ("D), where Ci, C2 > are uniform constants. 

However, the crude bound (5.32) hides considerable cancellation, a fact that will 
be crucial in our analysis. Relabel the {(/'^} so that, at A = 0, 

(5.33) ^i,_^^^^^- = {d/d5,)(^^s^, j = !:■■■ J. 

With convention (5.33), we have the sharpened bounds: 

Lemma 5.10. Let (H0)-(H4) and (V) (equivalently, (D1)-(D2)) hold, and let (p^ 
be labeled as in (5.33). Then, for |A| sufficiently small, there hold 



(5.34) \M%l\d%Ud%\ < C 



4- , , ^ , , -4- , \ ^ ^ for 7 = 1, ■ ■ • , £, 

^'^ ^'^ ^'^ "1 otherwise. 



where M"^ , d^ are as defined in Proposition 3.11. Moreover, 
(5.35) Residue a=o Af^_^_^^ ^^ = Residue a=o t^^fc, 

for l<j<e, all k. 

That is, blowup in Mjk occurs to order X^~^\D~^\ rather than |-D~^|, and, more 
importantly, only in (fast-decaying) stationary modes {d/d5j){u^ ,v^). Moreover, 
each stationary mode {d/dSj){u^ , v^) has consistent scattering coefficients with dual 
mode •0^ , to lowest order, across all of the various representations of Gx given in 
Proposition 3.11. 
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Proof. Formula (5.31) may be rewritten as 

(5.36) C+ = det 





jk 

from which we easily obtain the desired cancellation in M"*" = C^D~^. For example, 
for j > £, we have 

^t + Av?+, + 

. ^ c'+=det ; 

(5-37) I ■■■. ,■■■, ^t' + x<' + 



yielding \Mjk\ = \Cjk\\D\~^ < C as claimed, by elimination of £ zero-order terms, 
using linear dependency among fast modes at A = 0. For 1 < j < ^, there is 
only an £ — 1 order dependency, and we obtain instead the bound \Cjk\ < C|A|^~^, 
or \Mjk\ = \Cjk\\D\~^ < CIAI""*^. The bounds on \djk\ and \djk\ follow similarly; 
likewise, (5.35) follows easily from the observation that, since columns ip^ and 
V'iv-j+i agree, and V'7 is held fixed, the expansions of the various representations by 
Kramer's rule yield determinants which at A~^ order differ only by a transposition 
of columns. | 

In the Lax or overcompressive case, we can say a bit more: 

Lemma 5.11. Let (H0)-(H4) and (V) (equivalently, (D1)-(D2)) hold, with |A| 
sufficiently small. Then, for Lax and overcompressive shocks, with appropriate basis 
at X = (i.e. slow dual modes taken identically constant), there hold 

(5.38) \M,k\Adjk\Adjk\ < C 
if i^k is a fast mode, and 

(5.39) \M^kl\djkl\djk\<C\\\ 

if, additionally, ipj is a slow mode. 

As we shall see in the following section, this result has the consequence that only 
slow, constant dual modes play a role in long-time behavior of G. 

Proof. Transversality, 7 7^ 0, follows from (D2), so that the results of Lemma 5.8 
hold. We first establish the bound (5.38). By Lemma 5.10, we need only consider 

j = 1, . . . ,£, for which (pj = (p'J . By Lemma 5.8, all fast-growing modes ip^, ip^ 
lie in the fast-decaying manifolds Span ((/?J~, ■ ■ ■ , ^~_), Span ((/?^, ■ ■ ■ , iff ) at —00, 
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+00 respectively. It follows that in the righthand side of (5.36), there is at A = a 
linear dependency between columns 

(5.40) ftr--, v/-!, V'fc , Vj"+i, vt and (fj = v?+ 
i.e. an £-fold dependency among columns 

(5.41) ipf,---,ij~,---,ip+^andip];,---, ipj. 

It follows as in the proof of Lemma 5.10 that \Cjk\ < CX for A near zero, giving 
bound (5.38) for Mjk- If ^j is a slow mode, on the other hand, then the same 
argument shows that there is a linear dependency in columns 



(5-42) ifj,--- ,ipj_^,t(j,^,ipj^^,ip 



»+ 



and an {£ + l)-fold dependency in (5.41), since the omitted slow mode (p^ plays 



no role in either linear dependence; thus, we obtain the bound (5.39), instead. 



Analogous calculations yield the result for d.^, djk as well. 



5.3. Proof of Proposition 5.2. The proof of Proposition 5.2 is now just 
a matter of collecting the bounds of Lemmas 5.3-5.11, and substituting in the 
representations of Proposition 3.11. More precisely, approximating fast-decaying 
dual modes (f)'^ , . . . ,(f)f and (j)]^_i_^_i, . . . , (j)]^ by stationary modes {d/d5j){u^ , v^), 

and slow dual modes by e~^^^^^'^V{X) as described in Lemma 5.6, truncating (u(A) at 
second order, and keeping only the lowest order terms in the Laurent expansions of 
scattering coefficients M='=, d^ we obtain E\ and S\, respectively, as order A~^ and 
order one terms, except for negligible (9(e~^l^~^l) terms which we have accounted 
for in error term Rf^. 

Besides the latter, we have accounted in term R^ for: (i) the truncation er- 
rors involved in the aforementioned approximations (difference between stationary 
modes and actual fast modes, and between constant-coefficient approximant and 
actual slow dual modes, plus truncation errors in exponential rates); and, (ii) ne- 
glected fast-decaying forward/slow-decaying dual mode combinations (of order A, 
by Lemma 5.10). A delicate point is the fact that term (9(e~^l^~^l)(9(e~^l^l) arising 
in the dual truncation of E through the A(9(e~^l^l) portion of estimate (5.21) for 
{d/dy)Q- , absorbs in the (leading) (9(e~^l^~^l) term of (5.12). It was to obtain 
this key reduction that we carried out expansion (5.21) to such high order. 

In i?f , we have accounted for truncation errors in the slow/slow pairings approxi- 
mated by term S\ (difference between constant-coefficient approximants and actual 
slow forward modes, corresponding to an additional e"^'^' term in the truncation 
error factor, and between constant-coefficient approximants and actual slow dual 
modes, plus truncation errors in exponential rates). We have also accounted for 
the remaining, slow-decaying forward/fast-decaying dual modes, which according 
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to Lemma 5.11 have scattering coefficient of order A and tlierefore may be grouped 
witfi tfie difference between dual modes and tfieir constant-coefficient approximants 
(of smaller order by a factor of e"^'^') in the term 0{X). | 

5.4. Extended spectral theory. Finally, we cite without proof the extended 
spectral theory of [ZH] that we shall need in order to establish necessity of condition 
(V). In the case that (V) holds, this reduces to simple relations (5.35). 

Let C^p denote the space of C°° functions decaying exponentially in all deriva- 
tives at some sufficiently high rate. 

Definition 5.12. Let L be a linear ordinary differential operator with bounded, 
C°° coefficients (so that L : C^^ — ^ C^^ ), and let G\ denote the Green's function 
of L — XI. Further, let O be an open, simply connected domain intersecting the 
resolvent set of L, on which G\ has a (necessarily unique) meromorphic extension. 
Then, for Aq G O, we define the effective eigenprojection V\^^ : C^p -^ C°° by 

/ + CX) 
P\o{x,y)f{y)dy, 
-oo 

where 

(5.43) Pxo{x,y) = Resx^,Gx{x,y) 

and ResAo denotes residue at Xq. We will refer to Pxf^{x,y) as the projection 
kernel. Likewise, we define the effective eigenspace E^ (L) by 

E^^(L)= Range (Pao), 

and the effective point spectrum c' (L) of L in Q to be the set of X E fl such 
that dim ^'x,,{L) ^0. 

The above definition is the natural one from the point of view of the spectral 
resolution of the identity, / = f-p{L — XI)~^dX, hence also from the point of view 
of asymptotic behavior of solutions of PDE (see Section 6, below). Away from the 
essential spectrum of L, it corresponds with the usual definition [Kat,Y]. 

Definition 5.13. Let L, O, Aq be as above, and K be the order of the pole of 
(L — A/)~^ at Xq. For Aq G O, and k any integer, we define Qxo,k '■ C'^p -^ C°° 
by 



/-l-oo 
Qxo,k{^:y)f{y)dy, 
-oo 



where 

Q\o,k{x, y) = ResAo (A - Ao)''G'A(a;, y)- 

For < k < K , we define the effective eigenspace of ascent k by 

SAo,fc(^) = Range {Qxo,K-k) 
With the above definitions, we obtain the following, modified Fredholm Theory. 
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Proposition 5.14 [ZH]. Let L, \q, O he as in Definition 5.12, and K he the order 
of the pole of G\ at Aq . Then, 

(i) The operators V\g, Qa,),A: : C'^p -^ C°° are L-invariant, with 

(5.44) Qxo,k+i = {L-XoI)Qx„k = QxaAL-KI) 
for all k 7^ —1, and 

(5.45) Qx„k = {L - XoI)''Vxo 

for k > 0. 

(a) The effective eigenspace of ascent k satisfies 

(5.46) ^'x„kiL) = {L- Xom'x,,,+^{L). 
for all < k < K , with 

(5.47) {0} = E',^,oW C Ko,iiL) C ■ ■ ■ C KoM^) = KoW- 

Moreover, each containment in (5.47) is strict. 

(Hi) On V^^{C^p ), Vx„, Qx„,k all commute (k > 0), and Vxo is a projection. More 
generally, VxJ = f for any f e Ex,{L : C^^), henceEx.AL : C^^) C S'^^ ,(L) 
for allO<k< K. 

(iv) The multiplicity of the eigenvalue Aq, defined as dim E'^ (L), is finite and 
hounded hy Kn. Moreover, for all < k < K , 

(5.48) dim E',^,,(L) = dim E'^.^L*). 
Further, the projection kernel can he expanded as 

(5.49) PA„ = 5^V',(x)7r,(y), 

j 

where {i^j}, {ttj} are hases for T,'y^ [L], E^.(L*), respectively. 



loo 
^exp 



(v) (Restricted Fredholm alternative) For g E CI 

(L - Ao/)/ = g 
is soluhle in C°° if and soluhle in C^p only if Qxo,k-i 9 = 0, or equivalently 

(5.50) geE'^^^r 



-*\± 
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Proof: See [ZH]. | 

Remarks: For Aq in the resolvent set, Vx^ agrees with the standard definition, 
hence E^ ^ (L) agrees with the usual L^ eigenspace of generalized eigenfunctions of 
ascent < /c, for all p < cxd, since C^^ is dense in L^, p > 1, and Ti'^ j^(L) is closed. In 
the context of stability of traveling waves (more generally, whenever the coefficients 
of L exponentially approach constant values at ±00), E'^ ^(L) lies between the L^ 
subspace Eao,A:(-^) and the corresponding i^Loc subspace. 

The modification from standard Fredholm Theory is primarily that here conclu- 
sions (iii) and (v) apply on restricted domains. This restriction is clearly necessary, 
since the various expressions occurring in their statements would otherwise not be 
defined. In this regard, note that 'P\^^ is, strictly speaking, not a projection in 
the case that A G a^ssiL), since its domain does not then match its range. How- 
ever, it maintains the projective structure (5.49), somewhat justifying our abuse of 
notation. 

Now, suppose further (as holds in the case under consideration) that, on O: 
(i) L is a nondegenerate operator of some order s, i.e. the coefficient matrix of 
the principal part of the symbol is nonsingular for all x^; and, (ii) there exists an 
analytic choice of bases for stable/unstable manifolds at +00/— 00 of the associated 
eigenvalue equation. 

In this case, we can define an analytic Evans function -D(A) as in (3.21), which 
we assume does not vanish identically. It is easily verified that the order to which 
D vanishes at any Aq is independent of the choice of analytic bases $^. The 
following result, established in [ZH], generalizes to the extended spectral framework 
the standard result of Gardner and Jones [GJ.1-2] in the classical setting. 

Proposition 5.15 [ZH]. Let L, Aq be as above. Then, 

(i) dim E^ (L) is equal to the order d to which the Evans function Dl vanishes 
at Aq. 

(ii) Px^{L) = '^j '^j{x)TVj{y), where ipj and ttj are in E^^(L) and E^^(L*), 
respectively, with ascents summing to < K + 1, where K is the order of the pole of 
G\ at Aq. 

Remark 5.16. Evidently, the theory of this section goes through also, with 
suitable modifications regarding regularity, in the case (as for relaxation shocks) 
that the coefficients of L have finite regularity C^^^ , i.e. decay exponentially in up 
to s derivatives, in which case L : C^^p -^ C'exp • 

Section 6. Pointwise Green's function bounds. 



'^This assumption was made in [ZH] for convenience in calculation; the result holds in consid- 
erably more generality (for example for systems with real viscosity). 
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In this section, we prove the main results of the paper, estabhshing pointwise 
bounds on the Green's function G and consequently sufficiency and necessity of 
condition (V) for linearized orbital stability. 



6.1. Proof of Theorem 1.10. We will establish Theorem 1.10 as a corollary 
of the following, more detailed result: 



Proposition 6.1. Assuming (H0)-(H4) and the spectral stability criterion (V) 
(equivalently, (D1)-(D2)), the Green's function G{x,t;y) associated with (1.17) may 
he decomposed as in Theorem 1.10, where, for y <0: 
(6.1) 

n 

+ J2o(^{t + 1)- V2e-'7-+ + e-"!^!) (t + i)-i/2e-(x-y-ar*)VMt 
k=i 

a*~>0,a*~<0 



(6.2) 

J J 

R^x, t;y) = J2 0{e-^')6,.,^t{-y) + J2 0{e-^^\^-y\+'^: 
j=i i=i 

n 

-{x-y-a*-tf/Mt 



+ ^ O ("(t + l)-ie-''"^ + e-'^l"l) (t + l)-^/2g 
k=i 

a* >0,a* <0 

K, J 



3/2g-(x-a;- {t-\y/a-\)y/Mt^-rjx^ 



1 



al >0,a* + >0 
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and 

(6.3) 

J J 

Ry{x, t;y) = J2 0{e-^')6^--a^t{-y) + J2 0{e-^^\^-y\+'^) 

n 

fc=i 

+ E X{i<-ti>i,i}0((^ + i)-^/'e-(— r(^-i^K-i))VMt,-..- 

a*->0,a*+>0 

/or some i], M > 0, where cij are as in Theorem 1.10, x^ denotes the posi- 
tive/negative part of X, and indicator function Xs\a*-t\>\ \\ '^^ ^''^^ f^''^ \'^*k~^\ — Ivl 
and zero otherwise. Symmetric bounds hold for y > 0. 

Proof. Our starting point is the representation 

(6.4) G{x, t- y) = — P.V. / e^*GA(x, y)d\ 

valid by Proposition 2.3 for any tj sufficiently large. 

Case I. \x — y\/t large. We first treat the simpler case that \x — y\/t > S, S 
sufficiently large, for which G is known a priori to vanish, by the general property of 
finite propagation speed for systems of hyperbolic type. More precisely, G = for 
X ^ [zi{y, t), -2j(y, t)], where Zj{-, ■) are the standard hyperbolic characteristic paths 
defined in (1.36), ai < • ■ ■ < aj denoting the eigenvalues of A := df{u,v). This 
may be seen, for example, by standard energy estimates as in, e.g., [Ev]; however, 
it is instructive to recover this result by direct calculation. 

Suppose, then, that x ^ \z\,zj\. Fixing x, y, t, set A = 77 + i^. Applying 
Proposition 4.3, we obtain, for 77 > sufficiently large, the decomposition 

(6.5) G{x, t; y) = ^-P-V. / e^'Hx{x, y)d\ + ^-P-V. / e^*eA(x, y)d\ 

where 

(6.6) 

-E/=K+i«i(2/)"^e"-'^!' ^''''^^^'^^rj{x)Cj{x,y)l]{y) x > y, 



H\{x,y) 
and 



EjLittjd/) ^e ■^y^''''^^^'^^rj{x)Cj{x,y)l]{y) x<y, 



(6.7) Qxix, y) = X-^B{x, y; A) + A"^ {x - y)C{x, y; A) + X-^D{x, y; A) 
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with (j, B, C, D as defined in Proposition 4.3. 

For definiteness, take x > y. Then, the first term of (6.5) contributes to integral 
(6.4) the exphcitly evaluable quantity 

(6.8) 

-^P.V. / e^'Hx{x,y)d\= ^ iaj{y)-^F.y. e'^y'-^y^^^''''^'^^'^') dn 

27n Jr]-ioo j=K+l ^ J -oo J 

xe^(*-/;(Va.W)^^)r,.(;,)C,.(^,^)/5(^) 

= E (a,(y)-i5(t- r(l/a,(^))rf^))r,(x)C(a:,y)/5(y) 
j=K+i \ Jy / 

J 

j=K+l 

= H{x,t;y), 

where dj, Q, and H are as defined in Theorems 1.10 and 1.11. Note: in the final 
equality of (6.8) we have used the general fact that 

(6.9) fy{x,y,t)6{f{x,y,t)) = 6h{x,t){y). 

provided fy ^ 0, where h{x,t) is defined by f{x,h{x,t),t) = 0. This term clearly 
vanishes for x outside [zi{y,t),zj{y,t)], as claimed. 

Similar calculations show that the second, error term in (6.5) also vanishes. For 
example, the term e^*\~^B{x,y\ A) contributes 
(6.10) 

— E fP-V- r"(^ + ^0-'e^^(*-^-^^'/"^(^)^'^)rfe')e''(*-^.^(i/"^(^»^^)6+(x,y). 

The factor 

P.V. j '" {r] + iO~^e^^'-^y^^''''^'^^'^') d^, 

J — oo 

though not absolutely convergent, is integrable and uniformly bounded as a principal 
value integral, for all ?] G M, by explicit computation. On the other hand, 

(6.11) g^(*-/y''(l/«j(2))ci2) <:_ (J^-Vd{x,[z^{y,t),z,,{y,t)])/ minj^:,\aj{x)\ ^q 

as ?7 ^ +00, for each K + 1 < j < J , since Oj > on this range. Thus, taking 
?7 ^ oo, we find that the product (6.10) goes to zero, giving the result. 

Similarly, the contributions of terms e^^C{x,y; X) and e^^D{x,y; X) split into 
the product of a convergent, uniformly bounded integral in ^, a (constant) factor 
depending only on (x,y), and a factor q;(x, y, t, ?]) going to zero as ry -^ at rate 
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(6.11). Thus, taking i] -^ +00, we find that each of these terms vanishes, and thus 
G vanishes also, as claimed. 

Case II. \x — y\/t bounded. We now turn to the critical case that \x — y\/t < S 
for some fixed S. In this regime, note that any contribution of order e^*, 9 > 0, 
may be absorbed in the residual (error) term R (resp. R^, Ry)', we shall use this 
observation repeatedly. 

Decomposition of the contour. We begin by converting contour integral (6.4) 
into a more convenient form decomposing high, intermediate, and low frequency 
contributions. 

Observation 6.2. Assuming (H0)-(H4), there holds the representation 
G{x,t;y) = Ia + Ib + IIa + Ih 



(6.12) 



:= 


1 

2m 


:P.V. 


J -q — ioo 


,{x,y)dX 












+ 


1 

2711 


P.V. 


' r-vi 


-iR 

+ 

-ioo J 


'-m+iR J 


gAt 


(Gx- 


H> 


J(x, 


,y)d\ 


+ 


1 

27ri 


!/• 


'Gx{x, 


y)d\- 


1 /--^^ 

27ri y_^^. 


_+iR 
-iR 


e^'HxiX; 


,y)dX, 



(6.13) r := [-?7i -iR.rj- iR] U[r}-iR,r} + iR] U [t] + iR, -r]i + iR], 

for any r] > such that (6.4) holds, R sufficiently large, and —rji < as in (3.3). 
(Note that we have not here assumed (V)). 

Proof. We first observe that, by Proposition 4.3, L has no spectrum on the 
portion of F := {A : Re A > —r]i < 0} lying outside of the rectangle 

(6.14) 7^ := {A : -r]i < Re A < 77, -i? < Im A < i?} 

for ?7 > 0, i? > sufficiently large, hence G\ is analytic on this region. Since, also, 
H\ is analytic on the whole complex plane, contours involving either one of these 
contributions may be arbitrarily deformed within r\7^ without affecting the result, 
by Cauchy's Theorem. 
Recalling, further, that 

\Gx-H,^\ = 0{\-^) 

by Proposition 4.3, we find that 

—P.V. / e'^\Gx-Hx){x,y)d\ 
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may be deformed to 

^:P.V. / e''\Gx-Hx){x,y)d\ 

ZTV'l 

Noting, finally, that 



27ri Jd{A\n) 



+ -^ / e^'Hx{x,y)dX = 0, 
27r« JdTz, 



by Cauchy's Theorem, we obtain, finally, 
(6.15) 

— P.V. / e^^Gx{x,y)dX = 



2m 

-- / e^'Hx{x,y)d\+—-.'PN. / e^\Gx- Hx){x,y)d\ 

+ 7^ / e^'Hx{x,y)d\ 
27rz J^T^ 

from which the result follows by combining the second and third contour integrals 
along their common edges F. | 

Remark 6.3. We have already shown in (6.8) that 2^P-V. P_l'^ e'^^H\{x^ y)d\ - 
H{x, t; y) independent of ry. Thus, we may deform the the entire contour integral 



—P.V./ e^*Gx{x,y)d\ 

"'' J n—ioo 



to 



^TTZ .'rj—ioo 



T^P.V. / e^'Gx{x,y)dX 

27r« Jd{A\n) 



to obtain (6.12) by a different route. 

Observation 6.4. Assuming (V) together with (H0)-(H4), we may replace (6.12) 
by 

G{x, t- y) = Ia + h + Ih + Ih + ///, 

where la, lb and Ih are as in (6.12), and 

1 / i'-ili-ir/2 p—Vl+iR^ 

Ih:=—i + e^'GX{x,y)dxi, 

^^ \J-r]i-iR J-r]T_+ir/2j 
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(6.16) /// := 7^ / e'-'Gxix, y)dX, 

and 

(6.17) f := [-rji - zr/2, r] - ir/2] U [ij - zr/2, r] + ir/2] U [r] + zr/2, -rji + zr/2], 

for any t], r > 0, and rji sufficiently small with respect to r. 

Proof. By assumption ("D), L has no spectrum on the region between contour 
r and the union of contour F and the contour of term I la, hence Gx is analytic on 
that region, and 

IIa = IIa+in 

by Cauchy's Theorem, giving the result. | 

Using the final decomposition (6.16), we shall estimate in turn the high frequency 
contributions /„ and lb, the intermediate frequency contributions I la and lib, and 
the low frequency contributions ///. 

High frequency contribution. We first carry out the straightforward esti- 
matation of the high-frequency terms /„ and /&. The principal term /„ has already 
been computed in (6.8) to be H{x,t;y). Likewise, calculations similar to those of 
(6.10)-(6.11) show that the error term 



1 



rji—R f — Tji+ioo 



/7r« \ .l-rii-ioo J-r]i+iR J 

= — P.V. / / e^*(A-5(x,,;A) 

+ X~^{x - y)C{x, y; A) + X~'^D{x, y; X))dX 

is time-exponentially small. 

For example, forx>y the term e^*A~^i?(x, y; A) contributes 

y — P.V. ( /""" r^]{-rji+i^)-h<'-h^^^/'^^^^^^''^) d^ 

(6.18) ,=K+i 2^^ V-^— -^R J 

xe-''i(*-/;(i/"^(^»^06+(x,y), 
where 

(6.19) —P.V. if [ ^] {rj + i^yh'^i'-i:^^/''^^'')')'^') dC < oo, 
2m \J-^ Jr / 
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and 

(6.20) e'''-fy^^/''^^'^^'^%{x,y) < c^e"' fy^^'^J^'^^'^'-'^^'^-y^ < C2, 

for r]i sufficiently small. This may be absorbed in the first term of R, (4.36). 

Likewise, the contributions of terms e^^C{x,y; X) and e^^D{x,y; X) split into 
the product of a convergent, uniformly bounded integral in ^, a bounded factor 
analogous to (6.20), and the factor e~^^^, giving the result. 

Derivative bounds. Derivatives {d/dx)Ib (resp. {d/dy)Ib) may be treated in 
identical fashion using (4.13) (resp. (4.14), to show that they are absorbable in 
the estimates given for Rj^ (resp. Ry). We point out that the integral arising from 
term B^ (resp. By) of (4.13) (resp. (4.14)) corresponds to the first, delta-function 
term in (6.2) (resp. (6.3)), while the integral arising from term {x — y)C^ (resp. 
{x — y)Cy) vanishes, as the product of (x — y) and a delta function dz{y) with z ^ x. 

Intermediate frequency contribution. Error term //^ is time-exponentially 
small for r\\ sufficiently small, by the same calculation as in (6.18)-(6.20), hence 
negligible. Likewise, term Ila by the basic estimate (3.20) is seen to be time- 
exponentially small of order e"''^* for any 771 > sufficiently small that the associ- 
ated contour lies in the resolvent set of L. 

Low frequency contribution. It remains to estimate the low-frequency term 
///, which is of essentially the same form as the low- frequency contribution analyzed 
in [ZH,Z.2,Z.4] in the strictly parabolic case, in that the contour is the same and 
the resolvent kernel G\ satisfies identical bounds in this regime. Thus, we may 
conclude from these previous analyses that /// gives contribution E + S -V R^ as 
claimed, exactly as in the strictly parabolic case. For completeness, we indicate the 
main features of the argument here. 

Case t < 1. First observe that estimates in the short-time regime t < 1 are 
trivial, since then \e^^G\{x,y)\ is uniformly bounded on the compact set F, and 
we have |G(a;,t;y)| < C < e~ for ^ > sufficiently small. But, likewise, \E\ 
and 5" are uniformly bounded in this regime, hence time-exponentially decaying. 
As observed previously, all such terms are negligible, being absorbable in the error 
term R. Thus, we may add E ^ S and subtract G to obtain the result 

Case t > 1. Next, consider the critical (long-time) regime t > 1. For definiteness, 
take y < a; < 0; the other two cases are similar. Decomposing 

(6.21) 

G{x, t;y) = —-(l e^'Exix, y)dX + 7-- * e^'Sx{x, y)dX + ^—. I e^'Rx{x, y)dX, 
Ztti Jy Ztti Jy Ztti Jy 

with E\, S\, and Rx as defined in Proposition 5.2, we consider in turn each of the 
three terms on the righthand side. 

Ex term. Let us first consider the dominant term 

(6.22) ^^j_e^'Ex{x,y)dX, 
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which, by (5.8), is given by 

(«-23) El<-l|-(''(t))^rV(x,i;.), 



a, >0 



where 

(6.24) ak(x,t;y) := — i A-^e^*e^^/<"-^'^rK"')y rfA. 

Using Cauchy's Theorem, we may move the contour F to obtain 

(6.25) 

1 /•+^/2 _ 2 - -3 

1 / /•-^T-/2 /.-?7i+ir/2\ 



2m \^J_^^_i^l2 Jir/2 

+ ^Residue ;,^o e^*e(^/<~-^'^r/<"')y, 

or, rearranging and evaluating the final, residue term: 
(6.26) 

/ 1 /• + 00 _ 2 - -3 1 \ 

ak{x,t;y)= (—P.y. (z^)-ie^«(*+^/"fe ^e^ ('^^ /"^ ^^ c/^ + - ) 

-I / /•-r/2 r + o^X n ■, 



27r \ ./_oo J^/2 



2TXi \ J_^^_ij./2 Jir/2 

The first term in (6.26) may be explicitly evaluated to give 

(6.27) errfn ( ^ + ^^ ^ 

\^4(3r\y/al- 

where 

1 r^ 

(6.28) errfn (z) := — / e"^'rf2/, 

27r y_oo 

whereas the second and third terms are clearly time-exponentially small for t < C\y\ 
and r]i sufficiently small relative to r (see detailed discussion of similar calculations 
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below, under Rx term). In the trivial case t > C\y\, C > sufficiently large, we can 
simply move the contour to [—r]i — zr/2, —r]i + ir/2] to obtain (complete) residue 1 
plus a time-exponentially small error corresponding to the shifted contour integral, 
which result again agrees with (6.27) up to a time-exponentially small error. 
Expression (6.27) may be rewritten as 



(6.29) errfn 

plus error 



y + a-l t 



I "^ — -A I / I ^ — -i- 

errfn I — , — errfn ' 



errfn' I ^±£ J (-^{y + a*,- 1)^4(31- 1)-'/^^ 



(6.30) 

~ errfn ' I — 

'4(3rt 

= 0(t-ie(^+<"*)'/^*), 
for M > sufficiently large, and similarly for x- and y- derivatives. Multiplying by 

J,o 1 9 fu^{x)\ ^_t 



we find that term (6.29) gives contribution 

(6.31) [cf_] I- ( "'("=; ) Lr'-rfn (^M±£L 

whereas term (6.30) gives a contribution absorbable in R (resp. R^, Ry)] see Remark 
6.5 below for detailed discussion of a similar calculation. 
Finally, observing that 

is time-exponentially small for t > 1, since a^~ > 0, y < 0, and {d/d6j){u^ ,v^) < 
Ce~^l^', 6* > 0, we may subtract and add this term to (6.31) to obtain a total of 
E{xjt;y) plus terms absorbable in R (resp. Rn^, Ry). 
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Remark 6.5. There is a slight discrepancy between the form of E used here, 
and the form 



errfn 



(6.33) 



E K 



n d fu^{x) 




errfn 



-X + y + al t 



^prt 



— errfn 



given in [Z.2]; however, this is only cosmetic, since the two expressions are equal 
modulo terms absorbable in the error terms R (resp. R^, Ry)- To see this, first 
note that the difference is time-exponentially small outside the critical region 

X e [-Clog(t + 1), +Clog(t + 1)], 

for C > sufficiently large. Then, estimating 



errfn 



—X + y ± al t 



^prt 



errfn 



y±al t 



u^rt 



< \x\ 



errfn ' 



-X + y ± at t 



< \x\t~^/^e'^''~y^<'^^^/^^*^~\ 

and recalling that {d/d6j){u^, v^) = (9(e~^'^l), we find that the error is bounded by 

(6.34) C|x|e-^l^lt-V2g(.-j,-<-t)V4/3r* < C^^-B\^\/2^-l/2^i.-y-al~tf/4(3:-t^ 

and thus absorbable in R. The calculation for the y-derivative is similar. We have 
chosen to present the result in the form used here, since it is more straightforward 
to prove, and is more convenient for the final stability analysis (in fact, this was 
the form that was actually used in the stability analysis of [Z.2], as well). It was 
presented as (6.33) in [Z.2] in order to capture the parabolic property of almost- 
finite propagation speed; here, we have dealt with propagation speed separately, 
which seems the preferable approach. 

Remark 6.6. As described in [Z.2], the time-exponentially small terms (6.32) 
that we have (artificially) included in E are useful for linear and nonlinear stability 



CORRADO MASCIA AND KEVIN ZUMBRUN 73 

analyses, serving to regularize behavior at t = 0. Note, also, that the resulting 
"balanced" form of E is the one that appears naturally in the explicit Green's 
function formula for viscous shock solutions of Burgers equation; see equation (2.13), 
[Z.2], and surrounding example. 

Sx term. Next, consider the second-order term 

(6.35) J- I e^'Sx{x,y)dX, 

Ztci Jy 

which, by (5.10), is given by 

(6.36) ^ i?*"L*" afc(x,t;y)+ ^ [4'I]i?*"L*" ^^^(x, t; y) 



a* >0,a. <0 



where 



(6.37) ak{x,t;y):=— (I e^'e^-^/< +^'^^ /< ')("-^) rfA. 

and 

(6.38) ajk{x,t;y) := — I e^'e^'^/^r+^'^r /-r')-+(^/<~ -^'^l~ /<~')y dX. 

2m Jf 

Similarly as in the treatment of the E\ term, just above, by deforming the contour 
f to 

(6.39) r' := [-?7i - ir/2, -ir/2\ U [-zr/2, +zr/2] U [+zr/2, -r]i + zr/2], 

these may be transformed modulo time-exponentially decaying terms to the ele- 
mentary Fourier integrals 

1 /• + 00 _ _ _3 

PV / e»C(*-(^-y)/"fc )eC (-/3fc /«*, )(^-y) c^^ 

(6.40) 27r ■ ■ J_^ 

= (47r/3-t)-i/2e-("-^-<"*)'/4^r* 



and 



^P,V, 



+ 00 



^^i{{t-\y/al-\)-./a-)^-e[{|i^ /ar')\y\HPr /-r')\-\) d^ 

(6.41) 27r J_^ 

= (47r;g*-t)-v2e-(--;.")'/4/3;."* 



respectively, where (3jk—~^ and z*/^ are as defined in (1.38) and (1.37). These 
correspond to the first and third terms in expansion (1.35), the latter of which has 
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an additional factor e~^/(e^ + e~^). Noting that the second and fourth terms of 
(1.35) are time-exponentially small for t > 1, y < x < 0, and that 

for some ^ > 0, so is absorbable in error term R, we find that the total contribution 
of this term, modulo terms absorbable in R, is 5". 

Rx term. Finally, we briefly discuss the estimation of error term 

(6.42) ^ (f e^'Rx{x,y)dX, 

which decomposes into the sum of integrals involving the various terms of R^ and 
i?f given in (5.12) and (5.14). Since each of these are separately analytic, they may 
be split up and estimated sharply via the Riemann saddlepoint method (method 
of steepest descent), as described at great length in [ZH,HoZ.2]. That is, for each 
summand a\{x^ y) ~ e^*^'^-'^"'"'''*^'^-*^ in R\ we deform the contour F to a new, contour 
in A that is a mini-max contour for the modulus 

m^{x,y,\) := \e^'ax{x,y)\ = e^' ^*+^<^ ^"+^^ T^ 

passing through an appropriate saddlepoint/critical point of ma^x, y, ■): necessarily 
lying on the real axis, by the underlying complex symmetry resulting from reality 
of operator L. 

Since terms of each type appearing in R have been sharply estimated in [ZH], 
we shall omit the details, only describing two sample calculations to illustrate the 
method: 

Example 1: e~^l^~^l. Contour integrals of form 

(6.43) -^ /e^*e-^l"-^'rfA, 

27rz Jp 

arising through the pairing of fast forward and fast dual modes, may be deformed 
to 

1 ^— 77i+ir/2 

(6.44) -— / e^*e-^l"-^lc^A 

and estimated as (9(e~^^*~^'^~^'), a negligible, time-exponentially decaying contri- 
bution. I 

Example 2: ye^~^f°'k +^ Pk /"-k )(^-y). Contour integrals of form 

(6.45) — /e^*A'^e(-^/<~+^'^rK"')(^-y)rfA, 
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a^~ > 0, arising through the pairing of slow forward and slow dual modes, may be 
deformed to contour 

(6.46) r' := [-7]! - zr/2, 7]^ - ir/2] U [77* - zr/2, 77* + ir/2] U [77* + zr/2, -771 + zr/2], 

where saddlepoint 77* is defined as 



s if|«|<e 
p I p I — 

±£ if ^ ^ £, 



(6.47) r],{x,y,t):= 

with 

(6.48) S:=^^i^, p := ^£<i-i^ > 0. 

2t (aj^t 

and ?7i, £ > are chosen sufficiently small with respect to r, to yield, modulo 
time-exponentially decaying terms, the estimate 

(6.49) 

/+00 
-00 

if 1-1 < £, and 



-00 



■(r+l)/2g-r;i^ 



if | — I > £, ^ > 0. In either case, the main contribution lies along the central portion 
[?7* — ir/2, ?7h< + ir/2] of contour V . 

To see this, note that |x — y| and t are comparable when \a/p\ is bounded, whence 
the evident spatial decay e~^^ \^-y\ of the integrand along contour A = 77* -|-i^ may 
be converted to the e"'' * decay displayed in (6.49)-(6.50); likewise, temporal growth 
in e'^* on the horizontal portions of the contour, of order < e'^l'^^i"'"''^!-**, is dominated 
by the factor e~^^ = e~^'^ '■*, provided [771 1 + |£| is sufficiently small with respect 
to r^. For \a/p\ large, on the other hand, 77* is uniformly negative, and also \x — y\ 
is negligible with respect to t, whence the estimate (6.50) holds trivially. 

We point out that 77^, is easily determined, as the minimal point on the real axis 
of the quadratic function 

(6.51) f,,y4X) := Xt{-X/ar + X'pr/ar'){x - y), 

the argument of the integrand of (6.45). | 

Other terms may be treated similarly: All "constant-coefficient" terms 4'k4'k* or 
il^k^l^k* are of either the form treated in Example 1 (fast modes) or in example 2 
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(slow modes). Scattering pairs involving slow forward and slow dual modes from 
different families (i.e., terms with coefficients Mjfc, (i-^) may be treated similarly 
as in Example 2. Scattering pairs involving two fast modes are of the form already 
treated in Example 1, since both modes of a scattering pair are decaying. Scattering 
pairs involving one fast and one slow mode may be factored as the product of a 
term of the form treated in example 2 and a term that is uniformly exponentially 
decaying in either x or y; factoring out the exponentially decay, we may treat such 
terms as in Example 2. 

Remark 6.7. In [ZH,HoZ.2], the saddlepoint contours followed were hyperbolae 
through ?7^, rather than the simpler, vertical contours used here. Any reasonable 
contour through the saddlepoint will suffice, yielding the same estimates. 

Scattering coefficients. The relation (1.39) may now be deduced, a posteriori, 
from conservation of mass in the u variable of the linearized flow (1.17). For, by 
inspection, all terms save E and S in the decomposition (1.32) of G decay in L^, 
whence these terms must, time-asymptotically, carry exactly the mass in u of the 
initial perturbation dy{x)In 



I.e., 



/ + 00 
{E{x,t;y) + S{x,t;y)) dx = I^+r. 
-oo 

Taking y < for definiteness, and right-multiplying both sides of (6.52) by (r^~*, 0)*, 
we thus obtain the result from definitions (1.34)-(1.35). 
Moreover, rewriting (1.39) as 

/ [cl--] \ 



(6.53) 



1 5 



■ • • ' <-i_ ' ^u ' • • • ' ^n^' mi,..., mi) 



J+'+l 



.'^.+1 






1,0 



£,0 



k ' 



V [4':-] / 



where 

(6.54) 



mj :-- 



+ 00 



{d/dSj)u^{x) dx, 



we find that, under assumption (V), it uniquely determines the scattering coeffi- 
cients [c;^'_]. For, the determinant 



(6.55) 



det {rl' 



".•••. <-i_ ^r*^^..., <+, mi,..., me) 
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of the matrix on the righthand side of (6.53) is exactly the inviscid stabihty coef- 
ficient A defined in [ZS,Z.4], hence is nonvanishing by the equivalence of (T>2) and 
(D2) (recall discussion of Section 1.2, just below Theorem 1.4). Note that (6.55) 
reduces to the Liu-Majda determinant (1.21) in the Lax case, ^ = 1, with u^{x) 
parametrized as u^{x) := u{x — 5ix), for which mi reduces to [u\. 
Relation (1.40) follows from the observation that 



P\{x,y) := Residue \=oGx{x,y) 

I 






with TTj defined by the first, or second expressions appearing in (1.40), according 
as y is less than or equal to/greater than or equal to zero. For, the extended 
spectral theory of Section 5.4 then implies that tij are the left effective eigenfunctions 
associated with right eigenfunctions {d/ddj){u^ ^v^). Alternatively, (1.40) may be 
deduced by linear-algebraic manipulation directly from (6.53) and its counterpart 
for y > (the same identity with k, — everywhere replaced by /c, +). 

Remark 6.8. As pointed out in [LZ.2], the left eigenfunctions (constant vectors) 
TTj = (ttj^^j, 0) may themselves be determined by the relation 

(6.56) ^*, Jrp , • • • , <!,_ , r*+, . . . , <+, mi, . . . , m,) = (0, ej), 

from which (1.40) is easily deduced, where Cj here denotes the jth standard basis 
element in M^. This may be deduced a priori through the extended spectral theory 
of Section 5.4; it would interesting also to derive (6.53) a priori, directly through 
the relations (3.54)-(3.55). 



Remark 6.9 (The undercompressive case). In the undercompressive case, 
the result of Lemma 5.8 is false, and consequently the estimates of Lemma 5.11 
do not hold. This fact has the implication that shock dynamics are not governed 
solely by conservation of mass, as in the Lax or overcompressive case, but by more 
complicated dynamics of front interaction; for related discussion, see [LZ.2,Z.2]. At 
the level of Proposition 5.2, it means that the simple representations of E\ and S\ 
in terms of slow dual modes alone (corresponding to equilibrium characteristics that 
are incoming to the shock) are no longer valid in the undercompressive case, and 
there appear new terms involving rapidly decaying dual modes ~ e"^'^' related to 
inner layer dynamics. Though precise estimates can nonetheless be carried out, we 
have not found a similarly compact representation of the resulting bounds as that of 
the Lax/overcompressive case, and so we shall not state them here. We mention only 
that this rapid variation in the y-coordinate precludes the L^ stability arguments 
used here and in [Z.2], requiring instead detailed pointwise bounds as in [HZ.2,Z.5]. 
See [LZ.1-2,ZH,Z.2,Z.5] for further discussion of this interesting case. 
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6.2. Inner layer dynamics. Similarly as in [ZH], we now investigate dynamics 
of the inner shock layer, in the case that {T>) does not necessarily hold, in the process 
establishing the necessity of {V) for linearized orbital stability. 

Proposition 6.10. Let (H0)-(H4) hold. Then, there exists rj > such that, for x, 
y restricted to any bounded set, and t sufficiently large, 

(6.57) G{x,y;t)= Yl e^'Y.^^{L-\lfPx{x,y) + 0{e-'^'), 

Aeo-;(L)n{i?e(A)>o} fc>o 

where P\{x,y) is the effective projection kernel described in Definition 5.12, and 
o'_ (L) the effective -point spectrum. 

Proof. Similarly as in the proof of Proposition 6.1, decompose G into terms /„, 
/fe, //a, and Ilh of (6.12). Then, the same argument (Case II. \x — y\/t bounded) 
yields that la = H{x, t; y), while /;, and //;, are time-exponentially small. Observing 
that H — ior X, y bounded, and t sufficiently large, we have reduced the problem 
to the study of 

(6.58) IIa:=:^ (f e^'Gx{x,y)dX, 

where, recall, 

(6.59) r := [-?7i -iR.T]- iR] U[rj-iR,rj + iR] U [r] + iR, -r]i + iR]. 

Note that the high-frequency bounds of Proposition 4.3 imply that L has no point 
spectrum in O outside of the rectangle TZ enclosed by F U [—rji — iR, —rj + iR] . 
Choosing r]i sufficiently small, therefore, we may ensure that no effective point 
spectrum lies within the strip —rji < Re A < 0, by compactness of TZ together with 
the fact that effective eigenvalues are isolated from one another, as zeroes of the 
analytic Evans function. 

Recall that Gx is meromorphic on O. Thus, we may express (6.58), using 
Cauchy's Theorem, as 

(6.60) — - / e^^ G x{x,y)dX + Residue xene^^Gx{x,t]y). 



27ri j_^^_ipi 
By Definition 5.12 and Proposition 5.14, 
ResAe{Re \>o}e^^Gx{x,y) 

J2 e^^'Resxoe^^-^°'>'Gx{x,y) 

Ao6cr;(L)n{Re A>0} 

E ^^°* J](tV^!)ResA„(A - Xo)''Gx{x, y) 

(6.61) Ao6(j;(L)n{Re A>0} A:>0 

E e^^'Y.^^/k\)Qx,A^.y) 

Ao6cr;,(L)n{Re A>0} A:>0 

E e^^^'Y.^^/k\){L-\olfPx,{x,y)). 

Aoecr;,(L)n{Re A>0} A:>0 
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On the other hand, for x, y bounded and t sufficiently large, t dominates \x\ and 
\y\ and we obtain from Propositions 3.4 and 5.2 that 

\Gx{x,t;y)\<C 

for all A G O, for r]i sufficiently small, hence 

-r]x+iR 
-r}i_—iR 

Combining (6.60), (6.61), and (6.62), we obtain the result. 



(6.62) —[ ' e^^Gx{x,y)d\<2CRe-'^^K 

2iii J_„ ifi 



Corollary 6.11. Let (H0)-(H4) hold. Then, (V) is necessary for linearized orbital 
stability with respect to compactly supported initial data, as measured in any L^ 
norm, 1 < p < oo. 

Proof: From (6.57), we find that {u,v){-) is linearly orbitally stable only if 
Px = for aU A G {Re A > 0} \ {0} and Range Vq = Span {{d/d6j){u^ ,v^}. By 
Propositions 5.14 and 5.15, this is equivalent to (V). | 



Section 7. Stability. 

We now establish our remaining results on linearized and nonlinear stability. 

7.1. Linearized stability. Linearized orbital stability follows immediately 
from the pointwise bounds we have established. Similarly as in [Z.2], define the 

linear instantaneous projection: 

/+00 
E{x,t;y)Uo{y)dy 

-^(t)(a/a5)(^,g 



where Uq is the initial data for the linearized evolution equation (1.17). The coor- 
dinates 5 G M^^^ may be expressed, alternatively, as 

K^) = / e{y,t)vo{y)dy, 



where 

(7.2) E{x,t;y)=:{d/d6)(fs[l] ] e{y,t), 
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I.e., 

(7.3) e(2/, t) := ^ ( errfn ( ^±£ ] - errfn ( ^^£ 1 1 L*" 



4/^rv V\/4/?r^ 



for y < 0, and symmetrically for y > 0. 
Then, the solution U of (1.17) satisfies 



/-too 
{H + G){x,t;0)Uo{y)dy, 
-oo 

where 

(7.5) G:=S + R 

is the regular part and H the singular part of the time-decaying portion of the 
Green's function G. 

Lemina 7.1. Under assumptions (H0)-(H4) and (V), there hold: 

/ + CX3 
G(-,t;2/)/(y)dyU.<C(l + t)-^(i-i/'^)|/U., 
-oo 

/ + 00 
G{-,t-ymirYf{y)dy\Lv<C{l + t)-'^^^-^'^)-y^\f\L., 
-oo 

Gy{;t-y)f{y)dy\L.<C{l+t)-^^^-^'^^-^'^\f\L.+Ce-^'\f\L.. 

-oo 

and 

/ + 00 
H{-,t-y)f{y)dy\Lv<Ce-'^'\f\Lv, 
-oo 

/or a// t > 0, some C, r] > 0, for any 1 < r < p and f E L'^ (resp. L'p), where 
l/r + l/q = 1 + 1/p. 

Proof. Bounds (7.6)-(7.8) follow by the Hausdorff- Young inequality together 
with bounds (1.35) and (6.1)-(6.3), precisely as in [Z.2]. Bound (7.9) follows by 
direct computation and the fact that particle paths Zj{y,t) satisfy uniform bounds 

l/C<\{d/dy)zj\<C, 

for all y, t, by the fact that characteristic speeds aj{x) converge exponentially as 
X -^ ±oo to constant states. ■ 
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Corollary 7.2. Let (H0)-(H4) hold. Then, (V) is sufficient for L^ n Lp ^ Lp 
linearized orbital stability, for any p > 1. More precisely, for initial data Uq G 
L^nLP, the solution U{x,t) of the linearized evolution equations (1.17) satisfies the 
linear decay bounds: 

(7.10) \U{;t) - ip{;t)\L. < C{l+t)-^^'-'/P\\Uo\L^ + \Uo\Lr'). 

Proof. Immediate, from (7.4) and bounds (7.6) and (7.9), with q = p- § 

Proof of Theorem 1.4. Theorem 1.4 now follows immediately from Proposition 
6.11 together with Corollary 7.2. | 

7.2. Nonlinear stability (Proof of Theorem 1.13). We conclude by estab- 
lishing nonlinear stability of relaxation profiles for simultaneously symmetrizable 
and discrete kinetic models, as described in Theorem 1.13. 

7.2.1. Simultaneously symmetrizable models. We first establish the result 
of Theorem 1.13 for relaxation profiles of simultaneously symmetrizable relaxation 
models, under the weak shock assumption 

(7.11) \iu,vy\L^<e, 

e > sufficiently small with respect to the parameters of system (1.1). Before 
beginning, we recall some needed preliminaries regarding the class of simultaneously 
symmetrizable systems: 

Definition 7.3. A relaxation model (1.1) is called simultaneously symmetrizable, 
if there exists a smooth A^{uj v) symmetric, positive definite, such that ^^(w, v)A{u, v) 
and A'^{ujV)Q{u,v) are symmetric, and A^{u,v)Q{u,v) is negative semidefinite, 
where, as in the introduction, A := [df, dgY and Q := (0, dqY- 

Lemma 7.4 ([SK]). Assuming simultaneous symmetrizability, condition (H3) is 
equivalent to either of: 

(Kl) There exists a smooth skew- symmetric matrix K{u) such that 

(7.12) Re{KA-A'^Q){u)>e>0. 

A^ as in Definition 7.3. 

(K2) There is no eigenvector of A{u) lying in the kernel of Q{u). 

Proof. These and other useful equivalent formulations are established in [SK] . | 

Next, we recall the basic iteration scheme of [Z.2]; for precursors of this scheme, 
see [Go.2,K.l-2,LZ.l-2,ZH,HZ.l-2]. Restricting now to Lax-type shocks, define the 
nonlinear perturbation 

(7.13) U{x,t):=h){x + S{t),t)-f''_){x), 
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where 5{t) (estimating shock location) is to be determined later; for definiteness, 
fix 5{0) = 0. Then, 

(7.14) Ut-LU = m{U, U), + (0, IrYN2{U, U) + 6{t){U, + U,), 



(7.15) U:=[^l 

where 

(7.16) Nj{U,U) = 0{\Uf) 

so long as |t/| remains bounded. By Duhamel's principle, and the fact that 

/oo 
G{x,t;y)UMdy = e'^'U^ix) = U^{x), 
-oo 

we have 

/oo 
G{x,t]y)UQ{y)dy 
-OO 
ft /'OO 

-/ / Gy{x,t-s-y){N^{U,U)+5U){y,s)dyds 
(7.18) Jo ^-oo 

pt /•oo 

+ / / G{x,t-s-y){QJrYN2{U,U){y,s)dyds 

Jo J-oo 

Defining, by analogy with the linear case, the nonlinear instantaneous projection: 
ip{x,t) := -6{t)U^ 

oo 

E{x,t;y)Uo{y)dy 

oo 
t /-oo 



(7.19) 

pi /'OO 

-/ / Ey{x,t-s;y)iN^{U,U) + SU){y,s)dy 

Jq J-oo 



'0 -'-oo 

or equivalently, the instantaneous shock location: 



oo 



6{t) = - e{y,t)Uo{y)dy 



(7.20) " °" 



t p + CXD 

+ 1 / ey{y,t-s){Ni{U,U)+6U){y,s)dyds, 



-'-oo 
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where E, e are defined as in (1.34), (7.3), and recaUing (7.5), we thus obtain the 

reduced equations: 

(7.21) 

/CXD 
{H + G){x,t;y)Uo{y)dy 
-oo 



-oo 
"t /•oo 



pi /'OO 

+ / / H{x,t-s;y){N^{U,U), + {0,IrYN2{U,U) + 6U,)iy,s)dyds 

Jo J -oo 



-'-oo 
t />oo 

/ Gy{x,t- s- y){Ni{U, U) + 6U){y, s)dyds 

J -oo 

t />oo 

t: 



+ / / G{x,t- s;y){0,lr) N2{U,U)dyds, 

Jo J-oo 

and, difi'erentiating (7.20) with respect to t, 

/oo 
et{y,t)Uo{y)dy 
-oo 

t /• + 00 

i-oo 



(7.22) 



/•I p-too 

+ / eyt{y,t-s){Ni{U,U) + 6U){y,s)dyds. 

Jo J-oo 



Notes: In deriving (7.22), we have used the fact that ey(y, s) ^ as s ^ 0, as 
the difference of approaching heat kernels, in evaluating the boundary term 

/ + 00 
ey{y,0){Ni{U,U) + 5U){y,t)dy = 0. 
-oo 

(Indeed, |ey(-,s)|ii — > 0, see Remark 2.6, below). We have also used the fact that 
E{0,lr) = 0, to eliminate the term corresponding to nonlinear source {0,Ir)N2 in 
(7.19)-(7.20). Moreover, we have in the first place used the fact that the dimension 
i of the stationary manifold is one, a property of Lax and undercompressive shocks, 
in order to factor out U{x) in the first line of (7.19); for this reason, our analysis 
is inherently limited to the case of Lax shocks (recall that the favorable Green's 
function bounds of Theorem 1.10 are available for Lax and overcompressive shocks). 

The defining relation 6{t)uj: := —ip in (7.19) can be motivated heuristically by 

u{x,t)-^{x,t)r^u=(^iyx+6{t),t)-(^''_yx) 

r^U{x,t) + 6{t)U^{x), 

where U denotes the solution of the linearized perturbation equations, and U the 
background profile. Alternatively, it can be thought of as the requirement that the 
instantaneous projection of the shifted (nonlinear) perturbation variable U be zero, 
[HZ. 1-2]. 
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Lemma 7.5 [Z.2]. The kernel e satisfies 

(7.24) \ey{;t)\L.,\et{;t)\L.<Ct-^^'-'/p\ 

(7.25) |et,(-,t)|L. <Ct-^(i-i/p)-i/2, 

for all t > 0. Moreover, for y < we have the pointwise bounds 

(7.26) \ey{y,t)l\et{y,t)\ < Ct-h-^^^^^rr^ , 



(7.27) \ety{y,t)\<Ct-^e m^, 

for M > sufficiently large (i.e. > Ah±), and symmetrically for y > 0. 
Proof. For definiteness, take y < 0. Then, (7.3) gives 

(7.28) ey{y, t) = ( ^ ) {K{y + a_t, t) - K{y - a.t, t)) , 

\u+ -u-j 

(7.29) et{y, t) = ( ^ ] {{K + Ky){y + a.t, t) - {K + Ky){y - a.t, t)) 

\u+ — U- J 

(7.30) 

ety{y,t) = I J {{Ky + Kyy){y + a-t,t) - {Ky + Kyy){y - a-t,t)) , 



where 



(7.31) K{y,t):='' 



V/4:b-t 



4nb-t 



denotes an appropriate heat kernel. The pointwise bounds (7.26)-(7.27) follow 
immediately for t > 1 by properties of the heat kernel, in turn yielding (7.24)- 
(7.25) in this case. The bounds for small time t < 1 follow from estimates 



ry-a-t 
\Ky{y + a-t,t)-Ky{y-a-t,t)\ = \ / Kyy{z,t)dz 

J v+a — t 



ry~°'-^ 2 
(7.32) < Ct-^/2 / e^dz 



' y+a-t 

< Ct~^'^e ^ 
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and, similarly, 

\Kyy{y + a-t,t) - Kyy{y -a-,t)\ 
(7.33) 



f 


a_ 


t 




= 1/ 




^yyy\^^ 


t)dz\ 


J- 


-a- 


-t 




<ct- 


-2 


rv-a-t 
/ 


Mt dz 


<ct- 


-1 


(y+a_t)^ 
e Mt 





The bounds for \ey\ are again immediate. Note that we have taken crucial account 
of cancellation in the small time estimates of ej, ety. □ 

Remark 7.6: For t < 1, a calculation analogous to that of (7.32) yields |ey (y, t) \ < 

(y+a_t)'^ 

Ce Mt ^ and thus |e(-, s)\li ^ as s ^ 0. 

With these preparations, we are ready to carry out our analysis: 

Preliminary result. To illustrate the method, we first carry out the analysis com- 
pletely for the simplest case of discrete kinetic models with the minimal regularity 
assumptions 

(7.34) |t/o|LinHi<Co, 

(still keeping the weak shock assumption (7.11)). That is, we establish the second 
assertion of Theorem 1.13. Define 

C(t):= sup {\U{.,s)\L^{l + s)'/'' + \5{sm + s)-^+\5{s)\ 

(7.35) °-'-* 

We shall establish: 

Claim. For all t > for which a solution exists with ( uniformly bounded by 
some fixed, sufficiently small constant, there holds 

(7.36) m<c2{Co+atf)- 

From this result, it follows by continuous induction that, provided Co < 1/4C2, 
there holds 

(7.37) at) < 2C2C0 

for all t > such that ( remains small. By standard short-time theory/continuation, 
we find that the solution (unique, in this regularity class) in fact remains in L'^CiH^ 
for all t > 0, with bound (7.37), yielding existence and the claimed L^ bounds. 
Thus, it remains only to establish the claim above. 
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Proof of Claim. We must establish bounds on |6''|i2, \S\, \S\, |t/x|L2, and \U\loo. 
Provided that we can estabhsh the others, the L°° bound follows by Sobolev esti- 
mate 

\TT\ -r < \1J\^^'^\1J 1"*"/^ 

Likewise, |t/x|L2 may be estimated in by now routine fashion using a nonstandard 
energy estimate of the type formalized by Kawashima [Kaw]. The origin of this 
approach goes back to [Ka,MN] in the context of gas dynamics; see, e.g., [HoZ.2] 
for further discussion/references. 

We first perform a standard, "Friedrichs-type" estimate for symmetrizable hy- 
perbolic systems. Expanding and differentiating (7.14), and using the key fact that 
Ni{U, U) = in the case of discrete kinetic models, we obtain 

(7.38) U,, - (AU),, - (QU), = (0, IrYN^iU, U), + 6{t){U,, + U,,), 
where 

(7.39) N2{U,U), = 0{\U\^ + \U\\U,\) 

provided |t/|Loo remains uniformly bounded as we have assumed. Let A^ = A^{x) 
denote the simultaneous symmetrizer of A, Q. Taking the L^ inner product of 
A'^{x)Ux against (7.38), we obtain after rearrangement/integration by parts, and 
applications of Young's inequality, the energy estimate 

(7.40) ^{A^Ux,Ux)t<{Ux,A^QUx)+e\Ux\l,+e-^0{\U\l, + \6{t)\^), 

where £ > is as in (7.11). Here, we have freely used assumption (7.11) and the 
consequent higher order bound \U"\ < Ce, plus the bounds \Uj:\lv, \Uxx\lp < C- 

Next, we perform a nonstandard, "Kawashima- type" derivative estimate. Let 
K = K{x) denote the skew-symmetric matrix described in Lemma 7.4, associated 
with A, Q. Taking the L^ inner product of Ux against K times the undifferentiated 
equation 

(7.41) Ut - {AU)x -QU= (0, IrYN2{U, U) + 5{t){lJx + t/,), 
and noting that (integrating by parts, and using skew-symmetry of K) 

]^{Ux,KU)t=\{Ux,KUt) + \{UxuKU) 

(7.42) = \{Ux,KUt) + -]^{UuKUx) - \{UuKxU) 

= {Ux,KUt) + ^{U,KxUt), 
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we obtain after rearrangement the auxiliary energy estimate: 

(7.43) ^(t/,, KU)t < -{U,, KAU^) + e\U^\l, + e-^0{\U\l, + \m\^). 

Adding (7.40) and (7.43), and recalling (7.12), we obtain, finally: 

(7.44) 

^((AOt/,, U^) + (t/,, KU))t < \{U^. (AOQ - KA)U^) + e-^0{\U\l, + \m?) 

<C{\U\l. + \m\'). 

provided e is sufficiently small that e\Uj:\'^j^2 may be absorbed in the negative definite 
termi(t/„(AOQ-KA)t/,). 

Integrating (7.44) from to t, and recalling definition (7.35), thus yields 

(7.45) ((A^t/^t/,) + {U^,KU))\l <C f {\U\l, + \5\^){s)ds. 

Jo 

Assuming, then, that claim (7.36) holds when restricted to \U\l2 and \d\, i.e., 

(7.46) \U\L^is) < CiCo + at?){l + «)-'/' 
and 

(7.47) \s\{s)<C{Co + at)'){l + s)-'/^ 
for all < s < t, we obtain 

{{A^U^, U^) + (t/., KU))\l < C{m + Citf? Al + s)-'/^ ds 

{(Ati) Jo 

<c{m+c{t)Y{i+t)^/\ 

Using Young's inequality to bound 

{U^,KU)it) < e\U,\l,{t) + Ce-^\U\l,{t) 
(7 49 ) 

and recalling, by assumption, that {{A'^UxjUx) + {UxjKU)){0) < C^o, we may 
rearrange (7.48) to obtain 

(A'^t/„t/,)<c(c(0) + c(t)2)(i + t)V2, 

which, by uniform positive definiteness of A'^, verifies (7.36) for IU^Il'^ as well. 
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Thus, in order to establish the result, we have only to verify claim (7.36) when 
restricted to |t/|L2, |5|, and \d\: that is, to establish (7.46), (7.47), and 

(7.50) \S\{s) < C {Co + m'), 

for all < s < t. By (7.21)-(7.22), and A^i = 0, we have 

oo 



(7.51) 



(7.52) 



and 



(7.53) 



\U\^,{t)< I / {H + G){x,t;y)Uo{y)\^,dy 

— oo 

ft POO 



+ \ / H{x,t-s;y){0JrYN2{U,U)dyds\ 

'o J-oo 

"t poo 



/•z poo 

+ \ H{x,t- s;y)dU^{y,s)dyds\^2 

Jo J-oo 

pt poo 

+ \ Gy{x,t-s;y)6U{y,s)dyds\^2 

Jo J-oo 

pt /"CXD 

+ 1/ / G{x,t-s-y){Q,IrYN2{U,U)dyds\ 

Jo J-oo 



--: Ia + h + Ic + Id + le 



oo 



\8\{t) = \ / et{y,t)Uo{y)dy\ 

J-oo 

pt P + CO 

+ \ / eyt{y,t- s)SU{y,s)dyds\ 

Jo J-oo 



Ila + III 



b- 



/oo 
e{y,t)Uo{y)dy\ 
-CXD 

+ 1/ / ey{y,t- s)dU{y,s)dyds\ 

Jo J-oo 



Ilia + nil 



b- 



We estimate each term in turn, following the approach of [Z.2]. 
The linear term /„ satisfies bound 

(7.54) Ia<CCo{l + t)-'/\ 
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as already shown in the proof of Corollary 7.2. Likewise, applying the bounds of 
Lemma 7.1 together with definition (7.35), we have: 

/ H{x,t- s;y){0,IrYN2{U,U)dyds\^, 

J -oo 

t 



(7.55) 



(7.56) 



(7.57) 



and 



(7.58) 



Jo 

Jo 



<catf{i+t) 



-1/4 



ft poo 



pi /•oo 

Ic=\ / H{x,t- s;y)6Uj:{y,s)dyds\^2 

Jo J-oo 

<C f e-^(*-")|(J||t/^U2(s)rfs 
Jo 

<CC{tf I e-''(*-^)(l + s)-i/4rfs 
Jo 



<CQ{tf{l + t)-^'\ 



Id = \ / Gy{x,t- s;y)6U{y,s)dyds\^2 

Jo J -oo 

<C [ {l + {t-s))~^^^\6\\U\L2{s)ds 
Jo 

<CC{tf f (l + (t-s))"'/'(l + s)-3/4rfs 
Jo 



<catni+t) 



2n i/^-V4 



ft poo 



pi fOO 

h = \ / G{x,t-s-y){Q,IrfN2{U,U)dyds\ 

Jo J -oo 

<C [ (l + (t-s))"'/Vli2(s)rfs 
Jo 

< Catf f (1 + (t - .))"'/'(! + s)-^/^ds 
Jo 



<CQ{tf{l+t)-^'\ 
Summing bounds (7.54)-(7.58), we obtain (7.46), as claimed. 
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Similarly, applying the bounds of Lemma 7.5 together with definition (7.35), we 
find that 



(7.59) 



and 



(7.60) 



while 



(7.61) 



and 



(7.62) 



/oo 
et{y,t)Uo{y)dy\ 
-oo 

<\et{y,t)\L^{t)\Uo\L^ 
<CCo(l+t)-V2 



IIi,= \ / eyt{y,t- s)6U{y,s)dyds\ 

Jo J -oo 

< I \eytW{t - s)\5\\U\L^{s)ds 
Jo 

Jo 

<C7C(t)2(l+t)-V2, 

/oo 
e{y,t)Uo{y)dy\ 
-oo 

<\e{y,t)\L^{t)\Uo\L^ 
<CCo 

j-t p + co 

IIh=\ / ey{y,t- s)6U{y,s)dyds\ 

Jo J -oo 

< I \ey\L^{t - s)\d\\U\L2{s)ds 
Jo 

< CC{tf [ {t- s)-^/\l + s)-^/^ds 

Jo 

< CQ{tf. 

Summing (7.59)-(7.60) and (7.61)-(7.62), we obtain (7.47) and (7.50), as claimed. 
This completes the proof of the claim, and the result. | 

Remark 7.7. The hypotheses Uq & L^ r\ H^ made in this subsection on initial 
perturbations are similar to but slightly weaker than those for the stability result 
announced in [S,SX.2] for weak profiles of the discrete kinetic Broadwell model. We 
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obtain sharp rates of decay in L? and smallness in L°°, with nonsharp rates of decay 
in L^, 2 < p < oo, whereas they assert decay in L^, aU p > 2, with no rates. (RecaU 
that spectral stabihty holds for weak Broadwell shocks, by the results of [KM], so 
that we obtain a complete stability result in this case). 

The general case. We next sharpen our analysis for simultaneously symmetrizable 
systems under the strengthened regularity assumption 

(7.63) \Uo\L^nH^ <Co, 

and assuming that q E C^ for discrete kinetic models as well as general systems. 
That is, we establish the first assertion of Theorem 1.13. 
Similarly as in (7.35), define 

C(t):= sup {\U{-,s)\Hi{l + s)'/^ + \5ism + s)'/' + \5{s)\ 

(7.64) °-^-* 

+ \U.A; s)\Ml + s)-^/^ + \U{; s)Uoo(l + s)i/2 + \U^{., s)Uc 

We shall establish, as before, the claim (7.36), from which the desired sharp L^ 
and L°° bounds immediately follow, and thus, by interpolation, the claimed sharp 
bounds in all L^ such that 2 < p < oo. 

Proof of claim. As in the previous proof, we first observe that the bound on 
\Ux\l°° follows by Sobolev embedding, provided that we can establish the claimed 
H^ and H^ bounds on U. 

As before, we next show by energy estimates that the H^ estimate follows pro- 
vided we can establish the claimed bounds on \U\h^, |j|, and \U\l°°. However, we 
must now take a bit more care in this step, due to the presence of the higher-order 
nonlinearity term Ni{U,U)x in (7.14) in the general case. To avoid this term, we 
rewrite (7.14) more strategically, using the quadratic Leibnitz relation 

(7.65) A2U2 - AiUi = A2iU2 - Ui) + {A2 - Ai)Ui, 
in the quasilinear form 

(7.66) Ut - AUx -QU = Mi{U)Ux + (0, IrYM2{U)U + 6{t){Ux + t/,), 
where 

A:={df,dgY{u,v){x + d{t),t), 
Q:={0,dqY{u,v){x + d{t),t), 

and 

Mi(U) = 0(\U\) := A(x, t) - Aix), 
(7.68) S 

M2iU) = Oi\U\) := Qix,t) - Q{x), 
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A{x) and Q{x) defined as usual as {df, dgY{u, v) and (0, dqY{u, v), or, differentiating 
twice: 

(7.69) 

t/xxt - (AU,).. - (QU),, = {M{U)U^)^^ + (0, Ir)\M2{U)U)^^+6{t) {U,.. + U,.,). 

Let A^ denote tlie simultaneous symmetrizer of A, Q. Taking the L^ inner 
product of A^{x)Uxx against (7.69), we obtain after rearrangement/integration by 
parts, and several applications of Young's inequality, the energy estimate 

{A'^Uxx,Uxx)t<{Uxx,A''QUxx)+e\Uxx\l2+e-^0{\U\l, + \6{t)\''), 

where £ > is as in (7.11). Here, we have freely used assumption (7.11) and the 
consequent higher order bounds \U"\ < Ce, \tJ"'\ < Ce, plus the bounds \Ux\lp^ 
\Uxx\lp- \Uxxx\lp < C, and also used the fact that {\U\l°° + \Ux\l°°) < C({t) is 
assumed to remain small, to bound the new term 

(7.70) {AlUxx. Uxx) = 0{\Uxx\LAUt\L^) = 0{\Uxx\l<\U\l^ + \Ux\l^) 

arising from the fact that A^ is now time-dependent. 

Similarly, taking the L^ inner product of Uxx against K times the singly differ- 
entiated equation 

(7.71) Uxt-{AUx)x-{QU)x = {M{U)Ux)x + {0,Ir)HM2{U)U)x+S{t){Ux. + Uxx), 

where K denotes the skew-symmetric matrix described in Lemma 7.4, associated 
with A, Q, and proceeding as before, we obtain the auxiliary energy estimate: 

(7.72) ]^{Uxx, KUx)t < -{Uxx, KAUxx) + epxxW^ + e-^0{\U\l, + \5{t)\'). 

From here, we may proceed as in the previous case to obtain the desired bound on 
\U\h^, assuming the claimed bounds on |t/|//i and \5\. 

But, the bounds on \U\h^, |5|, and \5\ follow exactly as in the previous case, 
once we observe that Gx breaks up into the sum of a time-exponentially decaying 
delta-function term Hi analogous to H and a kernel Gi obeying the same bounds 
as G. This term therefore gives the same bounds as in the L^ estimates of the 
previous case. The term involving H, on the other hand, may be estimated (see 
comment just beneath (7.75), below) as 

1/ / H{x,t-s-y)0{\Uxx\\U\ + \Ux\^){y,s)dyds\^, 

<C [ e-^^'-'^\\{\Uxx\LAU\L^ + \Ux\LAUx\L^){s)ds 
Jo 

<cc{tf I e-'^^'-'Xi + sy^/^ds 

Jo 

<cc(t)2(i + t)-V4, 
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again giving the desired result. 

Finally the claimed bound on |t/|Loo follows, for terms involving the regular kernel 
G, by the exactly the computations used to estimate |^|, and for terms involving 
the singular kernel H, using the new bound |t/^x|L°° < CC,{t) coming from definition 
(7.64), and proceeding similarly as in the estimation of Uh'^. We omit the details, 
which are tedious but straightforward. | 

7.2.2. Discrete kinetic models. Finally, we briefly indicate the extension 
of our arguments to general L^, 1 < p < oo, and strong relaxation profiles, for 
discrete kinetic models with sufficiently regular coefficients q & C^ , not necessar- 
ily simultaneously symmetrizable, and perturbations that are sufficiently small in 
W^'^ n W^'°° . That is, we establish the third, and last assertion of Theorem 1.13, 
completing the proof of the theorem. 

The main idea in this subsection is to take advantage of the special property 
of discrete kinetic models that nonlinear terms, in the original, unshifted equa- 
tions (1.1), appear without derivatives. Likewise, the basic, unshifted perturbation 
variable 

(7.73) U ■.= {u,v){x,t)-{u,v){x) 
satisfies an equation with this same property: 

(7.74) Ut - {AU)^ -QU= (0, /.)*iV2(t/, U), 
where N{U, U) = (9(|f/p), and, more generally, 

(7.75) \{d/dx)''N{U,U)\<C J2 \{d/dx)'^U{d/dxfU\, 

0<a+P<k: a,p>0 

k = 0, . . . , 3, provided |f/|p^i,oo remains bounded, as we shall assume it does 
throughout our argument. (Note: a straightforward Leibnitz calculation shows that, 
provided coefficients q are sufficiently regular, bound (7.75) extends to k = 3K + 2, 
provided |C/|^x,oo remains bounded). 

Using this observation, and the previously observed fact that x-derivative bounds 
on the Green's function, modulo time-exponentially decaying singular terms (i.e., 
delta- functions and their derivatives), are essentially independent of the order of 
differentiation, we may use Green's function bounds instead of energy estimates to 
bound derivatives of the solution. That is, the integral equations for U close, with 
no loss of derivatives, in the special case of discrete kinetic models. On the other 
hand, the variable U is not expected to decay, since it includes a possible stationary 
component corresponding to translation of the original profile. Thus, the bounds 
we obtain by this technique are quite crude, namely 

(7.76) \Uxxx\lp < C\Uj:xx\lp ~ ^, 
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where U as usual denotes the shifted perturbation variable U{x,t) := {u,vY{x + 
d{t),t) — {u,vY{x). However, this is sufficient for our argument, since as we have 
seen in previous subsections, \{d/dx)^U\LP bounds improve as k decreases, by order 
t~^'^ for each successive derivative. At level \Ux\lpi we therefore obtain a uniform 
(9(1) bound, which, as in the proof of the previous case, is sufficient to establish 
sharp bounds in all L^, 1 < p < oo. Note that the resulting order-one |t/||yi,oo 
bound, which implies a corresponding bound on |t/|^i,oo, is consistent with the 
assumption under which we obtained bound (7.75), an important point. 

Remark 7.8. The equations for the more favorable (i.e., decaying) shifted 
perturbation U feature the nonlinear forcing term d{t)Ux, which does involve a 
derivative. Thus, the integral equations for U do not close, involving a loss of one 
derivative, and so we cannot directly estimate jt/xxiUp by this approach. 

We now describe the argument in somewhat further detail. Define 

(7.77) 

C(t):= sup f\U{-,s)\L^{l + sy/' + \U{-,s)\L^ + \d{s)\{l + sf^+\d{s)\ 

0<s<t ^ 

+ |t/'x(-,s)|LinL°° + |^xx(-,s)|LinL°°(l + S)'-^^^ + |t/'a;^x(-,s)|LinL°°(l + ■§)" 

As usual, our goal is to establish 

(7.78) at) < c2{Co + atf) 

for all t > for which a solution exists with ({t) sufficiently small, from which we 
may conclude as before both global existence and the bound 

(7.79) at) < CCo 

for Co sufficiently small. From (7.79) and the definition (7.77), we then immediately 
obtain the claimed sharp bounds in L^ and L^, and, by interpolation, sharp bounds 
for all 1 < p < oo, completing the proof. 

Observation 7.9. Assuming q E C^, (H0)-(H4), and (P), there holds for any 
< /c < 3 the decomposition: 



(7.80) {d/dx)''{G-E) = Gk + J2^^ 

j=o 



k 

J' 



where E is as defined in (1.34); kernel Gk is a measurable function in C^{y), satis- 
fying the same bounds asserted for G in Lemma 7.1; and 

(7.81) Hk{x, t- y) = ^k{x, t; y) {d/dy)^ 
where /Ua;(x, t; ■) are measures satisfying 

(7.82) |//fc(x,t;-| = 0(e-^*). 



CORRADO MASCIA AND KEVIN ZUMBRUN 95 

Proof. The regularity q E C^ yields Q E C^ for the coefficients of the linearized 
operator L, whence we obtain differentiability of order C^{x) (resp. C^{y)) for 
solutions of the eigenvalue (resp. adjoint eigenvalue) ODE from which the resolvent 
kernel Gx is constructed: more precisely, the functions V{x;X) (resp. V{y;X)) 
from which modes e^'^^V{x; X) (resp. e~^^^V{x; X)) are constructed are uniformly 
bounded in C^{x) (resp. C^{y)). 

With this observation, the claimed bounds follow in straightforward fashion, by 
the same procedure used to estimate G, Gn^, Gy in previous sections. Indeed, the 
proof is somewhat simpler, since we only require modulus bounds here, so need not 
carry out the detailed bookkeeping of previous sections. That bounds for the regular 
part Gk do not degrade with increasing k is clear from the elementary calculation 

{d/dx)e^^''V{x; A) = fie^^''V{x; X) + e^^"" {d / dx)V {x- A), 

together with the fact that, modulo negligible error terms, bounds on the regular 
part come entirely from a bounded set of the complex plane, on which ^ is bounded. 
Likewise, the singular terms H^ derive from large-A bounds by explicit computation, 
similarly as in previous sections. | 

We first apply Observation 7.9 in crude fashion to obtain the crucial bound 
\U\w^>^nw^'°° that will be used to close our iteration; this step replaces the energy 
estimates of the previous subsections. Clearly, it is sufficient to bound instead 
|t/|vi/3.inVK3.°°7 U defined as in (7.73) above, using 

(7.83) U{x - 6{t),t) - U{x, t) = U{x) - U{x - 6{t)) ~ 6{t)U{x) 

and assuming the desired order one bound on d{t). 

By (7.74), Observation 7.9, and Duhamel's formula, we have the representation: 

(7.84) 



U^^Ax,t)= / ({E^^^ + Gs){x,t;y)Uo{y) + J2Hk{x,t;y){d/dy)''Uo{y))dy 
-^-°° fe=0 

/•t /•OO 

+ / / {E,^^ + Gs){x,t-s;y){0JrYN2{U,U){y,s)dyds 

Jo J -OO 

3 I't />oo 

+ V/ / Hk{x,t-y){d/dy)\QJrYN2{U.U){y,s)dyds, 



from which we readily obtain the desired bound using the estimates of Observation 
7.9, together with lower derivative bounds 



\{d/dxYU\LvnL^{s) < , 
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< s < t, which follow from the U bounds of definition (7.77) by an argument 
similar to that given in (7.83) above. 
For example, we may bound the term 

t POO 

/ Hsix,t;y)0{\U\\U^^^\){y,s)dyds,\^,^^^ 

J -co 

arising in the expansion via (7.75) of the worst-case singular term 

(7.85) / / Hs{x,t;y)id/dyf{0,IrYN2iU,U){y,s)dyds, 

Jo J -oo 

by 

C [ e-^^'-'^\U\Lo^\U,.,\LinL^{y,s)dyds<Catf [ e-^^'-'\l + s)ds 
Jo Jo 

and similarly for all other terms arising in (7.85), since by homogeneity of our 
derivative estimates the bound on the nonlinear term depends only on the total 
number of derivatives. 

Likewise, the worst-case regular term 

(7.86) 1/ / E,^^{x,t-s-y){QJrfN2{U,U){y,s)dyds\^,^^^ 

Jo J-oo 

may be bounded using the triangle inequality by 

C [ \E,.^\L^nL^{t - s)\N\Li{s)ds < CCit)^ [ ds 
Jo Jo 

<cat)\i+t), 

completing the verification of the claimed bound for \Uxxx\LpnL°°- 

From this highest-derivative bound, the lower-derivative bounds on the shifted 
perturbation U may be obtained in routine fashion from the corresponding integral 
representations for U, using the bounds of Observation 7.9 and estimates precisely 
analogous to those of the previous subsections. We omit these lengthy but straight- 
forward calculations. This completes the proof of the claim, and the Theorem | 

Remark 7.10. In this subsection, we have made little attempt to minimize 
regularity requirements. It seems likely that the assumption W^'^ n VF^'°° on initial 
perturbation could be substantially reduced. 
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A. Appendices: Auxiliary Results/Calculations. 
Appendix Al. Structure and existence of Relaxation Profiles. 

In this appendix, we prove the results cited in the introduction concerning struc- 
ture and (small-amplitude) existence of relaxation profiles. 

Proof of Lemma 1.2. First, note that matrix (1.13) can have no nonzero purely 
imaginary eigenvalues i^, by (H3). Observing that 

det (Q„,Q„)f^" ^A fM=detf^" fv\ffu /.' 

\9u 9v J \Ir J \qu qv J \9u 9v 

Ju Jv 



det df* det o^det 

' 9u 9v 

is nonzero, by (HI), (H2), and (1.2), we find that g+ is a hyperbolic rest point of 
(1.11), from which (1.14) follows. | 

Proof of Lemma 1.3. In contrast to the analogous result of [MP] in the 

viscous case, which can be phrased a,s i± = d±, n = r, independently for each state 
ti±, the present result must be stated jointly, and in fact depends in an essential 
way on the fact that the two states are joined by an admissible connecting profile. 
Moreover, the proof emerges in a natural way from the investigation of stability 
properties through the study of the linearized eigenvalue ODE. A similar argument 
was sketched in Remark 2.3, [ZH], for the parabolic case. 

Precisely, recall from Lemma 3.1, Section 3, that the (n -|- r)-dimensional eigen- 
value ODE satisfies the consistent splitting hypothesis at spatial plus/minus infinity, 
as a consequence of (H3): that is, the dimensions of the stable/unstable subspaces 
of the "frozen," limiting coefficient matrices at plus and minus infinity have common 
values, for all frequencies Re A > 0. On the other hand, by Lemma 5.1, Section 
5, the eigenvalues of the frozen coefficient matrices at A = consist of r "fast" 
modes Re fXj ^ 0, corresponding to the eigenvalues of the reduced traveling wave 
ODE linearized at ('U±,i'±), and n "slow" modes \ij = corresponding to the n 
constants of motion of the full traveling wave ODE. At plus infinity, the former 
contribute (i_(_ decaying values, while the latter bifurcate as Re A increases into z_(_ 
growing modes and (n — i_|_) decaying modes, accounting for a total of (i_(_ -|- (n — i-i-) 
negative real part eigenvalues of the coefficient matrix for small A with Re A > 0. 
Likewise, at minus infinity, we find that there are (i_ -t- (n — z_) positive real part 
eigenvalues for small A with Re A > 0. Invoking consistent splitting, we have that 
(i_i_ + (n — Z-I-) + (i_ + (n — z_) = n + r, yielding the result. | 

Remark. In the strictly parabolic case, the argument above yields instead that 

d±-^ {n — i±) = n, or d± = i±, the result of Majda and Pego [MP]. 

Existence of weak profiles. For completeness, we present also a general result 
on existence of small-amplitude profiles, encompassing and slightly extending that 
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of [YoZ]; in particular, we do not require the "genuine coupling condition" (2.5) of 
[YoZ] , which was induced by the particular choice of coordinates made there. As in 
[YoZ], our analysis is based on a center manifold construction analogous to the one 
of [MP] in the viscous, strictly parabolic case; however, our argument is somewhat 
simpler than that of [YoZ], and reveals more clearly the analogy to [MP]. 

Proposition 8.1. Let (HO)-(Hl), 

(H2) ex {df*{u±)) real and distinct, 

and (H3) hold for {u,v) in a neighborhood of an equilibrium state {u±,v±)o ■= 
{uo,v*{uo)), let So :— ap^uo), where a*{u) denotes the pth eigenvalue a*(w) (char- 
acteristic speed) of the equilibrium flux Jacobian A* := df*{u), and let U be a 
sufficiently small neighborhood of {uq,v*{uq)). Then, for each equilibrium states 
{u±,v*{u±), s) in a sufficiently small neighborhood of {uq^vq, sq), such that the 
Rankine-Hugoniot condition f*{u+) — /*(«_) = s{u+ — U-) and the noncharac- 
teristic condition a*(w-|-) 7^ s for the equilibrium system (1-4) are satisfied, i.e., 
(f^i, s) forms a noncharacteristic shock triple for (I.4), there exists a traveling wave 
connection (1.7) lying in V if and only if the triple {u±^ s) satisfies the Liu-Oleinik 
entropy condition [L.3] for the equilibrium system (1.4)- 

Such a local connecting orbit, if it exists, is unique. Moreover, for triples that 
are uniformly noncharacteristic in the sense that |a*(w±) — s\/\u+ —U-\ > 9, 9 > 
a fixed constant (note: this holds for all triples when characteristic a* is genuinely 
nonlinear at uq), the local profiles satisfy uniform bounds 

(8.1) {u-u±,v-v±){z) < C\u+ - U-\e-K(^^)-'\\^\/(^ 

for z > 0, z < 0, respectively. 

Proof. For simplicity of exposition, we shall establish the result with («_, v*{u^) 
held fixed, without loss of generality taking («_, t'*(w), sq) = (0,0,0), and taking 
z < only in (8.1); The full result then follows, with z still held nonpositive in 
(8.1), by compactness, and the observation that all of our estimates are uniform in 
the various parameters of the problem; finally, we may extend (8.1) to 2; > by 
symmetry with respect to U- and u.^.. Alternatively, one may as in [MP] work with 
the "translated variable" (u^v) — (u-,V-), with (u-jV-) carried as an additional 
{n + r) -dimensional parameter, to obtain the full result through a single, larger 
center manifold construction. 

Including now the new parameter s as a supplementary variable, we may write 
the traveling wave ODE as 

{f{u,v)-suy = 0, 

(8.2) {g{u, v) - sv)' = q{u, v), 

s' = 0. 
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The map from (w, v) to (/, g) := (/ — su, g — sv) (with s held fixed) has Jacobian 

Ju S Jv 



(8.3) I JU ■^ JV 

\ 9u 9v - sv 

at (uo^vq) = (0,0) So = 0, hence is locaUy invertible by the noncharacteristic 
assumption (HI). 

Thus, we may rewrite (8.2) in the tilde variables as 

f = 0, 

(8.4) ~g' = q{u,v){fr9), 

s' = 0, 

or, eliminating / = 0, in the convenient form 

(8.5) fg^fQiu^vM-g) 

Linearizing (8.5) about base point (gojSo) := ((7(0,0), 0), and using (8.3), we 
obtain after a brief calculation the linearized equations 

<-) (f)'-(:^)(f 

where 

Evidently, the center subspace of the coefficient matrix of (8.6) is the direct sum of 
(Or, 1) and (C, 0), where C is the center subspace of M. 

Claim. The center subspace C of M is one-dimensional, consisting of its kernel. 
More precisely, it is spanned by (7^r*(0), where r* is the right eigenvector of the 
equilibrium fiux A* associated with the principal characteristic speed a*, and (7* 
analogously to /* is defined as gu — (lv^Qu9v, the variation of g along the equilibrium 
manifold: equivalently, 

(8 8) (^- ^"VV ^ ] = ( f^ 

Likewise, the left kernel of M is spanned by — /p/wQ^^(0), where /* is the left zero 
eigenvector of /* dual to r*. 
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Proof of Claim. Clearly, the eigenvalues of M are contained among the eigenval- 
ues of 

(8.9) M,:=f° OU/. /. 

which has no nonzero pure imaginary eigenvalues, by (H3). Thus, C consists entirely 
of the zero eigenspace of M. 

Next, supposing that x lies in the kernel of M, we find that 

fu fvV' fO 



9u 9v J \yX 

lies in the kernel of {qu,Qv), hence must be of form (r, —Qy^QufY- From the first 
coordinate of the equation 

fu fv\ f r \ _ fO 



9u 9v J \-Qv^Qur J \x 



we find that r must lie in the one-dimensional kernel of /.j^(O), whence we obtain 
(8.8) and the claimed description of the kernel of M. 
A direct computation then shows that 

-i;fvq-'M = -i;{fyq-^q^ /„ ) f f " ^A f^A 

\ 9u 9v / \ 1^ J 

-^lu: o)ff- [A f°- 

\ 9u 9v J \ ^r 

i.e., X := —IpfvQv^ satisfies xM = 0, and more generally, 
(8.10) 

Ju Jv 



-i;fvq-'{qu q.)({; {;) =-i;{f.q-\u f.){- 

^ \9u 9v J ^ \[ 

= -i;{{-f: o) + (^ fy 



Ju Jv 

9u 9v 



-11 {In 0, 



from which we may conclude that a; 7^ 0. This verifies the asserted description of 
the left kernel of M. 
Observing that 

~ 7* p — 1 * * 7* r?* * 

XX ~i^pjvqv 9u'^p *p-" '^p-i 

where B* is the Chapman-Enskog viscosity defined in (1.6) (Note: in the final 
equality, we have used the fact that f^fp = 0), we thus find that a necessary and 
sufficient condition for existence of a generalized eigenvector y. My = x, is that 

(8.11) XX = CB*rl = 0. 
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But, the low-frequency Fourier expansion in the foUowing appendix (A2) shows 
that the stabihty condition (H3) imphes l*B*r* > 0, and so no such generahzed 
eigenvector occurs. | 

Remark. The conclusions of the Claim may be obtained in a lengthier but more 
transparent fashion by the observation that zero eigenvalues/generalized eigenvalues 
of M correspond with those of 

Ju Jv \ I Ju Jv 



Qu Qv I \ Qu 9v 



With these preparations, the result now follows similarly as in [MP] . By standard 
center manifold construction, there exists a two-dimensional manifold tangent to 
the center subspace of M, which contains all locally bounded orbits of (8.5): in 
particular, all rest points (^±, s). For each fixed s, this reduces to a one-dimensional 
fiber 

(8.12) ~9 = m, 

9 E^, containing both of the rest points g±. Thus, there is a connection between 
them if and only if there is no rest point of (8.5) lying between them. 

By the Lax structure theorem [La] , locally, all rest points lie on the pth Hugoniot 
curve through -U- = of the equilibrium system, mapped through g{u,v*{u)) — 
sv*{u) to the g coordinates, and both Hugoniot curve and fiber are smooth curves 
lying approximately in the fi'^r*(0) direction. Thus, the ordering of the rest points 
must be the same along the center manifold fiber as along the Hugoniot curve. But, 
the analysis of [L.3] shows that two noncharacteristic rest points are consecutive 
along the Hugoniot curve if and only if the Liu-Oleinik admissibility condition or 
its reverse is satisfied. Thus, it remains only to check that the sense of a connection 
must agree with that of the (forward) Liu-Oleinik condition, or, equivalently, that 
there hold the one-sided Lax characteristic condition a*{u-) > s; for details, see 
[L.3]. 

To see that the one-sided Lax condition holds, we may use invariance of (8.12) 
under the flow of (8.5) to obtain 

(8.13) h{Orge = q{u,v){0,~g), 

where 

e' = h{e) 

describes the reduced flow along the center manifold. Differentiating (8.13) with 
respect to 6* at 6* = 0, we obtain 

he{0)ge{0) = M{s)g0{O), 
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where 

fu - SU f, 



(8.14) V 9u ^. - « V u,_o,o V ^' 

is a one-parameter perturbation of M above, with 
U5) M2:={qu q.) { {" {^ ^ '^ °" 

That is, /i6>(0), ^6i(0) are a right eigenvalue, eigenvector pair for M(s), bifurcating 
from the values 0, (7^r*(0) for M at s = 0. Since we have already ascertained 
that is an isolated eigenvalue of M, with left and right eigenvectors l*fyq~^{0) 
and fi'u^p(O), respectively, we may conclude by standard matrix perturbation theory 
[Kat] that hg{0) varies with respect to s as 

,.. IJJvq^Hm^g^r^ 

hgis) = s— -, h Ois ) 

(8.16) =s^i^^i^ + 0{s 



-h- I / 1 _>< 

/;s*r;(o) 



^ -'2^ 



where in the second equality we have used (8.10) and (8.8) in simplifying the nu- 
merator of the first-order coefficient. With (8.11), this gives 

sgn hg{0) = sgn s = sgn (s — a*(tt_)) 

for s sufficiently small, verifying that ^ = is a repelling point of the reduced flow 
if and only if the one-sided Lax condition holds. At the same time, it gives the 
bound (8.1) for -z < by straightforward estimates on the scalar system (7.51). 
This completes the proof of the Proposition. | 



Appendix A2. Expansion of the Fourier Symbol. 

In this appendix, we carry out the Taylor expansions about ^ = and ^ = oo of 
the Fourier symbol —i^A{x) + Q{x) of the frozen, constant-coefficient operator at 
a given x. 

High frequency expansion. By straightforward matrix perturbation theory 
[Kat], the ffist-order expansion at infinity of P{iC) with respect to (iO~^ i^ J^^st 

P = -iCA + Rdiag{r]i,...,r]j}L, 
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diag{?7i,...,?7j} := LQR, 

where L, R are composed of left, right eigenvectors of A, as described in the intro- 
duction. 

Low frequency expansion. We next carry out the expansion of P in ^ about 
zero. RecaU that 



L17) 



Q:- 







A 



Ju Jv 

9u 9v 



Claim. To second order, "slow" dispersion relations 

(8.18) X{tO = a{Q-iCA), A(0) = 
are given by 

(8.19) X,{0 = -tCa*-l3*e + ---, j = l,---n, 
with corresponding eigenvectors 

R*{iO = i y . \ + 

and 






^.20) 



/3* := l*B*r*. 



(8.21) 



B* := fv Qu - gvQv ^Qu + fuQv ^Q^ 



where a*, r*, /* are eigenvalues and associated right and left eigenvectors of A* := 
/*. The remaining r "fast", or strictly stable roots are ~ ^{qy). 

Remark. This expansion may be viewed as a rigorous justification of the 
Chapman-Enskog expansion in the constant-coefficient case. 

Proof. Setting 

(Q - iCA - i-i^a, - ^,e) ■■■)iy^+Vii+V2e+---)=Q. 

and collecting terms of successive orders in ^, we obtain: 
(Order 1): 



.22) 



QVo = Q 



Vn 



r 



3 
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(Order ^): 

(8.23) QVi^i{A~aj)Vo = 0, 
=^ (looking at first coordinate) 

(8.24) {A, - a,)r, = 0, A* := fl = /„ - Uq-\u 

=^ aj, Tj are eigenvalue a*, right eigenvector r* of A* . 

Further, comparing with the second coordinate in (8.23), we have 

(8.25) -{qu,qv)Vi = i{gu-:gv-o*)Vo, 

T 



s 



9u - {Qv - a*)q^ ^Qu 



rj) 



giving 

(8.26) Vi = 

(Order ^2): 

QV2 - i{A - aj)Vi + (3jVo = 0. 

Examining the first coordinate, we have 

PjTj =i{fu-a*,fv)Vi 

= i{A* - a*)s - Tfyq~^[gu - {Qv - a*)q~^qu]r*. 

Taking inner product with /*, left e-vector of A*, we obtain 



(8.28) 

= -rl*fvqy [gu-9vqy QufuV^j, 



completing the result. | 

Remark. The "strict asymptotic parabolicity" condition 

(8.29) Re/?*>0, j = l,...,n, 

reduces in the case n = r = 1 treated by Liu [L.2] to the classical subcharacteristic 
condition [Wh]; however, these are not equivalent. As discussed in [Yo.4], the 
relation between the subcharacteristic condition and stability is quite analogous 
to that between the CFL condition and stability in numerical analysis. In the 
symmetrizable case, condition (8.29) is equivalent (see [Kaw,Ze.2]) to Q < and 
no vector in ker{Q) is an eigenvector of A. 



CORRADO MASCIA AND KEVIN ZUMBRUN 105 

Appendix A3. The Gap and Tracking Lemmas. 

For completeness, we present here two technical lemmas referred to throughout 
the text, relating behavior of variable- and constant-coefficient ODE in, respectively, 
the asymptotically constant and the slowly varying coefficient case. The resulting 
bounds are crucial in our study of solutions of the generalized eigenvalue equations, 
in the former case on bounded domains in frequency space A, and in the latter as 
|A| -^ oo. 

A3.1. The Gap Lemma. 

A common problem arising in the asymptotic study of eigenvalue ODE is to relate 
behavior near x = ±cxo of solutions of an asymptotically constant-coefficient eigen- 
value equation to that of solutions of the corresponding limiting, constant-coefficient 
equations, in a manner that is smooth with respect to spectral parameters. More 
generally, consider a general ODE with parameter 

(8.30) W' = A{X,x)W, 

"/" denoting (d/dx), where the coefficient A is C^ in the evolution variable x and 
analytic (resp. C^) in parameter A, and converges as x ^ ±00 to limiting values 
A±. It is well known (see [Co], Thm. 4, p. 94) that, provided that 

/•±CSO 

(8.31) / |A- A±|rfx < +00, 



is a generalized eigenvector of order i) and certain solutions W- of (8.30) having 



Jo 

there is a one-to-one correspondence between the normal modes V- e^i ^ of the 
asymptotic systems 

(8.32) W'^K±{\)W, 

where V- , // are eigenvector and eigenvalue of A_ (alternatively, V^ x^e^i ^, if V^ 
is a generalized eigenvector of ord 
the same asymptotic behavior, i.e. 

(8.33) Wf{\x) = Vj^e^'f^l + o(l)) as a; ^ ±00 

(alternatively, W^{X,x) = V^x^e^^ "^(1 + o(l)). That is, the flows near ±00 of 
(8.30) and (8.32) are homeomorphic. 

Such a correspondence is of course highly nonunique, since (8.33) determines 
W- only up to faster decaying modes. However, provided that Re (u • ) is strictly 
separated from all other Re {p,^ ), i.e., that there is a spectral gap, the choice defined 
in [Co], Theorem 4 by fixed point iteration is in fact analytic (resp. C^) in A, as 
the uniform limit of an analytic (resp. C^) sequence of iterates. The argument 
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breaks down at points Aq where Re (fXj) = Re (nk) for some k 7^ j, since in this 
case (Re (//j) — Re (nk)) does not have a definite sign, and the definition of the 
fixed point iteration is determined by the signs of aU (Re (nj) — Re (/Ufc))- 

The purpose of the present section is to point out that analyticity (resp. smooth- 
ness) in A can be recovered in the absence of a spectral gap, by virtuaUy the same 
argument as in [Co] if we substitute for (8.31) the stronger hypothesis: 

(8.34) |A-A±| =C(e-°l^l) asx^ioo. 

This observation is a special case of the "Gap Lemma of [GZ] , also proved inde- 
pendently in [KS]. The original version was phrased in terms of the projectivized 
fiow associated with (8.30). Here, we give an alternative statement and derivation 
directly in terms of (8.30), a form more convenient for our needs. 

Proposition 8.2. In (8.30), let A be C^ in x and analytic (resp. C^ ) in X, with 
|A — A_(A)| = C(e~"l^l) as X ^ —00 for a > 0, and a < a. If V^~(A) is an 
eigenvector of A- with eigenvalue /u(A), both analytic (resp. C^ ) in X, then there 
exists a solution W{X,x) of (8.30) of form 

Vr(A,x) = y(x;A)e'^(^)", 

where V (hence W ) is C^ in x and locally analytic (resp. C^ ) in X and for each 
J = 0, 1, . . . satisfies 

(8.35) {d/dXyV{x;X) = {d/dXyV-{X) + 0{e-''\''^\V-{X)\), x < 0, 

Moreover, ifKe fx{X) > Re /u(A) — a for all (other) eigenvalues Ji of A-, then W is 
uniquely determined by (8.35), and (8.35) holds for a = a. 

Proof. Setting W{x) = e^^V{x), we can rewrite W = AW as 

v = (A_ -fii)v + ev, 

(8.36) 

^:=(A-A_) = 0(e-"l^l), 

and seek a solution V{x] A) -^ V~{x) as x ^ 00. 
Set 

(8.37) a — di < ai < a2 < a2 < en- 
Fixing a base point Aq, we can define on some neighborhood of Aq to the comple- 
mentary A_-invariant projections -P(A) and Q{X) where P projects onto the direct 
sum of all eigenspaces of A_ with eigenvalues Ji satisfying 

(8.38) Re {]!) < Re (//) + ^2, 
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and Q projects onto the direct sum of the remaining eigenspaces, with eigenvalues 
satisfying 

(8.39) Re (Jl) > Re (jj,) + 0^2 > Re (/u) + en. 

By basic matrix perturbation theory (eg. [Kat]) it follows that P and Q are analytic 
(resp. C"~) in a neighborhood of Aq, with 



^.40) 



g(A_-Ati')xp 
g(A_-M/)xg 



C(e"2^), X > 0, 
C(e"i^), x<0. 



Thus, for M > sufficiently large, the map T defined by 

(8.41) 

/x I'-M 

ei^--f'i')(^-y') pe{y)V{y)dy - / e^^--''^'^^''-y'^Qe{y)V{y)dy 
-00 J X 

is a contraction on L°°(— 00, — M]; for, applying (8.40), we have 
(8.42) 

J — 00 J X 

= 0{l)\Vi - y2|oc( e"^"e("-"^)^|^^ + e^i-e("-"^)^|-^ 
= (9(1)1^1 -y2|oc(e"" + e«") 
= (9(1)1^1 -y2|ooe-^'^<^. 

By iteration, we thus obtain a solution V G L°°(— 00, — M] of y = TV with 
F = (9(|y~|); since T clearly preserves analyticity (resp. smoothness), V{XjX) is 
analytic (resp. C^) in A as the uniform limit of analytic (resp. smooth) iterates 
(starting with Vq = 0). Differentiation shows that y is a bounded solution of 
V = TV iff it is a bounded solution of (8.36). Further, taking Vi = V , V2 = Q in 
(8.42), we obtain from the second to last equality that 

(8.43) \V -V-\ = \T{V) - T(0)| = C(l)e"^|y| = C(e"^)|y-|, 

giving (8.35) for j = 0. Derivative bounds, j > 0, follow by standard interior 
estimates, or, alternatively, by differentiating (8.41) with respect to A and repeating 
the same argument. Analyticity (resp. smoothness), and the bounds (8.35), extend 
to a; < by standard analytic dependence for the initial value problem at x = —M. 
Finally, if Re (//(A)) > Re (//(A)) — ^ for all other eigenvalues, then P = I, 
(5 = 0, and V = TV must hold for any V satisfying (8.35), by Duhamel's principle. 
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Further, the only term appearing in (8.42) is the first integral, giving bound (8.43) 
for a = a. I 

Proposition 8.2 extends also to subspaces of solutions. This can be seen most 
easily by associating to a /c-plane of solutions. Span {Wi{x), . . . , Wk{x)}, the cor- 
responding /c-form r] = Wi A ■ ■ ■ A Wk- The equations (8.30) induce a linear fiow 

(8.44) ?7' = A('=)(x;A)?7, 
on the space of /c-forms via the Leibnitz rule, 

(8.45) A('=)(Wi A---AWk) = {AWi A ■ ■ ■ AWk) + ■ ■ ■ + {Wi A ■ ■ ■ A AWk). 

The evolution of the /c-plane of solutions of (8.30) is clearly determined by that of 
r]{x; A). It is easily seen that for a given (constant) matrix A, the eigenvectors of 
A*^*^) are of form Fi A ■ ■ ■ A Vfc, where Span {l/i, ■ ■ • , Vfc} is an invariant subspace of 
A, and that the corresponding eigenvalue is the trace of A on that subspace. 

Definition 8.3. Let C = Span {Vk+i, . . . , Vn} and E = Span {Vi, . . . , Vk} be 
complementary A-invariant subspaces. We define their spectral gap to be the dif- 
ference (3 between the real part of the eigenvalue of minimal real part of A restricted 
to C and the real part of the eigenvalue of maximal real part of A restricted to E. 

If ?7 is a /c-form associated with an A-invariant subspace E as in the definition 
above, then the spectral gap (3 is the minimum difference between the real part 
of the eigenvalue (i of A*^'^-' associated with rj and the real part of the eigenvalue 
associated with any other eigenvector of A*^'^-'. Combining these observations with 
the result of Proposition 8.2, we obtain a complete version of the Gap Lemma of 
[GZ]: 

Corollary 8.4 (The Gap Lemma). Let A(x; A) be C^ in x, analytic (resp. C^ ) 
in X, with A{x; A) -^ A±{X) as x ^ ±cxd at exponential rate e""'^', a > 0, and let 
rj~{X) and (~ be analytic (resp. C^ ) k andn — k-forms associated to complementary 
A_{X) -invariant subspaces C~ and E~ as in Definition 3.2, with arbitrary spectral 
gap 13, and let tc-{X) be analytic, where tc-{X) is the trace o/A_(A) restricted to 
C~ . Then, there exists a solution W(A, x) of (8.44) of form 

W{X,x) = ?7(A,x)e^c-^ 

where r] (hence W) is C^ in x locally analytic (resp. C^ ) in X, and for each j — 
0, 1, . . . satisfies 

(8.46) {d/dXyr]{x;X) = {d/dXyr]-{X) + 0{e-'^\''\\r]-{X)\), x<0, 

for all a < a. Moreover, if 13 > —a, then rj is uniquely determined by (8.46) and 
(8.46) holds for a = a. 

Remark 8.5. Note that the construction of Corollary 8.4 determines at the same 
time analytic subspaces C{x; A), E{x; A) asymptotic to C~, E~ , and thereby corre- 
sponding complementary analytic projections Pc{x; A), Pe{x; A). Fixing xq and Aq, 
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and applying projections Pc{xo, A), Pe{xq, A) to fixed bases of C(xo, Aq), E{xq, Aq), 
we may thus obtain a locally analytic choice of individual solutions W~ [x; A) span- 
ning C(x; A), E{x;X), in a neighborhood of any Aq. Moreover, if the eigenvalues 
associated with C~ have real parts lying strictly between ryi and 772, then the stan- 
dard homeomorphism between solutions of limiting, constant coefficient equations 
and their variable-coefficient counterparts (described, e.g. in Section 3) yields that 
the flow on C(A, ■) decays/grows exponentially with |a:| on x < 0, with rate uni- 
formly bounded between r]i and 772- Alternatively, such an analytic choice may 
be obtained in more straightforward fashion by a standard "stable-manifold type" 
construction, as described in [Z.l], from which exponential decay/growth may be 
deduced at the same time. 

Remark 8.6. Bounds (8.35), (8.46) assert no information about the direction of 
V{-; A), ??(■; A) at x = 0. A review of our flxed-point construction reveals that in gen- 
eral no such information is available, i.e., T is a contraction only on L°°(— 00, — M], 
for M > sufficiently large, and not on L°°(— 00, 0]. Specifically, the second integral 
in definition (8.41) becomes for x > — M an expansive term with expansion coeffi- 
cient growing exponentially in a; + M with rate at least ai and possibly arbitrarily 
larger; at x = 0, this becomes at least order e"^^, and typically greater than one. 

In particular, we obtain from the bounds of the Gap Lemma no lower bound on 
the Evans function |iI>(A)|, i.e., no information on existence or nonexistence of point 
spectrum at A. Indeed, this could not be so, since we know a priori in many cases 
that -D(O) = even though the Gap Lemma estimates are uniformly valid up to A = 
0. To obtain lower bounds on the Evans function, for example in the high-frequency 
limit, we will use instead an alternative construction described in the following 
subsection, based on the projectivized flow of the eigenvalue ODE. By factoring out 
variations in modulus, this allows more accurate estimates of direction; in particular, 
the mapping analogous to T in projectivized coordinates ((8.54) below) will be seen 
to be a global L°°-contraction on the whole real line.^ 

A3. 2. The Tracking Lemma. 

Another general situation that arises in the asymptotic study of eigenvalue equa- 
tions is an ODE 



^.47) w' = {A{x,5) + 5e{x,5))w, weC 



N 



^In this regard, we point out a key error in the proof of Lemma 5 of [R], specifically of lower 
bound (20) asserted on the Evans function (bounds (18)-(19) clearly extend from a;o to a; = 
by Abel's formula). Namely, the claimed bounds on operator A^ of (34), following directly from 
contractivity of the mapping T\ (our T, above) are valid only on the interval [a;o,oo) where T\ 
is a contraction, and not all of [0, +00) as claimed. The error results from an incorrect citation 
of the Gap Lemma construction of [GZ,ZH], stating that mapping T\ itself is a contraction on 
[0,00)). As discussed in Remark 8.6 above, such information cannot be obtained by the Gap 
Lemma construction, but rather requires the projectivized version described below (indeed, the 
required lower bound in [R] is a trivial consequence of the machinery developed here, in Section 4 
and Appendix A3. 2 below). 
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with a small parameter 5^0, satisfying 

(8.48) |A| + |e| < C 
and 

(8.49) |A'| < Cd, 

where '"" denotes differentiation with respect to x: that is, an ODE with slowly 
varying coefficients. This situation arises in the limit as the frequency rather than 
the spatial variable goes to infinity, after rescaling to a length scale on which the 
resulting rapid oscillations in the solution have period of order one with 5 — i> as 
|frequency| — > oo. Thus, it is dual to the previous case. 

A3. 2.1. The basic estimate. To introduce ideas, we first present an analysis 
of the "standard" case treated in [GZ,ZH,Z.1,Z.4]. Suppose in addition to (8.48)- 

(8.49) that the spectrum of the matrix A(a:, 6) divides into two spectral groups, 

(8.50) ai{x , 6) , . . . ai{x , 5) < a{x) < a{x) < q;^_(_i(x, 5), . . .a]^{x, 5), 
where < denotes ordering with respect to real parts, with a uniform spectral gap 

(8.51) Q < 2r] < a{x) - a{x); 

in the setting of asymptotic eigenvalue equations, this is achieved by restricting to 
an appropriately small subset of the region of consistent splitting, or normal set, of 
the operator L under investigation. By standard matrix perturbation theory [Kat] , 
therefore, we have A(-) -invariant projections P{x) and Q{x) onto the eigenspaces 
associated with (cti, . . . , a^) and (a^+i, . . . a„), respectively, satisfying 

(8.52) \P'\ + \Q'\<C25, |P|,|Q|<C; 

we assume that these bounds hold uniformly over all — cxo < a; < -|-cxd (this can be 
guaranteed, for example, by the assumption that A varies within a compact set, as 
is the case in the applications we have in mind). 

Under these general assumptions, we will show that the "stable" / "unstable" 
manifolds of solutions of (8.47), decaying at +oo/— oo with rates ~ e-^/e"^, respec- 
tively, approximately track the corresponding subspaces of the principal coefficient 
matrix A(a;, 5) as they vary with x, lying always within angle 0{d/r]) and decreas- 
ing/increasing with uniform rate ~ ^9.^ j^olx £qj, ^^ ^ ^ max^: a, a < min^; a. 

Remark 8.7. The most common applications of this result concern the situation 
that 

a(x) = a < < a = (x{x^, 

in which case the manifolds described are truly stable/unstable and a = niiux a, 
a = maxx a. However, neither our results nor their useful application are limited 
to this case. 

We will follow an approach based on energy estimates/invariant cones, though 
other approaches are certainly possible; see [AGJ,GZ], for example, for a revealing 
alternative formulation of Proposition 8.22 (just below) in terms of projectivized 
differential forms. For further discussion, we refer the reader to [ZH], Section 7. 
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Proposition 8.8. For C sufficiently large, 5/r] sufficiently small, the cone 

K_(x):={«;:^^<CVr7} 

^ \Q[x)w\ ^ 

is positively invariant under the flow of (8. 47) (i.e., invariant in forward time 
X /"), and exponentially attracting on JI_(a:) := {w : - q/^) - < rj/C5}. Here, 



Q{x)w 

C depends only on the hounds (8.48)-(8.49) and (8.52) and r] as defined in (8.51) 
measures spectral gap. 

The motivation behind Proposition 8.8 is clear: the P-component of w is in some 
sense growing exponentially slower than the Q-component, by the amount of the 
spectral gap. To quantify this observation, we use the following standard linear 
algebraic result: 

Lemma 8.9. For matrices M s.t. (M) < C and Re 5{M) > rj > 0, and for any 
< rj < r], there exists a C°° choice of coordinate transformation S{M) such that 
M := SMS~^ is real positive definite, satisfying 

Re (M) := -{M + M*) >7j>0. 

Proof. See, e.g.. Proposition A. 9, p. 361 of [St].| 

It is worthwhile to consider the import of the Lemma in the context of a constant 
coefficient ODE 

(8.53) w' = Mw, 

namely, the existence of a coordinate change w = Sw for which the flow of (8.53), 
now of form 

(8.54) w' = Mw, 
strictly increases \w\. For, 

{-\w\'^y = -{{w,Mw) + {Mw,w)) = {w,Re Mw) 

(8.55) 2 2 

- 2' ' 
Related topics are Lyapunov theory and the Kreiss matrix theorem. 

Proof of Proposition 8.8. By rescaling if necessary hy y = e~^''^^—'^^'^W^ we 
can reduce without loss of generality to the symmetric case 

(8.56) tti, ...,«£< -?7 < < ?7 < a|?+i, . . . , aN- 



112 STABILITY OF RELAXATION SHOCKS 

By (8.49), there exists a change of coordinates w -^ S{x)w reducing (8.47) further 
to the case that A{x, 6) has block diagonal form, 

(8.57) A(x-,«)=(^^' ^0^ 

For, by standard matrix theory [Kat], (8.52) implies the existence of bases span- 
ning the left and right eigenspaces associated with projections P, Q, and depending 
smoothly on P, Q; it follows that these bases, and thus the corresponding diagonal- 
izing transformation S have derivatives of order 6, whence error term S~^S'w may 
be absorbed in G. Finally, using Lemma 8.9, we can reduce by a further coordinate 
change w -^ Tw to the case that 

(8.58) Re Ml < -?j, Re M2 > rj. 
Expressing (8.47) coordinate-wise, we have 

(8.59) w'^ = MjWj + dQjW, 

where \Qj\ < C, from which we obtain the growth/decay estimates 

\wj\' = Re {wj/\wj\,w'^) 

(8.60) = («'j/|w'il,Re {Mj)wj) + Re {wj/\wj\,Sejw) 

< ^rj\wj\ + Cd\w\. 
Defining r := -^^, we thus find, after some simplification, that 



W2 

, |'W2r|w'l| - |'W2||'Wl|' 
T = 

\W2\ 

im) -2^W2\\wi\ +C5\w\'^ 

I 19 

P2 



= -2^r + C5(l-fr2). 

It is thus clear that K_ := {w : r < ^} is invariant, provided 5 is small enough 
that ^ < 1. More generally, we have, for ^ < r < ^^, that r' < —rjr. Thus, K_ 
is exponentially attracting on J_ := {ty : r < 4^}-i 

Corollary 8.10. For C sufficiently large, 5/r] sufficiently small, solutions in 
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increase as x ^> +00 at exponential rate e"^, for any a < lim inf^^_|_oo a (a as 
defined in (8.50)). 

Proof. By Proposition 8.8, solutions in J_(x) eventually enter cone 

^- ■= {"^ ■ ^rrr ^ ^'^Z^^- 

\Q[x)w\ 
Thus, (8.60), translated back to original coordinates, gives 

/gg2) \W2\' >a\w2\-{C'^5'^/ri)\w2\ 

> a\w2\: 

for any 5< a, a <a —{C'^5'^/rf), from which the result follows.| 

Corollary 8.11 (The Tracking Lemma). For bjr] sufficiently small, solutions 
w~^/w~ of (8.47) decaying at +00/— 00 at rate e-^/e"^ lie always within cones 

(8.63) K_W :={..: »<CV,}. 



(8.64) K+(x) := {«, : M^^ < OS/v}. 

respectively, for any a < lini inf ^^^+00 a, a > linisup^,^.^ a. Moreover, there hold 
the uniform decay/growth rates 



(8.65) fJfW 

F (y)l 

for all X > y, and symmetrically for x < y, for any a > max^^a, a < mmj^a. 
(Note: C depends in part upon the choice of a and a). 

Proof. Without loss of generality, consider the +00 case. Bounds (8.63)-(8.64) 
follow immediately by contradiction from the previous Corollary, whence (8.65) 
follows from (8.62).| 

Remark 8.12. In the "standard" case that a < < a, a and a constant, the 
statement of Corollary 8.11 considerably simplifies; in particular, the conclusion 
applies to solutions merely decaying at +00/— 00, with no specified rate. In the case 
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of varying a and a, the rates given in (8.65) can clearly be sharpened to e-'y —^^' 

r X — / \ I 

and e-'y "^^^ ^, without any change in the argument. 

A3. 2. 2. Extensions, (i) The case r] ^ 0. The treatment above is essentially 
that given in [ZH,Z.1,Z.4], and suffices for the sectorial operators studied there. It 
suffices also for the nonsectorial, dispersive-diffusive operators studied in [HZ. 2], 
or even purely dispersive operators, since in either case a uniform spectral gap is 
maintained up to the essential spectrum boundary. However, it does not directly 
apply to the nonsectorial operators arising in the study of relaxation or real viscosity 
models, for which a, a, and rj typically vary with 5, with ry, d/rj -^ as 5 ^ 0, and 
derivative A' is typically of order rj » d. 

Here, we only point out that all conclusions of the above subsection remain valid 
also in this more general setting, provided that we can find a coordinate change 
reducing (8.47) to form (8.57)-(8.58). Such a reduction may in fact be carried out 
under rather general circumstances, as discussed in Section 4, above: for example, 
if the distance between all eigenvalues of A remains uniformly bounded from below 
even as the difference between their real parts goes to zero (as occurs for strictly 
hyperbolic relaxation models), or if the distance between groups {cti, . . .ai}{x, d) 
and {ct^-i-i, . . . aN}{x, S) remains uniformly bounded from below, while at the same 
time \A'\ < tj/C with C > sufficiently large, (as occurs for general relaxation 
models, away from the shock layer). 

(ii) Contraction mapping formulation. The analysis of this subsection can equiv- 
alently be phrased as an integral equation/contraction mapping argument giving 
existence of stable/unstable manifolds at the same time that estimates are ob- 
tained. This has the advantage of slightly further generality: for example, one can 
obtain complete short-time theory/pointwise bounds (which depend only on high- 
frequency estimates) by this technique with no assumptions other than smoothness 
on the coefficients of L. Moreover, it disassociates regularity in A from smoothness 
in X, allowing high accuracy in A with minimal regularity assumptions on coeffi- 
cients. For example, in the application to relaxation systems in Section 4, it allows 
us to terminate our diagonalization process at the first order, with 5 ~ 0(|A|~^), 
Tj ~ |A|~^, while still achieving the necessary regularity in A to order A~^, with 
error estimates to order |A|~^, whereas with the invariant cone approach we should 
have had to diagonalize to one further order, 5 ~ (9(|A|~"^, to obtain the needed 
estimates. This would have increased our smoothness requirements in hypothesis 
(HO) to /, (7, Q G C"^ for general, strictly hyperbolic relaxation systems, C^ for 
strictly hyperbolic discrete kinetic systems. 

Specifically, after first reducing to form (8.57)-(8.58), where Mj are assumed 
further to be diagonal, we may consider the "lifted equations" to the space of l- 
(resp. (A^ — t)-) forms, as in (8.44), Section A3.1, in this case appearing as 

(8.66) f]' = K^^\x;\)r] + 5er], 

where £i := ei A ■ ■ ■ Ae^ is the unique eigenvector of A*^^^ corresponding to the eigen- 
value of minimum real part, which moreover is separated from all other eigenvalues 
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with spectral gap 2r], Cj as usual denoting the standard jth basis element in C". 

Enumerating the standard basis elements of the space of £- forms as Sj, and 
defining r]j to be the coordinate representation of t] with respect to this basis, we 
obtain from structure (8.57)-(8.58) equations for rjj with the same, block-diagonal 
structure, where the first block, A^i, corresponding to rji is scalar, with Re A^i < 
— ?7, and (by a straightforward calculation using definition (8.45)) the second block 
again satisfies Re AI2 > V- Thus, defining fij := r]j/rji for j > 2, fj := {fj2, . . .), we 
obtain projectivized equations of form 

(8.67) fi={M2-MiI){x)fj + 5Q{f],x) 
where 

(8.68) Re {M2 - Mil) < -2r] 

for all j > 2 and Q(r, x) is a quadratic polynomial in r, with uniformly bounded 
coefficients, similarly as in (8.61). 

Define, similarly as in Section A3.1, the map 

POO 

(8.69) rfi{x):=6 J^'^''Q{f]{z),z)dz, 



X 



where T^^^ denotes the fiow from ^ to x of the approximate equation 

(8.70) f]={^M2-MxI){x)f\. 
From (8.68), we obtain the bound 

(8.71) IJ^^^^I < Ce-2^(^-^) 

for all z > X, similarly as in (8.60). From this, together with the quadratic form 
of Q, we readily find that T is a contraction on any ball in L°°, with constant 
C5/ri, provided that 5/ri is sufficiently small. From this observation, we obtain 
both existence of the desired fixed point ?) = T{ff) and the bound 

|r7-T(0)| = |T(r))-T(0)|<(C5/r7)|r}|, 

yielding \t)\ < C2|T(0)| < Csd/i]: the same rate of tracking as obtained before using 
the invariant cone approach. On the other hand, we have at the same time a means 
to generate the solution fj to any desired order in d/rj, by fixed point iteration: in 
particular, we have the first-order expansion 

(8.72) fi = T{0) + O{\dM'), 

where the first order term is explicitly given as an (finite) exponential integral 
against the constant term of Q. 
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Remark 8.13. Note that block-diagonal structure of A is necessary in order that 
the approximate flow (8.70) be linear; in general, it could be an arbitrary quadratic 
function, in which case Duhamel's principle would not apply. 

Remark 8.14. The estimate (8.72) in terms of £- (resp. (A^ — £)) forms yields 
a corresponding estimate in terms of subspaces: that is, a description of the stable 
manifold at +00 (resp. unstable manifold at —00) as a graph over W2 = 0, smooth 
in X and analytic in 6/r]. 

Remark 8.15. Clearly, the global decay bound (8.71) suffices for the above 
construction. The pointwise decay afforded by negativity of A42 — -Mil ^s not 
strictly necessary, and may be allowed to fail by order a on intervals with total 
length of order 1/a. This observation allows us to treat strong relaxation shocks, 
for which negativity may fail by order 77 on a shock layer of order l/ry. Likewise, in 
the tracking construction of A3. 2.1 and A3.2.2(i), negativity/positivity of M1/M2 
may be allowed to fail by order a on intervals with total length of order 1/a, during 
which solutions may move to a larger cone (by a bounded multiple) , which may be 
chosen to be invariant on "good" intervals. 

A3. 2. 3. The reduced flow. Finally, suppose that we have successfully carried 
out the reduction (8.57) of (8.47) to block-diagonal form in such a way that (8.48) 
remains valid; that we have verified positivity condition (8.58) (except perhaps on 
a finite collection of intervals of total length order l/a, where positivity may be 
allowed to fail to order a, see Remark 8.15), or by some other means established 
uniform exponential decay/growth 

(8.73) iJ^^n^C'e^^'^'""''', 

6* > 0, of fiows J-'j associated with the approximate (decoupled) equations 

(8.74) w'^ = MjWj] 

and that d/r] ^ as 5 — > 0, so that the entire stable/unstable manifold construction 
described in Sections A3. 2. 2 (i)-(ii) can be carried out (alternatively, the estimates 
of Sections A3.2.1-A3.2.2(i)). Then, the "reduced fiow" restricted to stable/unstable 
manifolds takes the simple form 

(8.75) w'j = MjWj + SQjjWj + {5^/r])nj, 

where Ojj, j = 1,2 are the diagonal blocks of error G, and terms TZj = 0{1) 
corresponding to off-diagonal blocks, arise through the basic estimate Iwd/I^Wjl < 

S/v, J T^J- 

Proposition 8.16. The flows (i.e., solution operator) ^J^^ of the reduced equa- 
tions (8.75) satisfy 

(8.76) fp"" = ^J^" + {5/r])£^{x, y, 5) + 0{5^ /rf)e-^^\^-y\ 
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for X ^ y, for any < 9 < 9, where 






satisfies the uniform exponential decay estimate 

(8.78) \£j{x,y,5)\ < Ci?7|x - yle"^''!"-^! < Cae"^"^!^"^!, 

is smooth in x and y and as smooth in d as is On, and ^]^^ are the approximate 
flows associated with (8.74) Cj and 0{-), depending only on 9 and the bounds (8.73), 
are of order C{9 — 9)~^ , where C is the constant of (8.73). 

Proof. Restrict to the case of the stable manifold, j = 1; the case j = 2 may be 
treated in symmetric fashion. Fixing y, define similarly as in previous subsections 
the map 



^.79) 



T^P" := J^y-- + /"Vr"(<5©ii + {d^/v)'J^i){z,d)fr"dz, 

Jy 



where jFf^^ denotes the fiow from ^ to x of the approximate equation (8.74). 
Then, from the decay estimate (8.73), it is readily established as in the argument 
of Section A3.2.2(ii) that T is a contraction on L°^[y, +cxd), with constant of order 
(5/?7, establishing the existence of a unique fixed-point solution. Carrying out two 
iterations of the fixed point iteration, we obtain the estimate 

(880) = ^r"" + f ^r-" (<^en + (5Vr7)7^l) (^, 5)^r^ dz + 0{5Vrf) 

= JT" + ^ r ^r^"©ii {z, s)j'r" dz + o(5 v^') , 

yielding (8.76) with ^ = and (8.77). 
The term Si may be estimated as 



PX PX 

\£i\<V \T^^''\\eii{z,6)\\J^f^'\dz<Cri e-^v\x-z\^-0v\y-z\^^ 

J V J V 



e 
y -^y 

sv\x-y\ ^ nr,(\^ — '„\p-{^-^)'n\^-y\\p-^'n\^-y\ 



= Cri\x - yle-'^"!"-^! < Cri{\x - y\e-^"-">'^\''-y\)e 
<C{9-9)-^e-^'^\''-y\ 

for any < 9 < 9, and likewise the term involving T^i in the second line of (8.80) 
may be estimated as (9(5^/?7^)e~^'''^~^'). By a similar calculation, we can show 
that T is actually a contraction in the exponentially weighted norm 



,~:=||/(y''(-^)| 



L°°[y, +00)7 
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yielding the improved bound 0{5'^ /rf)e~^'^^^~'^\ on the third, 0{-) term in the 
second hne of (8.80). SpecificaUy, we have 



2 He 



< e^'^'i^-^iii^i - T2\\e8 r I J-riKeii + (V^)7^l)(^, 5)1 e-'>\'-y\ dz 

Jy 

= C5||J^i - J^2|le / e-(^-^')^l"-"l dz 
Jy 

Combining these observations, we obtain (8.76) as claimed. | 

Remark 8.17. Represention (8.76) becomes particularly simple in the case that 
Mj are diagonal, and ^J^^ = e-'y ^^^'dz. This occurs, for example, in the case 
of strictly hyperbolic relaxation systems, or in the case of scalar (third and) higher 
order equations considered in [HZ. 2]. However, for more general applications, it 
is important that we do not limit ourselves to this case. For nonstrictly hyperbolic 
relaxation systems, and for higher order systems, Mj are block diagonal, with blocks 
of form 

mj,k = (3j,k + 0{5), 

f3j^k scalar. Thus, it is possible again to express the dominant (typically highly 
oscillatory) effects in each block as a simple exponential integral factor multiplying 
a higher-order, coupled flow. We describe such a calculation in detail in Section 4, 
above. 

Note that we do not require in Proposition 8.16 any spectral gap between different 
modes (blocks) of Mj in order to obtain an approximately block-diagonal flow. This 
is important, for example, in the case of dispersive-diffusive equations considered 
in [HZ. 2]. Indeed, Proposition 8.16 and the contraction mapping construction of 
Section A3.2.2(ii) together repair an omission in the analogous constructions carried 
out in [HZ. 2]: speciflcally in the proof of Lemma 3.2, [HZ. 2], wherein a nonexistent 
gap was implicitly assumed.^ 
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